Probabilistic boundary element method

NASA Technical Reports Server (NTRS)

Cruse, T. A.; Raveendra, S. T.

1989-01-01

The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.

Developments in boundary element methods - 2

NASA Astrophysics Data System (ADS)

Banerjee, P. K.; Shaw, R. P.

This book is a continuation of the effort to demonstrate the power and versatility of boundary element methods which began in Volume 1 of this series. While Volume 1 was designed to introduce the reader to a selected range of problems in engineering for which the method has been shown to be efficient, the present volume has been restricted to time-dependent problems in engineering. Boundary element formulation for melting and solidification problems in considered along with transient flow through porous elastic media, applications of boundary element methods to problems of water waves, and problems of general viscous flow. Attention is given to time-dependent inelastic deformation of metals by boundary element methods, the determination of eigenvalues by boundary element methods, transient stress analysis of tunnels and caverns of arbitrary shape due to traveling waves, an analysis of hydrodynamic loads by boundary element methods, and acoustic emissions from submerged structures.

Compact Fuel Element Environment Test

NASA Technical Reports Server (NTRS)

Bradley, D. E.; Mireles, O. R.; Hickman, R. R.; Broadway, J. W.

2012-01-01

Deep space missions with large payloads require high specific impulse (I(sub sp)) and relatively high thrust to achieve mission goals in reasonable time frames. Conventional, storable propellants produce average I(sub sp). Nuclear thermal rockets (NTRs) capable of high I(sub sp) thrust have been proposed. NTR employs heat produced by fission reaction to heat and therefore accelerate hydrogen, which is then forced through a rocket nozzle providing thrust. Fuel element temperatures are very high (up to 3,000 K) and hydrogen is highly reactive with most materials at high temperatures. Data covering the effects of high-temperature hydrogen exposure on fuel elements are limited. The primary concern is the mechanical failure of fuel elements that employ high melting point metals, ceramics, or a combination (cermet) as a structural matrix into which the nuclear fuel is distributed. It is not necessary to include fissile material in test samples intended to explore high-temperature hydrogen exposure of the structural support matrices. A small-scale test bed designed to heat fuel element samples via noncontact radio frequency heating and expose samples to hydrogen for typical mission durations has been developed to assist in optimal material and manufacturing process selection without employing fissile material. This Technical Memorandum details the test bed design and results of testing conducted to date.

Discontinuous Finite Element Quasidiffusion Methods

DOE PAGES

Anistratov, Dmitriy Yurievich; Warsa, James S.

2018-05-21

Here in this paper, two-level methods for solving transport problems in one-dimensional slab geometry based on the quasi-diffusion (QD) method are developed. A linear discontinuous finite element method (LDFEM) is derived for the spatial discretization of the low-order QD (LOQD) equations. It involves special interface conditions at the cell edges based on the idea of QD boundary conditions (BCs). We consider different kinds of QD BCs to formulate the necessary cell-interface conditions. We develop two-level methods with independent discretization of the high-order transport equation and LOQD equations, where the transport equation is discretized using the method of characteristics and themore » LDFEM is applied to the LOQD equations. We also formulate closures that lead to the discretization consistent with a LDFEM discretization of the transport equation. The proposed methods are studied by means of test problems formulated with the method of manufactured solutions. Numerical experiments are presented demonstrating the performance of the proposed methods. Lastly, we also show that the method with independent discretization has the asymptotic diffusion limit.« less

Discontinuous Finite Element Quasidiffusion Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Anistratov, Dmitriy Yurievich; Warsa, James S.

Here in this paper, two-level methods for solving transport problems in one-dimensional slab geometry based on the quasi-diffusion (QD) method are developed. A linear discontinuous finite element method (LDFEM) is derived for the spatial discretization of the low-order QD (LOQD) equations. It involves special interface conditions at the cell edges based on the idea of QD boundary conditions (BCs). We consider different kinds of QD BCs to formulate the necessary cell-interface conditions. We develop two-level methods with independent discretization of the high-order transport equation and LOQD equations, where the transport equation is discretized using the method of characteristics and themore » LDFEM is applied to the LOQD equations. We also formulate closures that lead to the discretization consistent with a LDFEM discretization of the transport equation. The proposed methods are studied by means of test problems formulated with the method of manufactured solutions. Numerical experiments are presented demonstrating the performance of the proposed methods. Lastly, we also show that the method with independent discretization has the asymptotic diffusion limit.« less

Evaluation of Test Methods for Triaxially Braided Composites using a Meso-Scale Finite Element Model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Chao

The characterization of triaxially braided composite is complicate due to the nonuniformity of deformation within the unit cell as well as the possibility of the freeedge effect related to the large size of the unit cell. Extensive experimental investigation has been conducted to develop more accurate test approaches in characterizing the actual mechanical properties of the material we are studying. In this work, a meso-scale finite element model is utilized to simulate two complex specimens: notched tensile specimen and tube tensile specimen, which are designed to avoid the free-edge effect and free-edge effect induced premature edge damage. The full fieldmore » strain data is predicted numerically and compared with experimental data obtained by Digit Image Correlation. The numerically predicted tensile strength values are compared with experimentally measured results. The discrepancy between numerically predicted and experimentally measured data, the capability of different test approaches are analyzed and discussed. The presented numerical model could serve as assistance to the evaluation of different test methods, and is especially useful in identifying potential local damage events.« less

REACTOR FUEL ELEMENTS TESTING CONTAINER

DOEpatents

Whitham, G.K.; Smith, R.R.

1963-01-15

This patent shows a method for detecting leaks in jacketed fuel elements. The element is placed in a sealed tank within a nuclear reactor, and, while the reactor operates, the element is sparged with gas. The gas is then led outside the reactor and monitored for radioactive Xe or Kr. (AEC)

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

NASA Astrophysics Data System (ADS)

Yang, Zailin; Wang, Yao; Hei, Baoping

2013-12-01

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of onedimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.

Multi-element microelectropolishing method

DOEpatents

Lee, Peter J.

1994-01-01

A method is provided for microelectropolishing a transmission electron microscopy nonhomogeneous multi-element compound foil. The foil is electrolyzed at different polishing rates for different elements by rapidly cycling between different current densities. During a first portion of each cycle at a first voltage a first element electrolyzes at a higher current density than a second element such that the material of the first element leaves the anode foil at a faster rate than the second element and creates a solid surface film, and such that the solid surface film is removed at a faster rate than the first element leaves the anode foil. During a second portion of each cycle at a second voltage the second element electrolyzes at a higher current density than the first element, and the material of the second element leaves the anode foil at a faster rate than the first element and creates a solid surface film, and the solid surface film is removed at a slower rate than the second element leaves the foil. The solid surface film is built up during the second portion of the cycle, and removed during the first portion of the cycle.

Multi-element microelectropolishing method

DOEpatents

Lee, P.J.

1994-10-11

A method is provided for microelectropolishing a transmission electron microscopy nonhomogeneous multi-element compound foil. The foil is electrolyzed at different polishing rates for different elements by rapidly cycling between different current densities. During a first portion of each cycle at a first voltage a first element electrolyzes at a higher current density than a second element such that the material of the first element leaves the anode foil at a faster rate than the second element and creates a solid surface film, and such that the solid surface film is removed at a faster rate than the first element leaves the anode foil. During a second portion of each cycle at a second voltage the second element electrolyzes at a higher current density than the first element, and the material of the second element leaves the anode foil at a faster rate than the first element and creates a solid surface film, and the solid surface film is removed at a slower rate than the second element leaves the foil. The solid surface film is built up during the second portion of the cycle, and removed during the first portion of the cycle. 10 figs.

Coupling Finite Element and Meshless Local Petrov-Galerkin Methods for Two-Dimensional Potential Problems

NASA Technical Reports Server (NTRS)

Chen, T.; Raju, I. S.

2002-01-01

A coupled finite element (FE) method and meshless local Petrov-Galerkin (MLPG) method for analyzing two-dimensional potential problems is presented in this paper. The analysis domain is subdivided into two regions, a finite element (FE) region and a meshless (MM) region. A single weighted residual form is written for the entire domain. Independent trial and test functions are assumed in the FE and MM regions. A transition region is created between the two regions. The transition region blends the trial and test functions of the FE and MM regions. The trial function blending is achieved using a technique similar to the 'Coons patch' method that is widely used in computer-aided geometric design. The test function blending is achieved by using either FE or MM test functions on the nodes in the transition element. The technique was evaluated by applying the coupled method to two potential problems governed by the Poisson equation. The coupled method passed all the patch test problems and gave accurate solutions for the problems studied.

Discrete Element Method and its application to materials failure problem on the example of Brazilian Test

NASA Astrophysics Data System (ADS)

Klejment, Piotr; Kosmala, Alicja; Foltyn, Natalia; Dębski, Wojciech

2017-04-01

The earthquake focus is the point where a rock under external stress starts to fracture. Understanding earthquake nucleation and earthquake dynamics requires thus understanding of fracturing of brittle materials. This, however, is a continuing problem and enduring challenge to geoscience. In spite of significant progress we still do not fully understand the failure of rock materials due to extreme stress concentration in natural condition. One of the reason of this situation is that information about natural or induced seismic events is still not sufficient for precise description of physical processes in seismic foci. One of the possibility of improving this situation is using numerical simulations - a powerful tool of contemporary physics. For this reason we used an advanced implementation of the Discrete Element Method (DEM). DEM's main task is to calculate physical properties of materials which are represented as an assembly of a great number of particles interacting with each other. We analyze the possibility of using DEM for describing materials during so called Brazilian Test. Brazilian Test is a testing method to obtain the tensile strength of brittle material. One of the primary reasons for conducting such simulations is to measure macroscopic parameters of the rock sample. We would like to report our efforts of describing the fracturing process during the Brazilian Test from the microscopic point of view and give an insight into physical processes preceding materials failure.

Apparatus for testing high pressure injector elements

NASA Technical Reports Server (NTRS)

Myers, William Neill (Inventor); Scott, Ewell M. (Inventor); Forbes, John C. (Inventor); Shadoan, Michael D. (Inventor)

1995-01-01

An apparatus for testing and evaluating the spray pattern of high pressure fuel injector elements for use in supplying fuel to combustion engines is presented. Prior art fuel injector elements were normally tested by use of low pressure apparatuses which did not provide a purge to prevent mist from obscuring the injector element or to prevent frosting of the view windows; could utilize only one fluid during each test; and had their viewing ports positioned one hundred eighty (180 deg) apart, thus preventing optimum use of laser diagnostics. The high pressure fluid injector test apparatus includes an upper hub, an upper weldment or housing, a first clamp and stud/nut assembly for securing the upper hub to the upper weldment, a standoff assembly within the upper weldment, a pair of window housings having view glasses within the upper weldment, an injector block assembly and purge plate within the upper weldment for holding an injector element to be tested and evaluated, a lower weldment or housing, a second clamp and stud/nut assembly for securing the lower weldment to the upper hub, a third clamp and stud/nut assembly for securing the lower hub to the lower weldment, mechanisms for introducing fluid under high pressure for testing an injector element, and mechanisms for purging the apparatus to prevent frosting of view glasses within the window housings and to permit unobstructed viewing of the injector element.

Apparatus for testing high pressure injector elements

NASA Technical Reports Server (NTRS)

Myers, William Neill (Inventor); Scott, Ewell M. (Inventor); Forbes, John C. (Inventor); Shadoan, Michael D. (Inventor)

1993-01-01

An apparatus for testing and evaluating the spray pattern of high pressure fuel injector elements for use in supplying fuel to combustion engines is presented. Prior art fuel injector elements were normally tested by use of low pressure apparatuses which did not provide a purge to prevent mist from obscuring the injector element or to prevent frosting of the view windows; could utilize only one fluid during each test; and had their viewing ports positioned one hundred eighty (180 deg) apart, thus preventing optimum use of laser diagnostics. The high pressure fluid injector test apparatus includes an upper hub, an upper weldment or housing, a first clamp and stud/nut assembly for securing the upper hub to the upper weldment, a standoff assembly within the upper weldment, a pair of window housings having view glasses within the upper weldment, an injector block assembly and purge plate within the upper weldment for holding an injector element to be tested and evaluated, a lower weldment or housing, a second clamp and stud/nut assembly for securing the lower weldment to the upper weldment, a lower hub, a third clamp and stud/nut assembly for securing the lower hub to the lower weldment, mechanisms for introducing fluid under high pressure for testing an injector element, and mechanisms for purging the apparatus to prevent frosting of view glasses within the window housings and to permit unobstructed viewing of the injector element.

COMPLEX VARIABLE BOUNDARY ELEMENT METHOD: APPLICATIONS.

USGS Publications Warehouse

Hromadka, T.V.; Yen, C.C.; Guymon, G.L.

1985-01-01

The complex variable boundary element method (CVBEM) is used to approximate several potential problems where analytical solutions are known. A modeling result produced from the CVBEM is a measure of relative error in matching the known boundary condition values of the problem. A CVBEM error-reduction algorithm is used to reduce the relative error of the approximation by adding nodal points in boundary regions where error is large. From the test problems, overall error is reduced significantly by utilizing the adaptive integration algorithm.

Acceleration of boundary element method for linear elasticity

NASA Astrophysics Data System (ADS)

Zapletal, Jan; Merta, Michal; Čermák, Martin

2017-07-01

In this work we describe the accelerated assembly of system matrices for the boundary element method using the Intel Xeon Phi coprocessors. We present a model problem, provide a brief overview of its discretization and acceleration of the system matrices assembly using the coprocessors, and test the accelerated version using a numerical benchmark.

Application of the control volume mixed finite element method to a triangular discretization

USGS Publications Warehouse

Naff, R.L.

2012-01-01

A two-dimensional control volume mixed finite element method is applied to the elliptic equation. Discretization of the computational domain is based in triangular elements. Shape functions and test functions are formulated on the basis of an equilateral reference triangle with unit edges. A pressure support based on the linear interpolation of elemental edge pressures is used in this formulation. Comparisons are made between results from the standard mixed finite element method and this control volume mixed finite element method. Published 2011. This article is a US Government work and is in the public domain in the USA. ?? 2012 John Wiley & Sons, Ltd. This article is a US Government work and is in the public domain in the USA.

Use of system identification techniques for improving airframe finite element models using test data

NASA Technical Reports Server (NTRS)

Hanagud, Sathya V.; Zhou, Weiyu; Craig, James I.; Weston, Neil J.

1991-01-01

A method for using system identification techniques to improve airframe finite element models was developed and demonstrated. The method uses linear sensitivity matrices to relate changes in selected physical parameters to changes in total system matrices. The values for these physical parameters were determined using constrained optimization with singular value decomposition. The method was confirmed using both simple and complex finite element models for which pseudo-experimental data was synthesized directly from the finite element model. The method was then applied to a real airframe model which incorporated all the complexities and details of a large finite element model and for which extensive test data was available. The method was shown to work, and the differences between the identified model and the measured results were considered satisfactory.

Design of horizontal-axis wind turbine using blade element momentum method

NASA Astrophysics Data System (ADS)

Bobonea, Andreea; Pricop, Mihai Victor

2013-10-01

The study of mathematical models applied to wind turbine design in recent years, principally in electrical energy generation, has become significant due to the increasing use of renewable energy sources with low environmental impact. Thus, this paper shows an alternative mathematical scheme for the wind turbine design, based on the Blade Element Momentum (BEM) Theory. The results from the BEM method are greatly dependent on the precision of the lift and drag coefficients. The basic of BEM method assumes the blade can be analyzed as a number of independent element in spanwise direction. The induced velocity at each element is determined by performing the momentum balance for a control volume containing the blade element. The aerodynamic forces on the element are calculated using the lift and drag coefficient from the empirical two-dimensional wind tunnel test data at the geometric angle of attack (AOA) of the blade element relative to the local flow velocity.

Discontinuous finite element method for vector radiative transfer

NASA Astrophysics Data System (ADS)

Wang, Cun-Hai; Yi, Hong-Liang; Tan, He-Ping

2017-03-01

The discontinuous finite element method (DFEM) is applied to solve the vector radiative transfer in participating media. The derivation in a discrete form of the vector radiation governing equations is presented, in which the angular space is discretized by the discrete-ordinates approach with a local refined modification, and the spatial domain is discretized into finite non-overlapped discontinuous elements. The elements in the whole solution domain are connected by modelling the boundary numerical flux between adjacent elements, which makes the DFEM numerically stable for solving radiative transfer equations. Several various problems of vector radiative transfer are tested to verify the performance of the developed DFEM, including vector radiative transfer in a one-dimensional parallel slab containing a Mie/Rayleigh/strong forward scattering medium and a two-dimensional square medium. The fact that DFEM results agree very well with the benchmark solutions in published references shows that the developed DFEM in this paper is accurate and effective for solving vector radiative transfer problems.

Three-dimensional finite element analysis of the shear bond test.

PubMed

DeHoff, P H; Anusavice, K J; Wang, Z

1995-03-01

The purpose of this study was to use finite element analyses to model the planar shear bond test and to evaluate the effects of modulus values, bonding agent thickness, and loading conditions on the stress distribution in the dentin adjacent to the bonding agent-dentin interface. All calculations were performed with the ANSYS finite element program. The planar shear bond test was modeled as a cylinder of resin-based composite bonded to a cylindrical dentin substrate. The effects of material, geometry and loading variables were determined primarily by use of a three-dimensional structural element. Several runs were also made using an axisymmetric element with harmonic loading and a plane strain element to determine whether two-dimensional analyses yield valid results. Stress calculations using three-dimensional finite element analyses confirmed the presence of large stress concentration effects for all stress components at the bonding agent-dentin interface near the application of the load. The maximum vertical shear stress generally occurs approximately 0.3 mm below the loading site and then decreases sharply in all directions. The stresses reach relatively uniform conditions within about 0.5 mm of the loading site and then increase again as the lower region of the interface is approached. Calculations using various loading conditions indicated that a wire-loop method of loading leads to smaller stress concentration effects, but a shear bond strength determined by dividing a failure load by the cross-sectional area grossly underestimates the true interfacial bond strength. Most dental researchers are using tensile and shear bond tests to predict the effects of process and material variables on the clinical performance of bonding systems but no evidence has yet shown that bond strength is relevant to clinical performance. A critical factor in assessing the usefulness of bond tests is a thorough understanding of the stress states that cause failure in the bond test

Space reactor fuel element testing in upgraded TREAT

NASA Astrophysics Data System (ADS)

Todosow, M.; Bezler, P.; Ludewig, H.; Kato, W. Y.

The testing of candidate fuel elements at prototypic operating conditions with respect to temperature, power density, hydrogen coolant flow rate, etc.; a crucial component in the development and qualification of nuclear rocket engines based on the Particle Bed Reactor (PBR); NERVA-derivative; and other concepts are discussed. Such testing may be performed at existing reactors, or at new facilities. A scoping study has been performed to assess the feasibility of testing PBR based fuel elements at the TREAT reactor. Initial results suggest that full-scale PBR elements could be tested at an average energy deposition of approximately 60-80 MW-s/L in the current TREAT reactor. If the TREAT reactor was upgraded to include fuel elements with a higher temperature limit, average energy deposition of approximately 100 MW/L may be achievable.

Space reactor fuel element testing in upgraded TREAT

NASA Astrophysics Data System (ADS)

Todosow, Michael; Bezler, Paul; Ludewig, Hans; Kato, Walter Y.

1993-01-01

The testing of candidate fuel elements at prototypic operating conditions with respect to temperature, power density, hydrogen coolant flow rate, etc., is a crucial component in the development and qualification of nuclear rocket engines based on the Particle Bed Reactor (PBR), NERVA-derivative, and other concepts. Such testing may be performed at existing reactors, or at new facilities. A scoping study has been performed to assess the feasibility of testing PBR based fuel elements at the TREAT reactor. Initial results suggests that full-scale PBR elements could be tested at an average energy deposition of ˜60-80 MW-s/L in the current TREAT reactor. If the TREAT reactor was upgraded to include fuel elements with a higher temperture limit, average energy deposition of ˜100 MW/L may be achievable.

Safety envelope for load tolerance of structural element design based on multi-stage testing

DOE PAGES

Park, Chanyoung; Kim, Nam H.

2016-09-06

Structural elements, such as stiffened panels and lap joints, are basic components of aircraft structures. For aircraft structural design, designers select predesigned elements satisfying the design load requirement based on their load-carrying capabilities. Therefore, estimation of safety envelope of structural elements for load tolerances would be a good investment for design purpose. In this article, a method of estimating safety envelope is presented using probabilistic classification, which can estimate a specific level of failure probability under both aleatory and epistemic uncertainties. An important contribution of this article is that the calculation uncertainty is reflected in building a safety envelope usingmore » Gaussian process, and the effect of element test data on reducing the calculation uncertainty is incorporated by updating the Gaussian process model with the element test data. It is shown that even one element test can significantly reduce the calculation uncertainty due to lacking knowledge of actual physics, so that conservativeness in a safety envelope is significantly reduced. The proposed approach was demonstrated with a cantilever beam example, which represents a structural element. The example shows that calculation uncertainty provides about 93% conservativeness against the uncertainty due to a few element tests. As a result, it is shown that even a single element test can increase the load tolerance modeled with the safety envelope by 20%.« less

Space reactor fuel element testing in upgraded TREAT

DOE Office of Scientific and Technical Information (OSTI.GOV)

Todosow, M.; Bezler, P.; Ludewig, H.

1993-01-14

The testing of candidate fuel elements at prototypic operating conditions with respect to temperature, power density, hydrogen coolant flow rate, etc., a crucial component in the development and qualification of nuclear rocket engines based on the Particle Bed Reactor (PBR), NERVA-derivative, and other concepts. Such testing may be performed at existing reactors, or at new facilities. A scoping study has been performed to assess the feasibility of testing PBR based fuel elements at the TREAT reactor. initial results suggest that full-scale PBR, elements could be tested at an average energy deposition of {approximately}60--80 MW-s/L in the current TREAT reactor. Ifmore » the TREAT reactor was upgraded to include fuel elements with a higher temperature limit, average energy deposition of {approximately}100 MW/L may be achievable.« less

Space reactor fuel element testing in upgraded TREAT

DOE Office of Scientific and Technical Information (OSTI.GOV)

Todosow, M.; Bezler, P.; Ludewig, H.

1993-05-01

The testing of candidate fuel elements at prototypic operating conditions with respect to temperature, power density, hydrogen coolant flow rate, etc., a crucial component in the development and qualification of nuclear rocket engines based on the Particle Bed Reactor (PBR), NERVA-derivative, and other concepts. Such testing may be performed at existing reactors, or at new facilities. A scoping study has been performed to assess the feasibility of testing PBR based fuel elements at the TREAT reactor. initial results suggest that full-scale PBR, elements could be tested at an average energy deposition of {approximately}60--80 MW-s/L in the current TREAT reactor. Ifmore » the TREAT reactor was upgraded to include fuel elements with a higher temperature limit, average energy deposition of {approximately}100 MW/L may be achievable.« less

Fast Boundary Element Method for acoustics with the Sparse Cardinal Sine Decomposition

NASA Astrophysics Data System (ADS)

Alouges, François; Aussal, Matthieu; Parolin, Emile

2017-07-01

This paper presents the newly proposed method Sparse Cardinal Sine Decomposition that allows fast convolution on unstructured grids. We focus on its use when coupled with finite element techniques to solve acoustic problems with the (compressed) Boundary Element Method. In addition, we also compare the computational performances of two equivalent Matlab® and Python implementations of the method. We show validation test cases in order to assess the precision of the approach. Eventually, the performance of the method is illustrated by the computation of the acoustic target strength of a realistic submarine from the Benchmark Target Strength Simulation international workshop.

Peridynamic Multiscale Finite Element Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Costa, Timothy; Bond, Stephen D.; Littlewood, David John

The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic andmore » local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and

Weak Galerkin finite element methods for Darcy flow: Anisotropy and heterogeneity

NASA Astrophysics Data System (ADS)

Lin, Guang; Liu, Jiangguo; Mu, Lin; Ye, Xiu

2014-11-01

This paper presents a family of weak Galerkin finite element methods (WGFEMs) for Darcy flow computation. The WGFEMs are new numerical methods that rely on the novel concept of discrete weak gradients. The WGFEMs solve for pressure unknowns both in element interiors and on the mesh skeleton. The numerical velocity is then obtained from the discrete weak gradient of the numerical pressure. The new methods are quite different than many existing numerical methods in that they are locally conservative by design, the resulting discrete linear systems are symmetric and positive-definite, and there is no need for tuning problem-dependent penalty factors. We test the WGFEMs on benchmark problems to demonstrate the strong potential of these new methods in handling strong anisotropy and heterogeneity in Darcy flow.

Adaptive Finite Element Methods for Continuum Damage Modeling

NASA Technical Reports Server (NTRS)

Min, J. B.; Tworzydlo, W. W.; Xiques, K. E.

1995-01-01

The paper presents an application of adaptive finite element methods to the modeling of low-cycle continuum damage and life prediction of high-temperature components. The major objective is to provide automated and accurate modeling of damaged zones through adaptive mesh refinement and adaptive time-stepping methods. The damage modeling methodology is implemented in an usual way by embedding damage evolution in the transient nonlinear solution of elasto-viscoplastic deformation problems. This nonlinear boundary-value problem is discretized by adaptive finite element methods. The automated h-adaptive mesh refinements are driven by error indicators, based on selected principal variables in the problem (stresses, non-elastic strains, damage, etc.). In the time domain, adaptive time-stepping is used, combined with a predictor-corrector time marching algorithm. The time selection is controlled by required time accuracy. In order to take into account strong temperature dependency of material parameters, the nonlinear structural solution a coupled with thermal analyses (one-way coupling). Several test examples illustrate the importance and benefits of adaptive mesh refinements in accurate prediction of damage levels and failure time.

Single element injector testing for STME injector technology

NASA Technical Reports Server (NTRS)

Hulka, J.; Schneider, J. A.; Davis, J.

1992-01-01

An oxidizer-swirled coaxial element injector is being developed for application in the liquid oxygen/gaseous hydrogen Space Transportation Main Engine (STME) for the National Launch System (NLS) vehicle. This paper reports on the first two parts of a four part single injector element study for optimization of the STME injector design. Measurements of Rupe mixing efficiency and atomization characteristics are reported for single element versions of injection elements from two multielement injectors that have been recently hot fire tested. Rather than attempting to measure a definitive mixing efficiency or droplet size parameters of these injector elements, the purpose of these experiments was to provide a baseline comparison for evaluating future injector element design modifications. Hence, all the experiments reported here were conducted with cold flow simulants to nonflowing, ambient conditions. Mixing experiments were conducted with liquid/liquid simulants to provide economical trend data. Atomization experiments were conducted with liquid/gas simulants without backpressure. The results, despite significant differences from hot fire conditions, were found to relate to mixing and atomization parameters deduced from the hot fire testing, suggesting that these experiments are valid for trend analyses. Single element and subscale multielement hot fire testing will verify optimized designs before committing to fullscale fabrication.

[Standard sample preparation method for quick determination of trace elements in plastic].

PubMed

Yao, Wen-Qing; Zong, Rui-Long; Zhu, Yong-Fa

2011-08-01

Reference sample was prepared by masterbatch method, containing heavy metals with known concentration of electronic information products (plastic), the repeatability and precision were determined, and reference sample preparation procedures were established. X-Ray fluorescence spectroscopy (XRF) analysis method was used to determine the repeatability and uncertainty in the analysis of the sample of heavy metals and bromine element. The working curve and the metrical methods for the reference sample were carried out. The results showed that the use of the method in the 200-2000 mg x kg(-1) concentration range for Hg, Pb, Cr and Br elements, and in the 20-200 mg x kg(-1) range for Cd elements, exhibited a very good linear relationship, and the repeatability of analysis methods for six times is good. In testing the circuit board ICB288G and ICB288 from the Mitsubishi Heavy Industry Company, results agreed with the recommended values.

Modal Testing of Seven Shuttle Cargo Elements for Space Station

NASA Technical Reports Server (NTRS)

Kappus, Kathy O.; Driskill, Timothy C.; Parks, Russel A.; Patterson, Alan (Technical Monitor)

2001-01-01

From December 1996 to May 2001, the Modal and Control Dynamics Team at NASA's Marshall Space Flight Center (MSFC) conducted modal tests on seven large elements of the International Space Station. Each of these elements has been or will be launched as a Space Shuttle payload for transport to the International Space Station (ISS). Like other Shuttle payloads, modal testing of these elements was required for verification of the finite element models used in coupled loads analyses for launch and landing. The seven modal tests included three modules - Node, Laboratory, and Airlock, and four truss segments - P6, P3/P4, S1/P1, and P5. Each element was installed and tested in the Shuttle Payload Modal Test Bed at MSFC. This unique facility can accommodate any Shuttle cargo element for modal test qualification. Flexure assemblies were utilized at each Shuttle-to-payload interface to simulate a constrained boundary in the load carrying degrees of freedom. For each element, multiple-input, multiple-output burst random modal testing was the primary approach with controlled input sine sweeps for linearity assessments. The accelerometer channel counts ranged from 252 channels to 1251 channels. An overview of these tests, as well as some lessons learned, will be provided in this paper.

Weak Galerkin finite element methods for Darcy flow: Anisotropy and heterogeneity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lin, Guang; Liu, Jiangguo; Mu, Lin

2014-11-01

This paper presents a family of weak Galerkin finite element methods (WGFEMs) for Darcy flow computation. The WGFEMs are new numerical methods that rely on the novel concept of discrete weak gradients. The WGFEMs solve for pressure unknowns both in element interiors and on the mesh skeleton. The numerical velocity is then obtained from the discrete weak gradient of the numerical pressure. The new methods are quite different than many existing numerical methods in that they are locally conservative by design, the resulting discrete linear systems are symmetric and positive-definite, and there is no need for tuning problem-dependent penalty factors.more » We test the WGFEMs on benchmark problems to demonstrate the strong potential of these new methods in handling strong anisotropy and heterogeneity in Darcy flow.« less

Method of lightening radiation darkened optical elements

DOEpatents

Reich, Frederich R.; Schwankoff, Albert R.

1980-01-01

A method of lightening a radiation-darkened optical element in wich visible optical energy or electromagnetic radiation having a wavelength in the range of from about 2000 to about 20,000 angstroms is directed into the radiation-darkened optical element; the method may be used to lighten radiation-darkened optical element in-situ during the use of the optical element to transmit data by electronically separating the optical energy from the optical output by frequency filtering, data cooling, or interlacing the optic energy between data intervals.

Analysis of random structure-acoustic interaction problems using coupled boundary element and finite element methods

NASA Technical Reports Server (NTRS)

Mei, Chuh; Pates, Carl S., III

1994-01-01

A coupled boundary element (BEM)-finite element (FEM) approach is presented to accurately model structure-acoustic interaction systems. The boundary element method is first applied to interior, two and three-dimensional acoustic domains with complex geometry configurations. Boundary element results are very accurate when compared with limited exact solutions. Structure-interaction problems are then analyzed with the coupled FEM-BEM method, where the finite element method models the structure and the boundary element method models the interior acoustic domain. The coupled analysis is compared with exact and experimental results for a simplistic model. Composite panels are analyzed and compared with isotropic results. The coupled method is then extended for random excitation. Random excitation results are compared with uncoupled results for isotropic and composite panels.

The nonconforming virtual element method for eigenvalue problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gardini, Francesca; Manzini, Gianmarco; Vacca, Giuseppe

We analyse the nonconforming Virtual Element Method (VEM) for the approximation of elliptic eigenvalue problems. The nonconforming VEM allow to treat in the same formulation the two- and three-dimensional case.We present two possible formulations of the discrete problem, derived respectively by the nonstabilized and stabilized approximation of the L 2-inner product, and we study the convergence properties of the corresponding discrete eigenvalue problems. The proposed schemes provide a correct approximation of the spectrum and we prove optimal-order error estimates for the eigenfunctions and the usual double order of convergence of the eigenvalues. Finally we show a large set of numericalmore » tests supporting the theoretical results, including a comparison with the conforming Virtual Element choice.« less

Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions.

PubMed

Bauler, Patricia; Huber, Gary A; McCammon, J Andrew

2012-04-28

Diffusion is often the rate determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. This paper proposes a new hybrid diffusion method that couples the strengths of each of these two methods. The method is derived for a general multidimensional system, and is presented using a basic test case for 1D linear and radially symmetric diffusion systems.

Material nonlinear analysis via mixed-iterative finite element method

NASA Technical Reports Server (NTRS)

Sutjahjo, Edhi; Chamis, Christos C.

1992-01-01

The performance of elastic-plastic mixed-iterative analysis is examined through a set of convergence studies. Membrane and bending behaviors are tested using 4-node quadrilateral finite elements. The membrane result is excellent, which indicates the implementation of elastic-plastic mixed-iterative analysis is appropriate. On the other hand, further research to improve bending performance of the method seems to be warranted.

The Relation of Finite Element and Finite Difference Methods

NASA Technical Reports Server (NTRS)

Vinokur, M.

1976-01-01

Finite element and finite difference methods are examined in order to bring out their relationship. It is shown that both methods use two types of discrete representations of continuous functions. They differ in that finite difference methods emphasize the discretization of independent variable, while finite element methods emphasize the discretization of dependent variable (referred to as functional approximations). An important point is that finite element methods use global piecewise functional approximations, while finite difference methods normally use local functional approximations. A general conclusion is that finite element methods are best designed to handle complex boundaries, while finite difference methods are superior for complex equations. It is also shown that finite volume difference methods possess many of the advantages attributed to finite element methods.

Improved finite-element methods for rotorcraft structures

NASA Technical Reports Server (NTRS)

Hinnant, Howard E.

1991-01-01

An overview of the research directed at improving finite-element methods for rotorcraft airframes is presented. The development of a modification to the finite element method which eliminates interelement discontinuities is covered. The following subject areas are discussed: geometric entities, interelement continuity, dependent rotational degrees of freedom, and adaptive numerical integration. This new methodology is being implemented as an anisotropic, curvilinear, p-version, beam, shell, and brick finite element program.

The Mimetic Finite Element Method and the Virtual Element Method for elliptic problems with arbitrary regularity.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Manzini, Gianmarco

2012-07-13

We develop and analyze a new family of virtual element methods on unstructured polygonal meshes for the diffusion problem in primal form, that use arbitrarily regular discrete spaces V{sub h} {contained_in} C{sup {alpha}} {element_of} N. The degrees of freedom are (a) solution and derivative values of various degree at suitable nodes and (b) solution moments inside polygons. The convergence of the method is proven theoretically and an optimal error estimate is derived. The connection with the Mimetic Finite Difference method is also discussed. Numerical experiments confirm the convergence rate that is expected from the theory.

Loading simulation of lumbar spine vertebrae during a compression test using the finite elements method and trabecular bone strength properties, determined by means of nanoindentations.

PubMed

Bouzakis, K D; Mitsi, S; Michailidis, N; Mirisidis, I; Mesomeris, G; Maliaris, G; Korlos, A; Kapetanos, G; Antonarakos, P; Anagnostidis, K

2004-06-01

The mechanical strength properties of lumbar spine vertebrae are of great importance in a wide range of applications. Herein, through nanoindentations and appropriate evaluation of the corresponding results, trabecular bone struts stress-strain characteristics can be determined. In the frame of the present paper, an L2 fresh cadaveric vertebra, from which posterior elements were removed, was subjected to compression. With the aid of developed finite elements method based algorithms, the cortical shell and the cancellous core bulk elasticity moduli and stresses were determined, whereas the tested vertebra geometrical model used in these algorithms was considered as having a compound structure, consisting of the cancellous bone surrounded by the cortical shell. Moreover nanoindentations were conducted and an appropriate evaluation method of the obtained results was applied to extract stress-strain curves of individual lumbar spine vertebra trabecular bone struts. These data were used in the mathematical description of the vertebrae compression test. The vertebral cancellous bone structure was simulated by a beam elements network, possessing an equivalent porosity and different stiffnesses in vertical and horizontal direction. Thus, the measured course of the compression load versus the occurring specimen deformation was verified.

Structural Anomaly Detection Using Fiber Optic Sensors and Inverse Finite Element Method

NASA Technical Reports Server (NTRS)

Quach, Cuong C.; Vazquez, Sixto L.; Tessler, Alex; Moore, Jason P.; Cooper, Eric G.; Spangler, Jan. L.

2005-01-01

NASA Langley Research Center is investigating a variety of techniques for mitigating aircraft accidents due to structural component failure. One technique under consideration combines distributed fiber optic strain sensing with an inverse finite element method for detecting and characterizing structural anomalies anomalies that may provide early indication of airframe structure degradation. The technique identifies structural anomalies that result in observable changes in localized strain but do not impact the overall surface shape. Surface shape information is provided by an Inverse Finite Element Method that computes full-field displacements and internal loads using strain data from in-situ fiberoptic sensors. This paper describes a prototype of such a system and reports results from a series of laboratory tests conducted on a test coupon subjected to increasing levels of damage.

High performance computation of radiative transfer equation using the finite element method

NASA Astrophysics Data System (ADS)

Badri, M. A.; Jolivet, P.; Rousseau, B.; Favennec, Y.

2018-05-01

This article deals with an efficient strategy for numerically simulating radiative transfer phenomena using distributed computing. The finite element method alongside the discrete ordinate method is used for spatio-angular discretization of the monochromatic steady-state radiative transfer equation in an anisotropically scattering media. Two very different methods of parallelization, angular and spatial decomposition methods, are presented. To do so, the finite element method is used in a vectorial way. A detailed comparison of scalability, performance, and efficiency on thousands of processors is established for two- and three-dimensional heterogeneous test cases. Timings show that both algorithms scale well when using proper preconditioners. It is also observed that our angular decomposition scheme outperforms our domain decomposition method. Overall, we perform numerical simulations at scales that were previously unattainable by standard radiative transfer equation solvers.

Elements configuration of the open lead test circuit

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fukuzaki, Yumi, E-mail: 14514@sr.kagawa-nct.ac.jp; Ono, Akira

In the field of electronics, small electronic devices are widely utilized because they are easy to carry. The devices have various functions by user’s request. Therefore, the lead’s pitch or the ball’s pitch have been narrowed and high-density printed circuit board has been used in the devices. Use of the ICs which have narrow lead pitch makes normal connection difficult. When logic circuits in the devices are fabricated with the state-of-the-art technology, some faults have occurred more frequently. It can be divided into types of open faults and short faults. We have proposed a new test method using a testmore » circuit in the past. This paper propose elements configuration of the test circuit.« less

INTER LABORATORY COMBAT HELMET BLUNT IMPACT TEST METHOD COMPARISON

DTIC Science & Technology

2018-03-26

HELMET BLUNT IMPACT TEST METHOD COMPARISON by Tony J. Kayhart Charles A. Hewitt and Jonathan Cyganik March 2018 Final...Report March 2016 – August 2017 Approved for public release; distribution is unlimited U.S. Army Natick Soldier Research ...INTER-LABORATORY COMBAT HELMET BLUNT IMPACT TEST METHOD COMPARISON 5a. CONTRACT NUMBER 5b. GRANT NUMBER 5c. PROGRAM ELEMENT NUMBER 6. AUTHOR

Active earth pressure model tests versus finite element analysis

NASA Astrophysics Data System (ADS)

Pietrzak, Magdalena

2017-06-01

The purpose of the paper is to compare failure mechanisms observed in small scale model tests on granular sample in active state, and simulated by finite element method (FEM) using Plaxis 2D software. Small scale model tests were performed on rectangular granular sample retained by a rigid wall. Deformation of the sample resulted from simple wall translation in the direction `from the soil" (active earth pressure state. Simple Coulomb-Mohr model for soil can be helpful in interpreting experimental findings in case of granular materials. It was found that the general alignment of strain localization pattern (failure mechanism) may belong to macro scale features and be dominated by a test boundary conditions rather than the nature of the granular sample.

PERFORMANCE TESTS OF SNAP 10A THERMOELECTRIC ELEMENTS

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bergdorf, C.G.

1961-08-30

Apparatus for the performanee testing of SNAP 10A thermoelectric elements was designed, constructed, and is now in operation. Elements may be tested for any desired length of tfme up to 1400 deg F and in a vacuum of 1 x 10/ sup -5/ of Hg. The equipment used for these tcsts may also be utilized for measuring Seebeck coefficient and resistance as a function of temperature. Element performance is derived from the data on voltages and temperatures. The performance variables which are reported in graphic form are as follows: loaded output voltage at any desired DELTA T; open circuit outputmore » voltage at any desired DELTA T; power output under optimum load conditions; current produced under matched load conditions; and internal resistance of the element. (auth)« less

Fabrication and Testing of a Tubular Thermoelectric Module Based on Oxide Elements

NASA Astrophysics Data System (ADS)

Merkulov, O. V.; Politov, B. V.; Chesnokov, K. Yu.; Markov, A. A.; Leonidov, I. A.; Patrakeev, M. V.

2018-02-01

A thermoelectric module with a tubular design, considered as particularly promising for energy recovery from hot exhaust gases of vehicles was fabricated and successfully tested under laboratory condition. The module included 24 functional oxide elements synthesized by conventional methods. Elements of p-type were made of cobalt-based oxide Ca2.7Bi0.3Co4O9, while n-type elements were produced from manganese-based oxide CaMn0.95Ta0.05O3. The module was assembled using silver paste and clamping mechanisms, with the thermoelectric elements electrically connected in series by metal plates (copper, silver) and placed around the hot tube. The output power of 138 mW was achieved under a temperature gradient of 430°C with a hot side temperature of 514°C. The tested module is characterized by an advantageous manufacturing factor of 0.87 at the contact resistance of 0.44 Ω, but a relatively low output power density per volume and per mass. Possible means by which the module's operational characteristics could be improved are considered.

A Summary of the Space-Time Conservation Element and Solution Element (CESE) Method

NASA Technical Reports Server (NTRS)

Wang, Xiao-Yen J.

2015-01-01

The space-time Conservation Element and Solution Element (CESE) method for solving conservation laws is examined for its development motivation and design requirements. The characteristics of the resulting scheme are discussed. The discretization of the Euler equations is presented to show readers how to construct a scheme based on the CESE method. The differences and similarities between the CESE method and other traditional methods are discussed. The strengths and weaknesses of the method are also addressed.

A finite element method for solving the shallow water equations on the sphere

NASA Astrophysics Data System (ADS)

Comblen, Richard; Legrand, Sébastien; Deleersnijder, Eric; Legat, Vincent

Within the framework of ocean general circulation modeling, the present paper describes an efficient way to discretize partial differential equations on curved surfaces by means of the finite element method on triangular meshes. Our approach benefits from the inherent flexibility of the finite element method. The key idea consists in a dialog between a local coordinate system defined for each element in which integration takes place, and a nodal coordinate system in which all local contributions related to a vectorial degree of freedom are assembled. Since each element of the mesh and each degree of freedom are treated in the same way, the so-called pole singularity issue is fully circumvented. Applied to the shallow water equations expressed in primitive variables, this new approach has been validated against the standard test set defined by [Williamson, D.L., Drake, J.B., Hack, J.J., Jakob, R., Swarztrauber, P.N., 1992. A standard test set for numerical approximations to the shallow water equations in spherical geometry. Journal of Computational Physics 102, 211-224]. Optimal rates of convergence for the P1NC-P1 finite element pair are obtained, for both global and local quantities of interest. Finally, the approach can be extended to three-dimensional thin-layer flows in a straightforward manner.

Vibration band gaps for elastic metamaterial rods using wave finite element method

NASA Astrophysics Data System (ADS)

Nobrega, E. D.; Gautier, F.; Pelat, A.; Dos Santos, J. M. C.

2016-10-01

Band gaps in elastic metamaterial rods with spatial periodic distribution and periodically attached local resonators are investigated. New techniques to analyze metamaterial systems are using a combination of analytical or numerical method with wave propagation. One of them, called here wave spectral element method (WSEM), consists of combining the spectral element method (SEM) with Floquet-Bloch's theorem. A modern methodology called wave finite element method (WFEM), developed to calculate dynamic behavior in periodic acoustic and structural systems, utilizes a similar approach where SEM is substituted by the conventional finite element method (FEM). In this paper, it is proposed to use WFEM to calculate band gaps in elastic metamaterial rods with spatial periodic distribution and periodically attached local resonators of multi-degree-of-freedom (M-DOF). Simulated examples with band gaps generated by Bragg scattering and local resonators are calculated by WFEM and verified with WSEM, which is used as a reference method. Results are presented in the form of attenuation constant, vibration transmittance and frequency response function (FRF). For all cases, WFEM and WSEM results are in agreement, provided that the number of elements used in WFEM is sufficient to convergence. An experimental test was conducted with a real elastic metamaterial rod, manufactured with plastic in a 3D printer, without local resonance-type effect. The experimental results for the metamaterial rod with band gaps generated by Bragg scattering are compared with the simulated ones. Both numerical methods (WSEM and WFEM) can localize the band gap position and width very close to the experimental results. A hybrid approach combining WFEM with the commercial finite element software ANSYS is proposed to model complex metamaterial systems. Two examples illustrating its efficiency and accuracy to model an elastic metamaterial rod unit-cell using 1D simple rod element and 3D solid element are

78 FR 33132 - Quality Verification for Plate-Type Uranium-Aluminum Fuel Elements for Use in Research and Test...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-06-03

... Fuel Elements for Use in Research and Test Reactors AGENCY: Nuclear Regulatory Commission. ACTION... Research and Test Reactors.'' This guide describes a method that the staff of the NRC considers acceptable... assurance program for verifying the quality of plate-type uranium-aluminum fuel elements used in research...

Single element injector cold flow testing for STME swirl coaxial injector element design

NASA Technical Reports Server (NTRS)

Hulka, J.; Schneider, J. A.

1993-01-01

An oxidizer-swirled coaxial element injector is being investigated for application in the Space Transportation Main Engine (STME). Single element cold flow experiments were conducted to provide characterization of the STME injector element for future analysis, design, and optimization. All tests were conducted to quiescent, ambient backpressure conditions. Spray angle, circumferential spray uniformity, dropsize, and dropsize distribution were measured in water-only and water/nitrogen flows. Rupe mixing efficiency was measured using water/sucrose solution flows with a large grid patternator for simple comparative evaluation of mixing. Factorial designs of experiment were used for statistical evaluation of injector geometrical design features and propellant flow conditions on mixing and atomization. Increasing the free swirl angle of the liquid oxidizer had the greatest influence on increasing the mixing efficiency. The addition of gas assistance had the most significant effect on reducing oxidizer droplet size parameters and increasing droplet size distribution. Increasing the oxidizer injection velocity had the greatest influence for reducing oxidizer droplet size parameters and increasing size distribution for non-gas assisted flows. Single element and multi-element subscale hot fire testing are recommended to verify optimized designs before committing to the STME design.

Topology optimization for three-dimensional electromagnetic waves using an edge element-based finite-element method.

PubMed

Deng, Yongbo; Korvink, Jan G

2016-05-01

This paper develops a topology optimization procedure for three-dimensional electromagnetic waves with an edge element-based finite-element method. In contrast to the two-dimensional case, three-dimensional electromagnetic waves must include an additional divergence-free condition for the field variables. The edge element-based finite-element method is used to both discretize the wave equations and enforce the divergence-free condition. For wave propagation described in terms of the magnetic field in the widely used class of non-magnetic materials, the divergence-free condition is imposed on the magnetic field. This naturally leads to a nodal topology optimization method. When wave propagation is described using the electric field, the divergence-free condition must be imposed on the electric displacement. In this case, the material in the design domain is assumed to be piecewise homogeneous to impose the divergence-free condition on the electric field. This results in an element-wise topology optimization algorithm. The topology optimization problems are regularized using a Helmholtz filter and a threshold projection method and are analysed using a continuous adjoint method. In order to ensure the applicability of the filter in the element-wise topology optimization version, a regularization method is presented to project the nodal into an element-wise physical density variable.

Topology optimization for three-dimensional electromagnetic waves using an edge element-based finite-element method

PubMed Central

Korvink, Jan G.

2016-01-01

This paper develops a topology optimization procedure for three-dimensional electromagnetic waves with an edge element-based finite-element method. In contrast to the two-dimensional case, three-dimensional electromagnetic waves must include an additional divergence-free condition for the field variables. The edge element-based finite-element method is used to both discretize the wave equations and enforce the divergence-free condition. For wave propagation described in terms of the magnetic field in the widely used class of non-magnetic materials, the divergence-free condition is imposed on the magnetic field. This naturally leads to a nodal topology optimization method. When wave propagation is described using the electric field, the divergence-free condition must be imposed on the electric displacement. In this case, the material in the design domain is assumed to be piecewise homogeneous to impose the divergence-free condition on the electric field. This results in an element-wise topology optimization algorithm. The topology optimization problems are regularized using a Helmholtz filter and a threshold projection method and are analysed using a continuous adjoint method. In order to ensure the applicability of the filter in the element-wise topology optimization version, a regularization method is presented to project the nodal into an element-wise physical density variable. PMID:27279766

Control system health test system and method

DOEpatents

Hoff, Brian D.; Johnson, Kris W.; Akasam, Sivaprasad; Baker, Thomas M.

2006-08-15

A method is provided for testing multiple elements of a work machine, including a control system, a component, a sub-component that is influenced by operations of the component, and a sensor that monitors a characteristic of the sub-component. In one embodiment, the method is performed by the control system and includes sending a command to the component to adjust a first parameter associated with an operation of the component. Also, the method includes detecting a sensor signal from the sensor reflecting a second parameter associated with a characteristic of the sub-component and determining whether the second parameter is acceptable based on the command. The control system may diagnose at least one of the elements of the work machine when the second parameter of the sub-component is not acceptable.

Nuclear Cryogenic Propulsion Stage (NCPS) Fuel Element Testing in the Nuclear Thermal Rocket Element Environmental Simulator (NTREES)

NASA Technical Reports Server (NTRS)

Emrich, William J., Jr.

2017-01-01

To satisfy the Nuclear Cryogenic Propulsion Stage (NCPS) testing milestone, a graphite composite fuel element using a uranium simulant was received from the Oakridge National Lab and tested in the Nuclear Thermal Rocket Element Environmental Simulator (NTREES) at various operating conditions. The nominal operating conditions required to satisfy the milestone consisted of running the fuel element for a few minutes at a temperature of at least 2000 K with flowing hydrogen. This milestone test was successfully accomplished without incident.

On Hybrid and mixed finite element methods

NASA Technical Reports Server (NTRS)

Pian, T. H. H.

1981-01-01

Three versions of the assumed stress hybrid model in finite element methods and the corresponding variational principles for the formulation are presented. Examples of rank deficiency for stiffness matrices by the hybrid stress model are given and their corresponding kinematic deformation modes are identified. A discussion of the derivation of general semi-Loof elements for plates and shells by the hybrid stress method is given. It is shown that the equilibrium model by Fraeijs de Veubeke can be derived by the approach of the hybrid stress model as a special case of semi-Loof elements.

[Progression on finite element modeling method in scoliosis].

PubMed

Fan, Ning; Zang, Lei; Hai, Yong; Du, Peng; Yuan, Shuo

2018-04-25

Scoliosis is a complex spinal three-dimensional malformation with complicated pathogenesis, often associated with complications as thoracic deformity and shoulder imbalance. Because the acquisition of specimen or animal models are difficult, the biomechanical study of scoliosis is limited. In recent years, along with the development of the computer technology, software and image, the technology of establishing a finite element model of human spine is maturing and it has been providing strong support for the research of pathogenesis of scoliosis, the design and application of brace, and the selection of surgical methods. The finite element model method is gradually becoming an important tool in the biomechanical study of scoliosis. Establishing a high quality finite element model is the basis of analysis and future study. However, the finite element modeling process can be complex and modeling methods are greatly varied. Choosing the appropriate modeling method according to research objectives has become researchers' primary task. In this paper, the author reviews the national and international literature in recent years and concludes the finite element modeling methods in scoliosis, including data acquisition, establishment of the geometric model, the material properties, parameters setting, the validity of the finite element model validation and so on. Copyright© 2018 by the China Journal of Orthopaedics and Traumatology Press.

Comparing the NIOSH Method 5040 to a Diesel Particulate Matter Meter for Elemental Carbon

NASA Astrophysics Data System (ADS)

Ayers, David Matthew

Introduction: The sampling of elemental carbon has been associated with monitoring exposures in the trucking and mining industries. Recently, in the field of engineered nanomaterials, single wall and muti-wall carbon nanotubes (MWCNTs) are being produced in ever increasing quantities. The only approved atmospheric sampling for multi-wall carbon nanotubes in NIOSH Method 5040. These results are accurate but can take up to 30 days for sample results to be received. Objectives: Compare the results of elemental carbon sampling from the NIOSH Method 5040 to a Diesel Particulate Matter (DPM) Meter. Methods: MWCNTs were transferred and weighed between several trays placed on a scale. The NIOSH Method 5040 and DPM sampling train was hung 6 inches above the receiving tray. The transferring and weighing of the MWCNTs created an aerosol containing elemental carbon. Twenty-one total samples using both meters type were collected. Results: The assumptions for a Two-Way ANOVA were violated therefore, Mann-Whitney U Tests and a Kruskal-Wallis Test were performed. The hypotheses for both research questions were rejected. There was a significant difference in the EC concentrations obtained by the NIOSH Method 5040 and the DPM meter. There were also significant differences in elemental carbon level concentrations when sampled using a DPM meter versus a sampling pump based upon the three concentration levels (low, medium and high). Conclusions: The differences in the EC concentrations were statistically significant therefore, the two methods (NIOSH Method 5040 and DPM) are not the same. The NIOSH Method 5040 should continue to be the only authorized method of establishing an EC concentration for MWCNTs until a MWCNT specific method or an instantaneous meter is invented.

Distinct Element Modeling of the Large Block Test

NASA Astrophysics Data System (ADS)

Carlson, S. R.; Blair, S. C.; Wagoner, J. L.

2001-12-01

The Yucca Mountain Site Characterization Project is investigating Yucca Mountain, Nevada as a potential nuclear waste repository site. As part of this effort, the Large Block, a 3m x 3m x 4.5m rectangular prism of Topopah Spring tuff, was excavated at Fran Ridge near Yucca Mountain. The Large Block was heated to a peak temperature of 145\\deg C along a horizontal plane 2.75m below the top of the block over a period of about one-year. Displacements were measured in three orthogonal directions with an array of six Multiple Point Borehole Extensometers (MPBX) and were numerically simulated in three dimensions with 3DEC, a distinct element code. The distinct element method was chosen to incorporate discrete fractures in the simulations. The model domain was extended 23m below the ground surface and, in the subsurface, 23m outward from each vertical face so that fixed displacement boundary conditions could be applied well away from the heated portion of the block. A single continuum model and three distinct element models, incorporating six to twenty eight mapped fractures, were tested. Two thermal expansion coefficients were tested for the six-fracture model: a higher value taken from laboratory measurements and a lower value from an earlier field test. The MPBX data show that the largest displacements occurred in the upper portion of the block despite the higher temperatures near the center. The continuum model was found to under-predict the MPBX displacements except in the east west direction near the base of the block. The high thermal expansion model over-predicted the MPBX displacements except in the north south direction near the top of the block. The highly fractured model under-predicted most of the MPBX displacements and poorly simulated the cool-down portion of the test. Although no model provided the single best fit to all of the MPBX data, the six and seven fracture models consistently provided good fits and in most cases showed much improvement over the

Compatible-strain mixed finite element methods for incompressible nonlinear elasticity

NASA Astrophysics Data System (ADS)

Faghih Shojaei, Mostafa; Yavari, Arash

2018-05-01

We introduce a new family of mixed finite elements for incompressible nonlinear elasticity - compatible-strain mixed finite element methods (CSFEMs). Based on a Hu-Washizu-type functional, we write a four-field mixed formulation with the displacement, the displacement gradient, the first Piola-Kirchhoff stress, and a pressure-like field as the four independent unknowns. Using the Hilbert complexes of nonlinear elasticity, which describe the kinematics and the kinetics of motion, we identify the solution spaces of the independent unknown fields. In particular, we define the displacement in H1, the displacement gradient in H (curl), the stress in H (div), and the pressure field in L2. The test spaces of the mixed formulations are chosen to be the same as the corresponding solution spaces. Next, in a conforming setting, we approximate the solution and the test spaces with some piecewise polynomial subspaces of them. Among these approximation spaces are the tensorial analogues of the Nédélec and Raviart-Thomas finite element spaces of vector fields. This approach results in compatible-strain mixed finite element methods that satisfy both the Hadamard compatibility condition and the continuity of traction at the discrete level independently of the refinement level of the mesh. By considering several numerical examples, we demonstrate that CSFEMs have a good performance for bending problems and for bodies with complex geometries. CSFEMs are capable of capturing very large strains and accurately approximating stress and pressure fields. Using CSFEMs, we do not observe any numerical artifacts, e.g., checkerboarding of pressure, hourglass instability, or locking in our numerical examples. Moreover, CSFEMs provide an efficient framework for modeling heterogeneous solids.

The aggregated unfitted finite element method for elliptic problems

NASA Astrophysics Data System (ADS)

Badia, Santiago; Verdugo, Francesc; Martín, Alberto F.

2018-07-01

Unfitted finite element techniques are valuable tools in different applications where the generation of body-fitted meshes is difficult. However, these techniques are prone to severe ill conditioning problems that obstruct the efficient use of iterative Krylov methods and, in consequence, hinders the practical usage of unfitted methods for realistic large scale applications. In this work, we present a technique that addresses such conditioning problems by constructing enhanced finite element spaces based on a cell aggregation technique. The presented method, called aggregated unfitted finite element method, is easy to implement, and can be used, in contrast to previous works, in Galerkin approximations of coercive problems with conforming Lagrangian finite element spaces. The mathematical analysis of the new method states that the condition number of the resulting linear system matrix scales as in standard finite elements for body-fitted meshes, without being affected by small cut cells, and that the method leads to the optimal finite element convergence order. These theoretical results are confirmed with 2D and 3D numerical experiments.

Applications of Taylor-Galerkin finite element method to compressible internal flow problems

NASA Technical Reports Server (NTRS)

Sohn, Jeong L.; Kim, Yongmo; Chung, T. J.

1989-01-01

A two-step Taylor-Galerkin finite element method with Lapidus' artificial viscosity scheme is applied to several test cases for internal compressible inviscid flow problems. Investigations for the effect of supersonic/subsonic inlet and outlet boundary conditions on computational results are particularly emphasized.

Application of the finite element method in orthopedic implant design.

PubMed

Saha, Subrata; Roychowdhury, Amit

2009-01-01

The finite element method (FEM) was first introduced to the field of orthopedic biomechanics in the early 1970s to evaluate stresses in human bones. By the early 1980s, the method had become well established as a tool for basic research and design analysis. Since the late 1980s and early 1990s, FEM has also been used to study bone remodeling. Today, it is one of the most reliable simulation tools for evaluating wear, fatigue, crack propagation, and so forth, and is used in many types of preoperative testing. Since the introduction of FEM to orthopedic biomechanics, there have been rapid advances in computer processing speeds, the finite element and other numerical methods, understanding of mechanical properties of soft and hard tissues and their modeling, and image-processing techniques. In light of these advances, it is accepted today that FEM will continue to contribute significantly to further progress in the design and development of orthopedic implants, as well as in the understanding of other complex systems of the human body. In the following article, different main application areas of finite element simulation will be reviewed including total hip joint arthroplasty, followed by the knee, spine, shoulder, and elbow, respectively.

40 CFR 61.125 - Test methods and procedures.

Code of Federal Regulations, 2014 CFR

2014-07-01

... Section 61.125 Protection of Environment ENVIRONMENTAL PROTECTION AGENCY (CONTINUED) AIR PROGRAMS (CONTINUED) NATIONAL EMISSION STANDARDS FOR HAZARDOUS AIR POLLUTANTS National Emission Standards for Radionuclide Emissions From Elemental Phosphorus Plants § 61.125 Test methods and procedures. (a) Each owner or...

Nonconforming mortar element methods: Application to spectral discretizations

NASA Technical Reports Server (NTRS)

Maday, Yvon; Mavriplis, Cathy; Patera, Anthony

1988-01-01

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Nonlocal and Mixed-Locality Multiscale Finite Element Methods

DOE PAGES

Costa, Timothy B.; Bond, Stephen D.; Littlewood, David J.

2018-03-27

In many applications the resolution of small-scale heterogeneities remains a significant hurdle to robust and reliable predictive simulations. In particular, while material variability at the mesoscale plays a fundamental role in processes such as material failure, the resolution required to capture mechanisms at this scale is often computationally intractable. Multiscale methods aim to overcome this difficulty through judicious choice of a subscale problem and a robust manner of passing information between scales. One promising approach is the multiscale finite element method, which increases the fidelity of macroscale simulations by solving lower-scale problems that produce enriched multiscale basis functions. Here, inmore » this study, we present the first work toward application of the multiscale finite element method to the nonlocal peridynamic theory of solid mechanics. This is achieved within the context of a discontinuous Galerkin framework that facilitates the description of material discontinuities and does not assume the existence of spatial derivatives. Analysis of the resulting nonlocal multiscale finite element method is achieved using the ambulant Galerkin method, developed here with sufficient generality to allow for application to multiscale finite element methods for both local and nonlocal models that satisfy minimal assumptions. Finally, we conclude with preliminary results on a mixed-locality multiscale finite element method in which a nonlocal model is applied at the fine scale and a local model at the coarse scale.« less

Nonlocal and Mixed-Locality Multiscale Finite Element Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Costa, Timothy B.; Bond, Stephen D.; Littlewood, David J.

In many applications the resolution of small-scale heterogeneities remains a significant hurdle to robust and reliable predictive simulations. In particular, while material variability at the mesoscale plays a fundamental role in processes such as material failure, the resolution required to capture mechanisms at this scale is often computationally intractable. Multiscale methods aim to overcome this difficulty through judicious choice of a subscale problem and a robust manner of passing information between scales. One promising approach is the multiscale finite element method, which increases the fidelity of macroscale simulations by solving lower-scale problems that produce enriched multiscale basis functions. Here, inmore » this study, we present the first work toward application of the multiscale finite element method to the nonlocal peridynamic theory of solid mechanics. This is achieved within the context of a discontinuous Galerkin framework that facilitates the description of material discontinuities and does not assume the existence of spatial derivatives. Analysis of the resulting nonlocal multiscale finite element method is achieved using the ambulant Galerkin method, developed here with sufficient generality to allow for application to multiscale finite element methods for both local and nonlocal models that satisfy minimal assumptions. Finally, we conclude with preliminary results on a mixed-locality multiscale finite element method in which a nonlocal model is applied at the fine scale and a local model at the coarse scale.« less

Optimal least-squares finite element method for elliptic problems

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Povinelli, Louis A.

1991-01-01

An optimal least squares finite element method is proposed for two dimensional and three dimensional elliptic problems and its advantages are discussed over the mixed Galerkin method and the usual least squares finite element method. In the usual least squares finite element method, the second order equation (-Delta x (Delta u) + u = f) is recast as a first order system (-Delta x p + u = f, Delta u - p = 0). The error analysis and numerical experiment show that, in this usual least squares finite element method, the rate of convergence for flux p is one order lower than optimal. In order to get an optimal least squares method, the irrotationality Delta x p = 0 should be included in the first order system.

An evaluation of objective rating methods for full-body finite element model comparison to PMHS tests.

PubMed

Vavalle, Nicholas A; Jelen, Benjamin C; Moreno, Daniel P; Stitzel, Joel D; Gayzik, F Scott

2013-01-01

Objective evaluation methods of time history signals are used to quantify how well simulated human body responses match experimental data. As the use of simulations grows in the field of biomechanics, there is a need to establish standard approaches for comparisons. There are 2 aims of this study. The first is to apply 3 objective evaluation methods found in the literature to a set of data from a human body finite element model. The second is to compare the results of each method, examining how they are correlated to each other and the relative strengths and weaknesses of the algorithms. In this study, the methods proposed by Sprague and Geers (magnitude and phase error, SGM and SGP), Rhule et al. (cumulative standard deviation, CSD), and Gehre et al. (CORrelation and Analysis, or CORA, size, phase, shape, corridor) were compared. A 40 kph frontal sled test presented by Shaw et al. was simulated using the Global Human Body Models Consortium midsized male full-body finite element model (v. 3.5). Mean and standard deviation experimental data (n = 5) from Shaw et al. were used as the benchmark. Simulated data were output from the model at the appropriate anatomical locations for kinematic comparison. Force data were output at the seat belts, seat pan, knee, and foot restraints. Objective comparisons from 53 time history data channels were compared to the experimental results. To compare the different methods, all objective comparison metrics were cross-plotted and linear regressions were calculated. The following ratings were found to be statistically significantly correlated (P < .01): SGM and CORrelation and Analysis (CORA) size, R (2) = 0.73; SGP and CORA shape, R (2) = 0.82; and CSD and CORA's corridor factor, R (2) = 0.59. Relative strengths of the correlated ratings were then investigated. For example, though correlated to CORA size, SGM carries a sign to indicate whether the simulated response is greater than or less than the benchmark signal. A further

Least-squares finite element methods for compressible Euler equations

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Carey, G. F.

1990-01-01

A method based on backward finite differencing in time and a least-squares finite element scheme for first-order systems of partial differential equations in space is applied to the Euler equations for gas dynamics. The scheme minimizes the L-sq-norm of the residual within each time step. The method naturally generates numerical dissipation proportional to the time step size. An implicit method employing linear elements has been implemented and proves robust. For high-order elements, computed solutions based on the L-sq method may have oscillations for calculations at similar time step sizes. To overcome this difficulty, a scheme which minimizes the weighted H1-norm of the residual is proposed and leads to a successful scheme with high-degree elements. Finally, a conservative least-squares finite element method is also developed. Numerical results for two-dimensional problems are given to demonstrate the shock resolution of the methods and compare different approaches.

Analysis of concrete beams using applied element method

NASA Astrophysics Data System (ADS)

Lincy Christy, D.; Madhavan Pillai, T. M.; Nagarajan, Praveen

2018-03-01

The Applied Element Method (AEM) is a displacement based method of structural analysis. Some of its features are similar to that of Finite Element Method (FEM). In AEM, the structure is analysed by dividing it into several elements similar to FEM. But, in AEM, elements are connected by springs instead of nodes as in the case of FEM. In this paper, background to AEM is discussed and necessary equations are derived. For illustrating the application of AEM, it has been used to analyse plain concrete beam of fixed support condition. The analysis is limited to the analysis of 2-dimensional structures. It was found that the number of springs has no much influence on the results. AEM could predict deflection and reactions with reasonable degree of accuracy.

Finite Elements Analysis of a Composite Semi-Span Test Article With and Without Discrete Damage

NASA Technical Reports Server (NTRS)

Lovejoy, Andrew E.; Jegley, Dawn C. (Technical Monitor)

2000-01-01

AS&M Inc. performed finite element analysis, with and without discrete damage, of a composite semi-span test article that represents the Boeing 220-passenger transport aircraft composite semi-span test article. A NASTRAN bulk data file and drawings of the test mount fixtures and semi-span components were utilized to generate the baseline finite element model. In this model, the stringer blades are represented by shell elements, and the stringer flanges are combined with the skin. Numerous modeling modifications and discrete source damage scenarios were applied to the test article model throughout the course of the study. This report details the analysis method and results obtained from the composite semi-span study. Analyses were carried out for three load cases: Braked Roll, LOG Down-Bending and 2.5G Up-Bending. These analyses included linear and nonlinear static response, as well as linear and nonlinear buckling response. Results are presented in the form of stress and strain plots. factors of safety for failed elements, buckling loads and modes, deflection prediction tables and plots, and strainage prediction tables and plots. The collected results are presented within this report for comparison to test results.

High speed inviscid compressible flow by the finite element method

NASA Technical Reports Server (NTRS)

Zienkiewicz, O. C.; Loehner, R.; Morgan, K.

1984-01-01

The finite element method and an explicit time stepping algorithm which is based on Taylor-Galerkin schemes with an appropriate artificial viscosity is combined with an automatic mesh refinement process which is designed to produce accurate steady state solutions to problems of inviscid compressible flow in two dimensions. The results of two test problems are included which demonstrate the excellent performance characteristics of the proposed procedures.

A New High-Order Spectral Difference Method for Simulating Viscous Flows on Unstructured Grids with Mixed Elements

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Mao; Qiu, Zihua; Liang, Chunlei

In the present study, a new spectral difference (SD) method is developed for viscous flows on meshes with a mixture of triangular and quadrilateral elements. The standard SD method for triangular elements, which employs Lagrangian interpolating functions for fluxes, is not stable when the designed accuracy of spatial discretization is third-order or higher. Unlike the standard SD method, the method examined here uses vector interpolating functions in the Raviart-Thomas (RT) spaces to construct continuous flux functions on reference elements. Studies have been performed for 2D wave equation and Euler equa- tions. Our present results demonstrated that the SDRT method ismore » stable and high-order accurate for a number of test problems by using triangular-, quadrilateral-, and mixed- element meshes.« less

Spectral/ hp element methods: Recent developments, applications, and perspectives

NASA Astrophysics Data System (ADS)

Xu, Hui; Cantwell, Chris D.; Monteserin, Carlos; Eskilsson, Claes; Engsig-Karup, Allan P.; Sherwin, Spencer J.

2018-02-01

The spectral/ hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate a C 0 - continuous expansion. Computationally and theoretically, by increasing the polynomial order p, high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/ hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use of the spectral/ hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order to use the spectral/ hp element method in more complex science and engineering applications are discussed.

Method and system for processing optical elements using magnetorheological finishing

DOEpatents

Menapace, Joseph Arthur; Schaffers, Kathleen Irene; Bayramian, Andrew James; Molander, William A

2012-09-18

A method of finishing an optical element includes mounting the optical element in an optical mount having a plurality of fiducials overlapping with the optical element and obtaining a first metrology map for the optical element and the plurality of fiducials. The method also includes obtaining a second metrology map for the optical element without the plurality of fiducials, forming a difference map between the first metrology map and the second metrology map, and aligning the first metrology map and the second metrology map. The method further includes placing mathematical fiducials onto the second metrology map using the difference map to form a third metrology map and associating the third metrology map to the optical element. Moreover, the method includes mounting the optical element in the fixture in an MRF tool, positioning the optical element in the fixture; removing the plurality of fiducials, and finishing the optical element.

Dynamic analysis of suspension cable based on vector form intrinsic finite element method

NASA Astrophysics Data System (ADS)

Qin, Jian; Qiao, Liang; Wan, Jiancheng; Jiang, Ming; Xia, Yongjun

2017-10-01

A vector finite element method is presented for the dynamic analysis of cable structures based on the vector form intrinsic finite element (VFIFE) and mechanical properties of suspension cable. Firstly, the suspension cable is discretized into different elements by space points, the mass and external forces of suspension cable are transformed into space points. The structural form of cable is described by the space points at different time. The equations of motion for the space points are established according to the Newton’s second law. Then, the element internal forces between the space points are derived from the flexible truss structure. Finally, the motion equations of space points are solved by the central difference method with reasonable time integration step. The tangential tension of the bearing rope in a test ropeway with the moving concentrated loads is calculated and compared with the experimental data. The results show that the tangential tension of suspension cable with moving loads is consistent with the experimental data. This method has high calculated precision and meets the requirements of engineering application.

The Blended Finite Element Method for Multi-fluid Plasma Modeling

DTIC Science & Technology

2016-07-01

Briefing Charts 3. DATES COVERED (From - To) 07 June 2016 - 01 July 2016 4. TITLE AND SUBTITLE The Blended Finite Element Method for Multi-fluid Plasma...BLENDED FINITE ELEMENT METHOD FOR MULTI-FLUID PLASMA MODELING Éder M. Sousa1, Uri Shumlak2 1ERC INC., IN-SPACE PROPULSION BRANCH (RQRS) AIR FORCE RESEARCH...MULTI-FLUID PLASMA MODEL 2 BLENDED FINITE ELEMENT METHOD Blended Finite Element Method Nodal Continuous Galerkin Modal Discontinuous Galerkin Model

Large-scale computation of incompressible viscous flow by least-squares finite element method

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Lin, T. L.; Povinelli, Louis A.

1993-01-01

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

A Multifunctional Interface Method for Coupling Finite Element and Finite Difference Methods: Two-Dimensional Scalar-Field Problems

NASA Technical Reports Server (NTRS)

Ransom, Jonathan B.

2002-01-01

A multifunctional interface method with capabilities for variable-fidelity modeling and multiple method analysis is presented. The methodology provides an effective capability by which domains with diverse idealizations can be modeled independently to exploit the advantages of one approach over another. The multifunctional method is used to couple independently discretized subdomains, and it is used to couple the finite element and the finite difference methods. The method is based on a weighted residual variational method and is presented for two-dimensional scalar-field problems. A verification test problem and a benchmark application are presented, and the computational implications are discussed.

Finite elements and the method of conjugate gradients on a concurrent processor

NASA Technical Reports Server (NTRS)

Lyzenga, G. A.; Raefsky, A.; Hager, G. H.

1985-01-01

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90 percent for sufficiently large problems.

Finite elements and the method of conjugate gradients on a concurrent processor

NASA Technical Reports Server (NTRS)

Lyzenga, G. A.; Raefsky, A.; Hager, B. H.

1984-01-01

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90% for sufficiently large problems.

Finite element method formulation in polar coordinates for transient heat conduction problems

NASA Astrophysics Data System (ADS)

Duda, Piotr

2016-04-01

The aim of this paper is the formulation of the finite element method in polar coordinates to solve transient heat conduction problems. It is hard to find in the literature a formulation of the finite element method (FEM) in polar or cylindrical coordinates for the solution of heat transfer problems. This document shows how to apply the most often used boundary conditions. The global equation system is solved by the Crank-Nicolson method. The proposed algorithm is verified in three numerical tests. In the first example, the obtained transient temperature distribution is compared with the temperature obtained from the presented analytical solution. In the second numerical example, the variable boundary condition is assumed. In the last numerical example the component with the shape different than cylindrical is used. All examples show that the introduction of the polar coordinate system gives better results than in the Cartesian coordinate system. The finite element method formulation in polar coordinates is valuable since it provides a higher accuracy of the calculations without compacting the mesh in cylindrical or similar to tubular components. The proposed method can be applied for circular elements such as boiler drums, outlet headers, flux tubes. This algorithm can be useful during the solution of inverse problems, which do not allow for high density grid. This method can calculate the temperature distribution in the bodies of different properties in the circumferential and the radial direction. The presented algorithm can be developed for other coordinate systems. The examples demonstrate a good accuracy and stability of the proposed method.

Hybrid finite element and Brownian dynamics method for charged particles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huber, Gary A., E-mail: ghuber@ucsd.edu; Miao, Yinglong; Zhou, Shenggao

2016-04-28

Diffusion is often the rate-determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. A previous study introduced a new hybrid diffusion method that couples the strengths of each of these two methods, but was limited by the lack of interactions among the particles; the force on each particle had to be from an external field. This study further develops the method to allow charged particles. The method is derived for a general multidimensional system and is presented usingmore » a basic test case for a one-dimensional linear system with one charged species and a radially symmetric system with three charged species.« less

Structural Element Testing in Support of the Design of the NASA Composite Crew Module

NASA Technical Reports Server (NTRS)

Kellas, Sotiris; Jackson, Wade C.; Thesken, John C.; Schleicher, Eric; Wagner, Perry; Kirsch, Michael T.

2012-01-01

In January 2007, the NASA Administrator and Associate Administrator for the Exploration Systems Mission Directorate chartered the NASA Engineering and Safety Center (NESC) to design, build, and test a full-scale Composite Crew Module (CCM). For the design and manufacturing of the CCM, the team adopted the building block approach where design and manufacturing risks were mitigated through manufacturing trials and structural testing at various levels of complexity. Following NASA's Structural Design Verification Requirements, a further objective was the verification of design analysis methods and the provision of design data for critical structural features. Test articles increasing in complexity from basic material characterization coupons through structural feature elements and large structural components, to full-scale structures were evaluated. This paper discusses only four elements tests three of which include joints and one that includes a tapering honeycomb core detail. For each test series included are specimen details, instrumentation, test results, a brief analysis description, test analysis correlation and conclusions.

A multidimensional finite element method for CFD

NASA Technical Reports Server (NTRS)

Pepper, Darrell W.; Humphrey, Joseph W.

1991-01-01

A finite element method is used to solve the equations of motion for 2- and 3-D fluid flow. The time-dependent equations are solved explicitly using quadrilateral (2-D) and hexahedral (3-D) elements, mass lumping, and reduced integration. A Petrov-Galerkin technique is applied to the advection terms. The method requires a minimum of computational storage, executes quickly, and is scalable for execution on computer systems ranging from PCs to supercomputers.

Method and apparatus for staking optical elements

DOEpatents

Woods, Robert O.

1988-01-01

A method and apparatus for staking two optical elements together in order to retain their alignment is disclosed. The apparatus includes a removable adaptor made up of first and second adaptor bodies each having a lateral slot in their front and side faces. The adaptor also includes a system for releasably attaching each adaptor body to a respective optical element such that when the two optical elements are positioned relative to one another the adaptor bodies are adjacent and the lateral slots therein are aligned to form key slots. The adaptor includes keys which are adapted to fit into the key slots. A curable filler material is employed to retain the keys in the key slots and thereby join the first and second adaptor bodies to form the adaptor. Also disclosed is a method for staking together two optical elements employing the adaptor of the present invention.

Method and apparatus for staking optical elements

DOEpatents

Woods, Robert O.

1988-10-04

A method and apparatus for staking two optical elements together in order to retain their alignment is disclosed. The apparatus includes a removable adaptor made up of first and second adaptor bodies each having a lateral slot in their front and side faces. The adaptor also includes a system for releasably attaching each adaptor body to a respective optical element such that when the two optical elements are positioned relative to one another the adaptor bodies are adjacent and the lateral slots therein are aligned to form key slots. The adaptor includes keys which are adapted to fit into the key slots. A curable filler material is employed to retain the keys in the key slots and thereby join the first and second adaptor bodies to form the adaptor. Also disclosed is a method for staking together two optical elements employing the adaptor of the present invention.

Elemental composition of edible nuts: fast optimization and validation procedure of an ICP-OES method.

PubMed

Tošić, Snežana B; Mitić, Snežana S; Velimirović, Dragan S; Stojanović, Gordana S; Pavlović, Aleksandra N; Pecev-Marinković, Emilija T

2015-08-30

An inductively coupled plasma-optical emission spectrometry method for the speedy simultaneous detection of 19 elements in edible nuts (walnuts: Juglans nigra; almonds: Prunus dulcis; hazelnuts: Corylus avellana; Brazil nuts: Bertholletia excelsa; cashews: Anacardium occidentalle; pistachios: Pistacia vera; and peanuts: Arachis hypogaea) available on the Serbian markets, was optimized and validated through the selection of instrumental parameters and analytical lines free from spectral interference and with the lowest matrix effects. The analysed macro-elements were present in the following descending order: Na > Mg > Ca > K. Of all the trace elements, the tested samples showed the highest content of Fe. The micro-element Se was detected in all the samples of nuts. The toxic elements As, Cd and Pb were either not detected or the contents were below the limit of detection. One-way analysis of variance, Student's t-test, Tukey's HSD post hoc test and hierarchical agglomerative cluster analysis were applied in the statistical analysis of the results. Based on the detected content of analysed elements it can be concluded that nuts may be a good additional source of minerals as micronutrients. © 2014 Society of Chemical Industry.

NOTE: Solving the ECG forward problem by means of a meshless finite element method

NASA Astrophysics Data System (ADS)

Li, Z. S.; Zhu, S. A.; He, Bin

2007-07-01

The conventional numerical computational techniques such as the finite element method (FEM) and the boundary element method (BEM) require laborious and time-consuming model meshing. The new meshless FEM only uses the boundary description and the node distribution and no meshing of the model is required. This paper presents the fundamentals and implementation of meshless FEM and the meshless FEM method is adapted to solve the electrocardiography (ECG) forward problem. The method is evaluated on a single-layer torso model, in which the analytical solution exists, and tested in a realistic geometry homogeneous torso model, with satisfactory results being obtained. The present results suggest that the meshless FEM may provide an alternative for ECG forward solutions.

Numerical sedimentation particle-size analysis using the Discrete Element Method

NASA Astrophysics Data System (ADS)

Bravo, R.; Pérez-Aparicio, J. L.; Gómez-Hernández, J. J.

2015-12-01

Sedimentation tests are widely used to determine the particle size distribution of a granular sample. In this work, the Discrete Element Method interacts with the simulation of flow using the well known one-way-coupling method, a computationally affordable approach for the time-consuming numerical simulation of the hydrometer, buoyancy and pipette sedimentation tests. These tests are used in the laboratory to determine the particle-size distribution of fine-grained aggregates. Five samples with different particle-size distributions are modeled by about six million rigid spheres projected on two-dimensions, with diameters ranging from 2.5 ×10-6 m to 70 ×10-6 m, forming a water suspension in a sedimentation cylinder. DEM simulates the particle's movement considering laminar flow interactions of buoyant, drag and lubrication forces. The simulation provides the temporal/spatial distributions of densities and concentrations of the suspension. The numerical simulations cannot replace the laboratory tests since they need the final granulometry as initial data, but, as the results show, these simulations can identify the strong and weak points of each method and eventually recommend useful variations and draw conclusions on their validity, aspects very difficult to achieve in the laboratory.

Structural Health Monitoring Using High-Density Fiber Optic Strain Sensor and Inverse Finite Element Methods

NASA Technical Reports Server (NTRS)

Vazquez, Sixto L.; Tessler, Alexander; Quach, Cuong C.; Cooper, Eric G.; Parks, Jeffrey; Spangler, Jan L.

2005-01-01

In an effort to mitigate accidents due to system and component failure, NASA s Aviation Safety has partnered with industry, academia, and other governmental organizations to develop real-time, on-board monitoring capabilities and system performance models for early detection of airframe structure degradation. NASA Langley is investigating a structural health monitoring capability that uses a distributed fiber optic strain system and an inverse finite element method for measuring and modeling structural deformations. This report describes the constituent systems that enable this structural monitoring function and discusses results from laboratory tests using the fiber strain sensor system and the inverse finite element method to demonstrate structural deformation estimation on an instrumented test article

Dynamic correction of the laser beam coordinate in fabrication of large-sized diffractive elements for testing aspherical mirrors

NASA Astrophysics Data System (ADS)

Shimansky, R. V.; Poleshchuk, A. G.; Korolkov, V. P.; Cherkashin, V. V.

2017-05-01

This paper presents a method of improving the accuracy of a circular laser system in fabrication of large-diameter diffractive optical elements by means of a polar coordinate system and the results of their use. An algorithm for correcting positioning errors of a circular laser writing system developed at the Institute of Automation and Electrometry, SB RAS, is proposed and tested. Highprecision synthesized holograms fabricated by this method and the results of using these elements for testing the 6.5 m diameter aspheric mirror of the James Webb space telescope (JWST) are described..

Ablative Thermal Response Analysis Using the Finite Element Method

NASA Technical Reports Server (NTRS)

Dec John A.; Braun, Robert D.

2009-01-01

A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.

Production test IP-376-D, Supplement B Irradiation of MGCR-HDR-3 Test Element

DOE Office of Scientific and Technical Information (OSTI.GOV)

Baars, R.E.

The objective of this supplement to PT-IP-376-D, Irradiation of MGCR-HDR-3 Test Element is to authorize 1000 hours of operation at a maximum test specimen surface temperature of 1700 F. The original production test authorized a test duration of four months at a maximum specimen surface temperature of 1500 F; supplement A authorized extension of the test duration to ten months. The desired increase in surface temperature is requested to demonstrate the general feasibility of operation of the fuel element at 1700 F, and to obtain specific information on the performance of Hastelloy-X cladding and fuel bodies. The increased temperature hasmore » been approved by the Atomic Energy Commission.« less

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

NASA Astrophysics Data System (ADS)

Gao, Kai; Fu, Shubin; Chung, Eric T.

2018-05-01

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss-Lobatto-Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

Application of the Finite Element Method to Rotary Wing Aeroelasticity

NASA Technical Reports Server (NTRS)

Straub, F. K.; Friedmann, P. P.

1982-01-01

A finite element method for the spatial discretization of the dynamic equations of equilibrium governing rotary-wing aeroelastic problems is presented. Formulation of the finite element equations is based on weighted Galerkin residuals. This Galerkin finite element method reduces algebraic manipulative labor significantly, when compared to the application of the global Galerkin method in similar problems. The coupled flap-lag aeroelastic stability boundaries of hingeless helicopter rotor blades in hover are calculated. The linearized dynamic equations are reduced to the standard eigenvalue problem from which the aeroelastic stability boundaries are obtained. The convergence properties of the Galerkin finite element method are studied numerically by refining the discretization process. Results indicate that four or five elements suffice to capture the dynamics of the blade with the same accuracy as the global Galerkin method.

Addressable test matrix for measuring analog transfer characteristics of test elements used for integrated process control and device evaluation

NASA Technical Reports Server (NTRS)

Buehler, Martin G. (Inventor)

1988-01-01

A set of addressable test structures, each of which uses addressing schemes to access individual elements of the structure in a matrix, is used to test the quality of a wafer before integrated circuits produced thereon are diced, packaged and subjected to final testing. The electrical characteristic of each element is checked and compared to the electrical characteristic of all other like elements in the matrix. The effectiveness of the addressable test matrix is in readily analyzing the electrical characteristics of the test elements and in providing diagnostic information.

Ground shake test of the UH-60A helicopter airframe and comparison with NASTRAN finite element model predictions

NASA Technical Reports Server (NTRS)

Howland, G. R.; Durno, J. A.; Twomey, W. J.

1990-01-01

Sikorsky Aircraft, together with the other major helicopter airframe manufacturers, is engaged in a study to improve the use of finite element analysis to predict the dynamic behavior of helicopter airframes, under a rotorcraft structural dynamics program called DAMVIBS (Design Analysis Methods for VIBrationS), sponsored by the NASA-Langley. The test plan and test results are presented for a shake test of the UH-60A BLACK HAWK helicopter. A comparison is also presented of test results with results obtained from analysis using a NASTRAN finite element model.

Fatigue Testing of Wing Beam by the Resonance Method

NASA Technical Reports Server (NTRS)

Bleakney, William M

1938-01-01

Preliminary fatigue tests on two aluminum-alloy wing-beam specimens subjected to reversed axial loading are described. The motion used consists in incorporating one or two reciprocating motors in a resonance system of which the specimen is the spring element. A description is given of the reciprocating motors, and of the method of assembling and adjusting the vibrating system. The results indicate that the method is well adapted to fatigue tests of not only uniform wing beams but also wing beams with asymmetrical local reinforcements.

Testing of the Defense Waste Processing Facility Cold Chemical Dissolution Method in Sludge Batch 9 Qualification

DOE Office of Scientific and Technical Information (OSTI.GOV)

Edwards, T.; Pareizs, J.; Coleman, C.

For each sludge batch that is processed in the Defense Waste Processing Facility (DWPF), the Savannah River National Laboratory (SRNL) tests the applicability of the digestion methods used by the DWPF Laboratory for elemental analysis of Sludge Receipt and Adjustment Tank (SRAT) Receipt samples and SRAT Product process control samples. DWPF SRAT samples are typically dissolved using a method referred to as the DWPF Cold Chemical or Cold Chem Method (CC), (see DWPF Procedure SW4- 15.201). Testing indicates that the CC method produced mixed results. The CC method did not result in complete dissolution of either the SRAT Receipt ormore » SRAT Product with some fine, dark solids remaining. However, elemental analyses did not reveal extreme biases for the major elements in the sludge when compared with analyses obtained following dissolution by hot aqua regia (AR) or sodium peroxide fusion (PF) methods. The CC elemental analyses agreed with the AR and PF methods well enough that it should be adequate for routine process control analyses in the DWPF after much more extensive side-by-side tests of the CC method and the PF method are performed on the first 10 SRAT cycles of the Sludge Batch 9 (SB9) campaign. The DWPF Laboratory should continue with their plans for further tests of the CC method during these 10 SRAT cycles.« less

Experimental validation of finite element and boundary element methods for predicting structural vibration and radiated noise

NASA Technical Reports Server (NTRS)

Seybert, A. F.; Wu, T. W.; Wu, X. F.

1994-01-01

This research report is presented in three parts. In the first part, acoustical analyses were performed on modes of vibration of the housing of a transmission of a gear test rig developed by NASA. The modes of vibration of the transmission housing were measured using experimental modal analysis. The boundary element method (BEM) was used to calculate the sound pressure and sound intensity on the surface of the housing and the radiation efficiency of each mode. The radiation efficiency of each of the transmission housing modes was then compared to theoretical results for a finite baffled plate. In the second part, analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise level radiated from the box. The FEM was used to predict the vibration, while the BEM was used to predict the sound intensity and total radiated sound power using surface vibration as the input data. Vibration predicted by the FEM model was validated by experimental modal analysis; noise predicted by the BEM was validated by measurements of sound intensity. Three types of results are presented for the total radiated sound power: sound power predicted by the BEM model using vibration data measured on the surface of the box; sound power predicted by the FEM/BEM model; and sound power measured by an acoustic intensity scan. In the third part, the structure used in part two was modified. A rib was attached to the top plate of the structure. The FEM and BEM were then used to predict structural vibration and radiated noise respectively. The predicted vibration and radiated noise were then validated through experimentation.

An enriched finite element method to fractional advection-diffusion equation

NASA Astrophysics Data System (ADS)

Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam

2017-08-01

In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.

Nuclear Cryogenic Propulsion Stage (NCPS) Fuel Element Testing in the Nuclear Thermal Rocket Element Environmental Simulator (NTREES)

NASA Technical Reports Server (NTRS)

Emrich, William J., Jr.

2017-01-01

To support the on-going nuclear thermal propulsion effort, a state-of-the-art non nuclear experimental test setup has been constructed to evaluate the performance characteristics of candidate fuel element materials and geometries in representative environments. The facility to perform this testing is referred to as the Nuclear Thermal Rocket Element Environment Simulator (NTREES). Last year NTREES was successfully used to satisfy a testing milestone for the Nuclear Cryogenic Propulsion Stage (NCPS) project and met or exceeded all required objectives.

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gao, Kai; Fu, Shubin; Chung, Eric T.

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructsmore » high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.« less

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

DOE PAGES

Gao, Kai; Fu, Shubin; Chung, Eric T.

2018-02-04

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructsmore » high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss–Lobatto–Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.« less

Stabilised finite-element methods for solving the level set equation with mass conservation

NASA Astrophysics Data System (ADS)

Kabirou Touré, Mamadou; Fahsi, Adil; Soulaïmani, Azzeddine

2016-01-01

Finite-element methods are studied for solving moving interface flow problems using the level set approach and a stabilised variational formulation proposed in Touré and Soulaïmani (2012; Touré and Soulaïmani To appear in 2016), coupled with a level set correction method. The level set correction is intended to enhance the mass conservation satisfaction property. The stabilised variational formulation (Touré and Soulaïmani 2012; Touré and Soulaïmani, To appear in 2016) constrains the level set function to remain close to the signed distance function, while the mass conservation is a correction step which enforces the mass balance. The eXtended finite-element method (XFEM) is used to take into account the discontinuities of the properties within an element. XFEM is applied to solve the Navier-Stokes equations for two-phase flows. The numerical methods are numerically evaluated on several test cases such as time-reversed vortex flow, a rigid-body rotation of Zalesak's disc, sloshing flow in a tank, a dam-break over a bed, and a rising bubble subjected to buoyancy. The numerical results show the importance of satisfying global mass conservation to accurately capture the interface position.

Transforming Mean and Osculating Elements Using Numerical Methods

NASA Technical Reports Server (NTRS)

Ely, Todd A.

2010-01-01

Mean element propagation of perturbed two body orbits has as its mathematical basis averaging theory of nonlinear dynamical systems. Averaged mean elements define the long-term evolution characteristics of an orbit. Using averaging theory, a near identity transformation can be found that transforms the mean elements back to the osculating elements that contain short period terms in addition to the secular and long period mean elements. The ability to perform the conversion is necessary so that orbit design conducted in mean elements can be converted back into osculating results. In the present work, this near identity transformation is found using the Fast Fourier Transform. An efficient method is found that is capable of recovering the osculating elements to first order

A method of selecting grid size to account for Hertz deformation in finite element analysis of spur gears

NASA Technical Reports Server (NTRS)

Coy, J. J.; Chao, C. H. C.

1981-01-01

A method of selecting grid size for the finite element analysis of gear tooth deflection is presented. The method is based on a finite element study of two cylinders in line contact, where the criterion for establishing element size was that there be agreement with the classical Hertzian solution for deflection. The results are applied to calculate deflection for the gear specimen used in the NASA spur gear test rig. Comparisons are made between the present results and the results of two other methods of calculation. The results have application in design of gear tooth profile modifications to reduce noise and dynamic loads.

[Application of finite element method in spinal biomechanics].

PubMed

Liu, Qiang; Zhang, Jun; Sun, Shu-Chun; Wang, Fei

2017-02-25

The finite element model is one of the most important methods in study of modern spinal biomechanics, according to the needs to simulate the various states of the spine, calculate the stress force and strain distribution of the different groups in the state, and explore its principle of mechanics, mechanism of injury, and treatment effectiveness. In addition, in the study of the pathological state of the spine, the finite element is mainly used in the understanding the mechanism of lesion location, evaluating the effects of different therapeutic tool, assisting and completing the selection and improvement of therapeutic tool, in order to provide a theoretical basis for the rehabilitation of spinal lesions. Finite element method can be more provide the service for the patients suffering from spinal correction, operation and individual implant design. Among the design and performance evaluation of the implant need to pay attention to the individual difference and perfect the evaluation system. At present, how to establish a model which is more close to the real situation has been the focus and difficulty of the study of human body's finite element.Although finite element method can better simulate complex working condition, it is necessary to improve the authenticity of the model and the sharing of the group by using many kinds of methods, such as image science, statistics, kinematics and so on. Copyright© 2017 by the China Journal of Orthopaedics and Traumatology Press.

Finite element methods on supercomputers - The scatter-problem

NASA Technical Reports Server (NTRS)

Loehner, R.; Morgan, K.

1985-01-01

Certain problems arise in connection with the use of supercomputers for the implementation of finite-element methods. These problems are related to the desirability of utilizing the power of the supercomputer as fully as possible for the rapid execution of the required computations, taking into account the gain in speed possible with the aid of pipelining operations. For the finite-element method, the time-consuming operations may be divided into three categories. The first two present no problems, while the third type of operation can be a reason for the inefficient performance of finite-element programs. Two possibilities for overcoming certain difficulties are proposed, giving attention to a scatter-process.

Time Triggered Ethernet System Testing Means and Method

NASA Technical Reports Server (NTRS)

Smithgall, William Todd (Inventor); Hall, Brendan (Inventor); Varadarajan, Srivatsan (Inventor)

2014-01-01

Methods and apparatus are provided for evaluating the performance of a Time Triggered Ethernet (TTE) system employing Time Triggered (TT) communication. A real TTE system under test (SUT) having real input elements communicating using TT messages with output elements via one or more first TTE switches during a first time interval schedule established for the SUT. A simulation system is also provided having input simulators that communicate using TT messages via one or more second TTE switches with the same output elements during a second time interval schedule established for the simulation system. The first and second time interval schedules are off-set slightly so that messages from the input simulators, when present, arrive at the output elements prior to messages from the analogous real inputs, thereby having priority over messages from the real inputs and causing the system to operate based on the simulated inputs when present.

Nuclear fuel elements and method of making same

DOEpatents

Schweitzer, Donald G.

1992-01-01

A nuclear fuel element for a high temperature gas nuclear reactor that has an average operating temperature in excess of 2000.degree. C., and a method of making such a fuel element. The fuel element is characterized by having fissionable fuel material localized and stabilized within pores of a carbon or graphite member by melting the fissionable material to cause it to chemically react with the carbon walls of the pores. The fissionable fuel material is further stabilized and localized within the pores of the graphite member by providing one or more coatings of pyrolytic carbon or diamond surrounding the porous graphite member so that each layer defines a successive barrier against migration of the fissionable fuel from the pores, and so that the outermost layer of pyrolytic carbon or diamond forms a barrier between the fissionable material and the moderating gases used in an associated high temperature gas reactor. The method of the invention provides for making such new elements either as generally spherically elements, or as flexible filaments, or as other relatively small-sized fuel elements that are particularly suited for use in high temperature gas reactors.

A weak Galerkin generalized multiscale finite element method

DOE PAGES

Mu, Lin; Wang, Junping; Ye, Xiu

2016-03-31

In this study, we propose a general framework for weak Galerkin generalized multiscale (WG-GMS) finite element method for the elliptic problems with rapidly oscillating or high contrast coefficients. This general WG-GMS method features in high order accuracy on general meshes and can work with multiscale basis derived by different numerical schemes. A special case is studied under this WG-GMS framework in which the multiscale basis functions are obtained by solving local problem with the weak Galerkin finite element method. Convergence analysis and numerical experiments are obtained for the special case.

A weak Galerkin generalized multiscale finite element method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Wang, Junping; Ye, Xiu

In this study, we propose a general framework for weak Galerkin generalized multiscale (WG-GMS) finite element method for the elliptic problems with rapidly oscillating or high contrast coefficients. This general WG-GMS method features in high order accuracy on general meshes and can work with multiscale basis derived by different numerical schemes. A special case is studied under this WG-GMS framework in which the multiscale basis functions are obtained by solving local problem with the weak Galerkin finite element method. Convergence analysis and numerical experiments are obtained for the special case.

Method and apparatus for testing a forward-moving strand

DOEpatents

Ducommun, Joel; Vulliens, Philippe

1980-01-01

In a method for testing a continuously forward-moving strand a light beam which passes along a plane that extends approximately perpendicularly to the longitudinal axis of the strand is introduced into the strand. The brightness value is measured on a place of the strand exterior which is distal from the light incidence place by means of at least one photoelectronic element disposed directly on the strand exterior and the measured result is evaluated in a gating circuit which is electrically connected to the photoelectronic element.

Domain decomposition methods for nonconforming finite element spaces of Lagrange-type

NASA Technical Reports Server (NTRS)

Cowsar, Lawrence C.

1993-01-01

In this article, we consider the application of three popular domain decomposition methods to Lagrange-type nonconforming finite element discretizations of scalar, self-adjoint, second order elliptic equations. The additive Schwarz method of Dryja and Widlund, the vertex space method of Smith, and the balancing method of Mandel applied to nonconforming elements are shown to converge at a rate no worse than their applications to the standard conforming piecewise linear Galerkin discretization. Essentially, the theory for the nonconforming elements is inherited from the existing theory for the conforming elements with only modest modification by constructing an isomorphism between the nonconforming finite element space and a space of continuous piecewise linear functions.

Application of finite elements heterogeneous multi-scale method to eddy currents non destructive testing of carbon composites material

NASA Astrophysics Data System (ADS)

Khebbab, Mohamed; Feliachi, Mouloud; El Hadi Latreche, Mohamed

2018-03-01

In this present paper, a simulation of eddy current non-destructive testing (EC NDT) on unidirectional carbon fiber reinforced polymer is performed; for this magneto-dynamic formulation in term of magnetic vector potential is solved using finite element heterogeneous multi-scale method (FE HMM). FE HMM has as goal to compute the homogenized solution without calculating the homogenized tensor explicitly, the solution is based only on the physical characteristic known in micro domain. This feature is well adapted to EC NDT to evaluate defect in carbon composite material in microscopic scale, where the defect detection is performed by coil impedance measurement; the measurement value is intimately linked to material characteristic in microscopic level. Based on this, our model can handle different defects such as: cracks, inclusion, internal electrical conductivity changes, heterogeneities, etc. The simulation results were compared with the solution obtained with homogenized material using mixture law, a good agreement was found.

Low Cost Nuclear Thermal Rocket Cermet Fuel Element Environment Testing

NASA Technical Reports Server (NTRS)

Bradley, David E.; Mireles, Omar R.; Hickman, Robert R.

2011-01-01

Deep space missions with large payloads require high specific impulse (Isp) and relatively high thrust in order to achieve mission goals in reasonable time frames. Conventional, storable propellants produce average Isp. Nuclear thermal rockets (NTR) capable of high Isp thrust have been proposed. NTR employs heat produced by fission reaction to heat and therefore accelerate hydrogen which is then forced through a rocket nozzle providing thrust. Fuel element temperatures are very high (up to 3000K) and hydrogen is highly reactive with most materials at high temperatures. Data covering the effects of high temperature hydrogen exposure on fuel elements is limited. The primary concern is the mechanical failure of fuel elements which employ high-melting-point metals, ceramics or a combination (cermet) as a structural matrix into which the nuclear fuel is distributed. It is not necessary to include fissile material in test samples intended to explore high temperature hydrogen exposure of the structural support matrices. A small-scale test bed designed to heat fuel element samples via non-contact RF heating and expose samples to hydrogen is being developed to assist in optimal material and manufacturing process selection without employing fissile material. This paper details the test bed design and results of testing conducted to date.

Introducing the Boundary Element Method with MATLAB

ERIC Educational Resources Information Center

Ang, Keng-Cheng

2008-01-01

The boundary element method provides an excellent platform for learning and teaching a computational method for solving problems in physical and engineering science. However, it is often left out in many undergraduate courses as its implementation is deemed to be difficult. This is partly due to the perception that coding the method requires…

Wave Scattering in Heterogeneous Media using the Finite Element Method

DTIC Science & Technology

2016-10-21

AFRL-AFOSR-JP-TR-2016-0086 Wave Scattering in Heterogeneous Media using the Finite Element Method Chiruvai Vendhan INDIAN INSTITUTE OF TECHNOLOGY...Scattering in Heterogeneous Media using the Finite Element Method 5a.  CONTRACT NUMBER 5b.  GRANT NUMBER FA2386-12-1-4026 5c.  PROGRAM ELEMENT NUMBER 61102F 6...14.  ABSTRACT The primary aim of this study is to develop a finite element model for elastic scattering by axisymmetric bodies submerged in a

Space-Time Conservation Element and Solution Element Method Being Developed

NASA Technical Reports Server (NTRS)

Chang, Sin-Chung; Himansu, Ananda; Jorgenson, Philip C. E.; Loh, Ching-Yuen; Wang, Xiao-Yen; Yu, Sheng-Tao

1999-01-01

The engineering research and design requirements of today pose great computer-simulation challenges to engineers and scientists who are called on to analyze phenomena in continuum mechanics. The future will bring even more daunting challenges, when increasingly complex phenomena must be analyzed with increased accuracy. Traditionally used numerical simulation methods have evolved to their present state by repeated incremental extensions to broaden their scope. They are reaching the limits of their applicability and will need to be radically revised, at the very least, to meet future simulation challenges. At the NASA Lewis Research Center, researchers have been developing a new numerical framework for solving conservation laws in continuum mechanics, namely, the Space-Time Conservation Element and Solution Element Method, or the CE/SE method. This method has been built from fundamentals and is not a modification of any previously existing method. It has been designed with generality, simplicity, robustness, and accuracy as cornerstones. The CE/SE method has thus far been applied in the fields of computational fluid dynamics, computational aeroacoustics, and computational electromagnetics. Computer programs based on the CE/SE method have been developed for calculating flows in one, two, and three spatial dimensions. Results have been obtained for numerous problems and phenomena, including various shock-tube problems, ZND detonation waves, an implosion and explosion problem, shocks over a forward-facing step, a blast wave discharging from a nozzle, various acoustic waves, and shock/acoustic-wave interactions. The method can clearly resolve shock/acoustic-wave interactions, wherein the difference of the magnitude between the acoustic wave and shock could be up to six orders. In two-dimensional flows, the reflected shock is as crisp as the leading shock. CE/SE schemes are currently being used for advanced applications to jet and fan noise prediction and to chemically

Experimental validation of boundary element methods for noise prediction

NASA Technical Reports Server (NTRS)

Seybert, A. F.; Oswald, Fred B.

1992-01-01

Experimental validation of methods to predict radiated noise is presented. A combined finite element and boundary element model was used to predict the vibration and noise of a rectangular box excited by a mechanical shaker. The predicted noise was compared to sound power measured by the acoustic intensity method. Inaccuracies in the finite element model shifted the resonance frequencies by about 5 percent. The predicted and measured sound power levels agree within about 2.5 dB. In a second experiment, measured vibration data was used with a boundary element model to predict noise radiation from the top of an operating gearbox. The predicted and measured sound power for the gearbox agree within about 3 dB.

Ultrasonic Method for Deployment Mechanism Bolt Element Preload Verification

NASA Technical Reports Server (NTRS)

Johnson, Eric C.; Kim, Yong M.; Morris, Fred A.; Mitchell, Joel; Pan, Robert B.

2014-01-01

Deployment mechanisms play a pivotal role in mission success. These mechanisms often incorporate bolt elements for which a preload within a specified range is essential for proper operation. A common practice is to torque these bolt elements to a specified value during installation. The resulting preload, however, can vary significantly with applied torque for a number of reasons. The goal of this effort was to investigate ultrasonic methods as an alternative for bolt preload verification in such deployment mechanisms. A family of non-explosive release mechanisms widely used by satellite manufacturers was chosen for the work. A willing contractor permitted measurements on a sampling of bolt elements for these release mechanisms that were installed by a technician following a standard practice. A variation of approximately 50% (+/- 25%) in the resultant preloads was observed. An alternative ultrasonic method to set the preloads was then developed and calibration data was accumulated. The method was demonstrated on bolt elements installed in a fixture instrumented with a calibrated load cell and designed to mimic production practice. The ultrasonic method yielded results within +/- 3% of the load cell reading. The contractor has since adopted the alternative method for its future production. Introduction

Parallel, adaptive finite element methods for conservation laws

NASA Technical Reports Server (NTRS)

Biswas, Rupak; Devine, Karen D.; Flaherty, Joseph E.

1994-01-01

We construct parallel finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. A posteriori estimates of spatial errors are obtained by a p-refinement technique using superconvergence at Radau points. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We compare results using different limiting schemes and demonstrate parallel efficiency through computations on an NCUBE/2 hypercube. We also present results using adaptive h- and p-refinement to reduce the computational cost of the method.

A fully-coupled discontinuous Galerkin spectral element method for two-phase flow in petroleum reservoirs

NASA Astrophysics Data System (ADS)

Taneja, Ankur; Higdon, Jonathan

2018-01-01

A high-order spectral element discontinuous Galerkin method is presented for simulating immiscible two-phase flow in petroleum reservoirs. The governing equations involve a coupled system of strongly nonlinear partial differential equations for the pressure and fluid saturation in the reservoir. A fully implicit method is used with a high-order accurate time integration using an implicit Rosenbrock method. Numerical tests give the first demonstration of high order hp spatial convergence results for multiphase flow in petroleum reservoirs with industry standard relative permeability models. High order convergence is shown formally for spectral elements with up to 8th order polynomials for both homogeneous and heterogeneous permeability fields. Numerical results are presented for multiphase fluid flow in heterogeneous reservoirs with complex geometric or geologic features using up to 11th order polynomials. Robust, stable simulations are presented for heterogeneous geologic features, including globally heterogeneous permeability fields, anisotropic permeability tensors, broad regions of low-permeability, high-permeability channels, thin shale barriers and thin high-permeability fractures. A major result of this paper is the demonstration that the resolution of the high order spectral element method may be exploited to achieve accurate results utilizing a simple cartesian mesh for non-conforming geological features. Eliminating the need to mesh to the boundaries of geological features greatly simplifies the workflow for petroleum engineers testing multiple scenarios in the face of uncertainty in the subsurface geology.

Correlation of finite element free vibration predictions using random vibration test data. M.S. Thesis - Cleveland State Univ.

NASA Technical Reports Server (NTRS)

Chambers, Jeffrey A.

1994-01-01

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Solution-adaptive finite element method in computational fracture mechanics

NASA Technical Reports Server (NTRS)

Min, J. B.; Bass, J. M.; Spradley, L. W.

1993-01-01

Some recent results obtained using solution-adaptive finite element method in linear elastic two-dimensional fracture mechanics problems are presented. The focus is on the basic issue of adaptive finite element method for validating the applications of new methodology to fracture mechanics problems by computing demonstration problems and comparing the stress intensity factors to analytical results.

A 3D finite element ALE method using an approximate Riemann solution

DOE PAGES

Chiravalle, V. P.; Morgan, N. R.

2016-08-09

Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problemmore » results are presented.« less

A 3D finite element ALE method using an approximate Riemann solution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chiravalle, V. P.; Morgan, N. R.

Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problemmore » results are presented.« less

Two-point method uncertainty during control and measurement of cylindrical element diameters

NASA Astrophysics Data System (ADS)

Glukhov, V. I.; Shalay, V. V.; Radev, H.

2018-04-01

The topic of the article is devoted to the urgent problem of the reliability of technical products geometric specifications measurements. The purpose of the article is to improve the quality of parts linear sizes control by the two-point measurement method. The article task is to investigate methodical extended uncertainties in measuring cylindrical element linear sizes. The investigation method is a geometric modeling of the element surfaces shape and location deviations in a rectangular coordinate system. The studies were carried out for elements of various service use, taking into account their informativeness, corresponding to the kinematic pairs classes in theoretical mechanics and the number of constrained degrees of freedom in the datum element function. Cylindrical elements with informativity of 4, 2, 1 and θ (zero) were investigated. The uncertainties estimation of in two-point measurements was made by comparing the results of of linear dimensions measurements with the functional diameters maximum and minimum of the element material. Methodical uncertainty is formed when cylindrical elements with maximum informativeness have shape deviations of the cut and the curvature types. Methodical uncertainty is formed by measuring the element average size for all types of shape deviations. The two-point measurement method cannot take into account the location deviations of a dimensional element, so its use for elements with informativeness less than the maximum creates unacceptable methodical uncertainties in measurements of the maximum, minimum and medium linear dimensions. Similar methodical uncertainties also exist in the arbitration control of the linear dimensions of the cylindrical elements by limiting two-point gauges.

Iterative methods for elliptic finite element equations on general meshes

NASA Technical Reports Server (NTRS)

Nicolaides, R. A.; Choudhury, Shenaz

1986-01-01

Iterative methods for arbitrary mesh discretizations of elliptic partial differential equations are surveyed. The methods discussed are preconditioned conjugate gradients, algebraic multigrid, deflated conjugate gradients, an element-by-element techniques, and domain decomposition. Computational results are included.

PWSCC Assessment by Using Extended Finite Element Method

NASA Astrophysics Data System (ADS)

Lee, Sung-Jun; Lee, Sang-Hwan; Chang, Yoon-Suk

2015-12-01

The head penetration nozzle of control rod driving mechanism (CRDM) is known to be susceptible to primary water stress corrosion cracking (PWSCC) due to the welding-induced residual stress. Especially, the J-groove dissimilar metal weld regions have received many attentions in the previous studies. However, even though several advanced techniques such as weight function and finite element alternating methods have been introduced to predict the occurrence of PWSCC, there are still difficulties in respect of applicability and efficiency. In this study, the extended finite element method (XFEM), which allows convenient crack element modeling by enriching degree of freedom (DOF) with special displacement function, was employed to evaluate structural integrity of the CRDM head penetration nozzle. The resulting stress intensity factors of surface cracks were verified for the reliability of proposed method through the comparison with those suggested in the American Society of Mechanical Engineering (ASME) code. The detailed results from the FE analyses are fully discussed in the manuscript.

Modeling of light dynamic cone penetration test - Panda 3 ® in granular material by using 3D Discrete element method

NASA Astrophysics Data System (ADS)

Tran, Quoc Anh; Chevalier, Bastien; Benz, Miguel; Breul, Pierre; Gourvès, Roland

2017-06-01

The recent technological developments made on the light dynamic penetration test Panda 3 ® provide a dynamic load-penetration curve σp - sp for each impact. This curve is influenced by the mechanical and physical properties of the investigated granular media. In order to analyze and exploit the load-penetration curve, a numerical model of penetration test using 3D Discrete Element Method is proposed for reproducing tests in dynamic conditions in granular media. All parameters of impact used in this model have at first been calibrated by respecting mechanical and geometrical properties of the hammer and the rod. There is a good agreement between experimental results and the ones obtained from simulations in 2D or 3D. After creating a sample, we will simulate the Panda 3 ®. It is possible to measure directly the dynamic load-penetration curve occurring at the tip for each impact. Using the force and acceleration measured in the top part of the rod, it is possible to separate the incident and reflected waves and then calculate the tip's load-penetration curve. The load-penetration curve obtained is qualitatively similar with that obtained by experimental tests. In addition, the frequency analysis of the measured signals present also a good compliance with that measured in reality when the tip resistance is qualitatively similar.

Feasibility of Combining Common Data Elements Across Studies to Test a Hypothesis.

PubMed

Corwin, Elizabeth J; Moore, Shirley M; Plotsky, Andrea; Heitkemper, Margaret M; Dorsey, Susan G; Waldrop-Valverde, Drenna; Bailey, Donald E; Docherty, Sharron L; Whitney, Joanne D; Musil, Carol M; Dougherty, Cynthia M; McCloskey, Donna J; Austin, Joan K; Grady, Patricia A

2017-05-01

The purpose of this article is to describe the outcomes of a collaborative initiative to share data across five schools of nursing in order to evaluate the feasibility of collecting common data elements (CDEs) and developing a common data repository to test hypotheses of interest to nursing scientists. This initiative extended work already completed by the National Institute of Nursing Research CDE Working Group that successfully identified CDEs related to symptoms and self-management, with the goal of supporting more complex, reproducible, and patient-focused research. Two exemplars describing the group's efforts are presented. The first highlights a pilot study wherein data sets from various studies by the represented schools were collected retrospectively, and merging of the CDEs was attempted. The second exemplar describes the methods and results of an initiative at one school that utilized a prospective design for the collection and merging of CDEs. Methods for identifying a common symptom to be studied across schools and for collecting the data dictionaries for the related data elements are presented for the first exemplar. The processes for defining and comparing the concepts and acceptable values, and for evaluating the potential to combine and compare the data elements are also described. Presented next are the steps undertaken in the second exemplar to prospectively identify CDEs and establish the data dictionaries. Methods for common measurement and analysis strategies are included. Findings from the first exemplar indicated that without plans in place a priori to ensure the ability to combine and compare data from disparate sources, doing so retrospectively may not be possible, and as a result hypothesis testing across studies may be prohibited. Findings from the second exemplar, however, indicated that a plan developed prospectively to combine and compare data sets is feasible and conducive to merged hypothesis testing. Although challenges exist in

Cylinder surface test with Chebyshev polynomial fitting method

NASA Astrophysics Data System (ADS)

Yu, Kui-bang; Guo, Pei-ji; Chen, Xi

2017-10-01

Zernike polynomials fitting method is often applied in the test of optical components and systems, used to represent the wavefront and surface error in circular domain. Zernike polynomials are not orthogonal in rectangular region which results in its unsuitable for the test of optical element with rectangular aperture such as cylinder surface. Applying the Chebyshev polynomials which are orthogonal among the rectangular area as an substitution to the fitting method, can solve the problem. Corresponding to a cylinder surface with diameter of 50 mm and F number of 1/7, a measuring system has been designed in Zemax based on Fizeau Interferometry. The expressions of the two-dimensional Chebyshev polynomials has been given and its relationship with the aberration has been presented. Furthermore, Chebyshev polynomials are used as base items to analyze the rectangular aperture test data. The coefficient of different items are obtained from the test data through the method of least squares. Comparing the Chebyshev spectrum in different misalignment, it show that each misalignment is independence and has a certain relationship with the certain Chebyshev terms. The simulation results show that, through the Legendre polynomials fitting method, it will be a great improvement in the efficient of the detection and adjustment of the cylinder surface test.

The L sub 1 finite element method for pure convection problems

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan

1991-01-01

The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

A fictitious domain approach for the Stokes problem based on the extended finite element method

NASA Astrophysics Data System (ADS)

Court, Sébastien; Fournié, Michel; Lozinski, Alexei

2014-01-01

In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by eXtended Finite Element Method and studied for Poisson problem in [Renard]. The method allows computations in domains whose boundaries do not match. A mixed finite element method is used for fluid flow. The interface between the fluid and the structure is localized by a level-set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf-sup condition between the spaces for velocity and the Lagrange multiplier. Convergence analysis is given and several numerical tests are performed to illustrate the capabilities of the method.

Dissipation-preserving spectral element method for damped seismic wave equations

NASA Astrophysics Data System (ADS)

Cai, Wenjun; Zhang, Huai; Wang, Yushun

2017-12-01

This article describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems, which has superior behaviors in long-time stability and dissipation preservation. To reveal the intrinsic dissipative properties of the model equations, we first reformulate the original systems in their equivalent conformal multi-symplectic structures and derive the corresponding conformal symplectic conservation laws. We thereafter separate each system into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed conformal symplectic method. A relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh wave in elastic wave propagation. More comprehensive numerical experiments are presented to investigate the long-time simulation, low dispersion and energy conservation properties of the conformal symplectic methods in both the attenuating

Ethnomathematics elements in Batik Bali using backpropagation method

NASA Astrophysics Data System (ADS)

Lestari, Mei; Irawan, Ari; Rahayu, Wanti; Wayan Parwati, Ni

2018-05-01

Batik is one of traditional arts that has been established by the UNESCO as Indonesia’s cultural heritage. Batik has varieties and motifs, and each motifs has its own uniqueness but seems similar, that makes it difficult to identify. This study aims to develop an application that can identify typical batik Bali with etnomatematics elements on it. Etnomatematics is a study that shows relation between culture and mathematics concepts. Etnomatematics in Batik Bali is more to geometrical concept in line of strong Balinese culture element. The identification process is use backpropagation method. Steps of backpropagation methods are image processing (including scalling and tresholding image process). Next step is insert the processed image to an artificial neural network. This study resulted an accuracy of identification of batik Bali that has Etnomatematics elements on it.

Low Cost Nuclear Thermal Rocket Cermet Fuel Element Environment Testing

NASA Technical Reports Server (NTRS)

Bradley, D. E.; Mireles, O. R.; Hickman, R. R.

2011-01-01

Deep space missions with large payloads require high specific impulse and relatively high thrust to achieve mission goals in reasonable time frames.1,2 Conventional storable propellants produce average specific impulse. Nuclear thermal rockets capable of producing high specific impulse are proposed. Nuclear thermal rockets employ heat produced by fission reaction to heat and therefore accelerate hydrogen, which is then forced through a rocket nozzle providing thrust. Fuel element temperatures are very high (up to 3000 K), and hydrogen is highly reactive with most materials at high temperatures. Data covering the effects of high-temperature hydrogen exposure on fuel elements are limited.3 The primary concern is the mechanical failure of fuel elements that employ high-melting-point metals, ceramics, or a combination (cermet) as a structural matrix into which the nuclear fuel is distributed. The purpose of the testing is to obtain data to assess the properties of the non-nuclear support materials, as-fabricated, and determine their ability to survive and maintain thermal performance in a prototypical NTR reactor environment of exposure to hydrogen at very high temperatures. The fission process of the planned fissile material and the resulting heating performance is well known and does not therefore require that active fissile material be integrated in this testing. A small-scale test bed designed to heat fuel element samples via non-contact radio frequency heating and expose samples to hydrogen is being developed to assist in optimal material and manufacturing process selection without employing fissile material. This paper details the test bed design and results of testing conducted to date.

Simulating Space Capsule Water Landing with Explicit Finite Element Method

NASA Technical Reports Server (NTRS)

Wang, John T.; Lyle, Karen H.

2007-01-01

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Global Sensitivity Applied to Dynamic Combined Finite Discrete Element Methods for Fracture Simulation

NASA Astrophysics Data System (ADS)

Godinez, H. C.; Rougier, E.; Osthus, D.; Srinivasan, G.

2017-12-01

Fracture propagation play a key role for a number of application of interest to the scientific community. From dynamic fracture processes like spall and fragmentation in metals and detection of gas flow in static fractures in rock and the subsurface, the dynamics of fracture propagation is important to various engineering and scientific disciplines. In this work we implement a global sensitivity analysis test to the Hybrid Optimization Software Suite (HOSS), a multi-physics software tool based on the combined finite-discrete element method, that is used to describe material deformation and failure (i.e., fracture and fragmentation) under a number of user-prescribed boundary conditions. We explore the sensitivity of HOSS for various model parameters that influence how fracture are propagated through a material of interest. The parameters control the softening curve that the model relies to determine fractures within each element in the mesh, as well a other internal parameters which influence fracture behavior. The sensitivity method we apply is the Fourier Amplitude Sensitivity Test (FAST), which is a global sensitivity method to explore how each parameter influence the model fracture and to determine the key model parameters that have the most impact on the model. We present several sensitivity experiments for different combination of model parameters and compare against experimental data for verification.

Extension and Validation of a Hybrid Particle-Finite Element Method for Hypervelocity Impact Simulation. Chapter 2

NASA Technical Reports Server (NTRS)

Fahrenthold, Eric P.; Shivarama, Ravishankar

2004-01-01

The hybrid particle-finite element method of Fahrenthold and Horban, developed for the simulation of hypervelocity impact problems, has been extended to include new formulations of the particle-element kinematics, additional constitutive models, and an improved numerical implementation. The extended formulation has been validated in three dimensional simulations of published impact experiments. The test cases demonstrate good agreement with experiment, good parallel speedup, and numerical convergence of the simulation results.

Comparison of Test and Finite Element Analysis for Two Full-Scale Helicopter Crash Tests

NASA Technical Reports Server (NTRS)

Annett, Martin S.; Horta,Lucas G.

2011-01-01

Finite element analyses have been performed for two full-scale crash tests of an MD-500 helicopter. The first crash test was conducted to evaluate the performance of a composite deployable energy absorber under combined flight loads. In the second crash test, the energy absorber was removed to establish the baseline loads. The use of an energy absorbing device reduced the impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to analytical results. Details of the full-scale crash tests and development of the system-integrated finite element model are briefly described along with direct comparisons of acceleration magnitudes and durations for the first full-scale crash test. Because load levels were significantly different between tests, models developed for the purposes of predicting the overall system response with external energy absorbers were not adequate under more severe conditions seen in the second crash test. Relative error comparisons were inadequate to guide model calibration. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used for the second full-scale crash test. The calibrated parameter set reduced 2-norm prediction error by 51% but did not improve impact shape orthogonality.

Using Response Surface Methods to Correlate the Modal Test of an Inflatable Test Article

NASA Technical Reports Server (NTRS)

Gupta, Anju

2013-01-01

This paper presents a practical application of response surface methods (RSM) to correlate a finite element model of a structural modal test. The test article is a quasi-cylindrical inflatable structure which primarily consists of a fabric weave, with an internal bladder and metallic bulkheads on either end. To mitigate model size, the fabric weave was simplified by representing it with shell elements. The task at hand is to represent the material behavior of the weave. The success of the model correlation is measured by comparing the four major modal frequencies of the analysis model to the four major modal frequencies of the test article. Given that only individual strap material properties were provided and material properties of the overall weave were not available, defining the material properties of the finite element model became very complex. First it was necessary to determine which material properties (modulus of elasticity in the hoop and longitudinal directions, shear modulus, Poisson's ratio, etc.) affected the modal frequencies. Then a Latin Hypercube of the parameter space was created to form an efficiently distributed finite case set. Each case was then analyzed with the results input into RSM. In the resulting response surface it was possible to see how each material parameter affected the modal frequencies of the analysis model. If the modal frequencies of the analysis model and its corresponding parameters match the test with acceptable accuracy, it can be said that the model correlation is successful.

Solving the MHD equations by the space time conservation element and solution element method

NASA Astrophysics Data System (ADS)

Zhang, Moujin; John Yu, S.-T.; Henry Lin, S.-C.; Chang, Sin-Chung; Blankson, Isaiah

2006-05-01

We apply the space-time conservation element and solution element (CESE) method to solve the ideal MHD equations with special emphasis on satisfying the divergence free constraint of magnetic field, i.e., ∇ · B = 0. In the setting of the CESE method, four approaches are employed: (i) the original CESE method without any additional treatment, (ii) a simple corrector procedure to update the spatial derivatives of magnetic field B after each time marching step to enforce ∇ · B = 0 at all mesh nodes, (iii) a constraint-transport method by using a special staggered mesh to calculate magnetic field B, and (iv) the projection method by solving a Poisson solver after each time marching step. To demonstrate the capabilities of these methods, two benchmark MHD flows are calculated: (i) a rotated one-dimensional MHD shock tube problem and (ii) a MHD vortex problem. The results show no differences between different approaches and all results compare favorably with previously reported data.

Method for recovering catalytic elements from fuel cell membrane electrode assemblies

DOEpatents

Shore, Lawrence [Edison, NJ; Matlin, Ramail [Berkeley Heights, NJ; Heinz, Robert [Ludwigshafen, DE

2012-06-26

A method for recovering catalytic elements from a fuel cell membrane electrode assembly is provided. The method includes converting the membrane electrode assembly into a particulate material, wetting the particulate material, forming a slurry comprising the wetted particulate material and an acid leachate adapted to dissolve at least one of the catalytic elements into a soluble catalytic element salt, separating the slurry into a depleted particulate material and a supernatant containing the catalytic element salt, and washing the depleted particulate material to remove any catalytic element salt retained within pores in the depleted particulate material.

A Discontinuous Galerkin Finite Element Method for Hamilton-Jacobi Equations

NASA Technical Reports Server (NTRS)

Hu, Changqing; Shu, Chi-Wang

1998-01-01

In this paper, we present a discontinuous Galerkin finite element method for solving the nonlinear Hamilton-Jacobi equations. This method is based on the Runge-Kutta discontinuous Galerkin finite element method for solving conservation laws. The method has the flexibility of treating complicated geometry by using arbitrary triangulation, can achieve high order accuracy with a local, compact stencil, and are suited for efficient parallel implementation. One and two dimensional numerical examples are given to illustrate the capability of the method.

Induction Heating Model of Cermet Fuel Element Environmental Test (CFEET)

NASA Technical Reports Server (NTRS)

Gomez, Carlos F.; Bradley, D. E.; Cavender, D. P.; Mireles, O. R.; Hickman, R. R.; Trent, D.; Stewart, E.

2013-01-01

Deep space missions with large payloads require high specific impulse and relatively high thrust to achieve mission goals in reasonable time frames. Nuclear Thermal Rockets (NTR) are capable of producing a high specific impulse by employing heat produced by a fission reactor to heat and therefore accelerate hydrogen through a rocket nozzle providing thrust. Fuel element temperatures are very high (up to 3000 K) and hydrogen is highly reactive with most materials at high temperatures. Data covering the effects of high-temperature hydrogen exposure on fuel elements are limited. The primary concern is the mechanical failure of fuel elements due to large thermal gradients; therefore, high-melting-point ceramics-metallic matrix composites (cermets) are one of the fuels under consideration as part of the Nuclear Cryogenic Propulsion Stage (NCPS) Advance Exploration System (AES) technology project at the Marshall Space Flight Center. The purpose of testing and analytical modeling is to determine their ability to survive and maintain thermal performance in a prototypical NTR reactor environment of exposure to hydrogen at very high temperatures and obtain data to assess the properties of the non-nuclear support materials. The fission process and the resulting heating performance are well known and do not require that active fissile material to be integrated in this testing. A small-scale test bed; Compact Fuel Element Environmental Tester (CFEET), designed to heat fuel element samples via induction heating and expose samples to hydrogen is being developed at MSFC to assist in optimal material and manufacturing process selection without utilizing fissile material. This paper details the analytical approach to help design and optimize the test bed using COMSOL Multiphysics for predicting thermal gradients induced by electromagnetic heating (Induction heating) and Thermal Desktop for radiation calculations.

Guided wave propagation and spectral element method for debonding damage assessment in RC structures

NASA Astrophysics Data System (ADS)

Wang, Ying; Zhu, Xinqun; Hao, Hong; Ou, Jinping

2009-07-01

A concrete-steel interface spectral element is developed to study the guided wave propagation along the steel rebar in the concrete. Scalar damage parameters characterizing changes in the interface (debonding damage) are incorporated into the formulation of the spectral finite element that is used for damage detection of reinforced concrete structures. Experimental tests are carried out on a reinforced concrete beam with embedded piezoelectric elements to verify the performance of the proposed model and algorithm. Parametric studies are performed to evaluate the effect of different damage scenarios on wave propagation in the reinforced concrete structures. Numerical simulations and experimental results show that the method is effective to model wave propagation along the steel rebar in concrete and promising to detect damage in the concrete-steel interface.

On finite element methods for the Helmholtz equation

NASA Technical Reports Server (NTRS)

Aziz, A. K.; Werschulz, A. G.

1979-01-01

The numerical solution of the Helmholtz equation is considered via finite element methods. A two-stage method which gives the same accuracy in the computed gradient as in the computed solution is discussed. Error estimates for the method using a newly developed proof are given, and the computational considerations which show this method to be computationally superior to previous methods are presented.

Simulation and evaluation of rupturable coated capsules by finite element method.

PubMed

Yang, Yan; Fang, Jie; Shen, Lian; Shan, Weiguang

2017-03-15

The objective of this study was to simulate and evaluate the burst behavior of rupturable coated capsules by finite element method (FEM). Film and coated capsules were prepared by dip-coating method and their dimensions were determined by stereomicroscope. Mechanical properties of the film were measured by tensile test and used as material properties of FEM models. Swelling pressure was determined by restrained expansion method and applied to the internal surface of FEM models. Water uptake of coated capsules was determined to study the formation of internal pressure. Burst test and in vitro dissolution was used to verify the FEM models, which were used to study and predict the coating burst behavior. Simulated results of coating burst behavior were well agreed with the experiment results. Swelling pressure, material properties and dimensions of coating had influence on the maximum stress. Burst pressure and critical L-HPC content were calculated for burst prediction and formulation optimization. FEM simulation was a feasible way to simulate and evaluate the burst behavior of coated capsules. Copyright © 2017 Elsevier B.V. All rights reserved.

Scattering Cross Section of Sound Waves by the Modal Element Method

NASA Technical Reports Server (NTRS)

Baumeister, Kenneth J.; Kreider, Kevin L.

1994-01-01

#he modal element method has been employed to determine the scattered field from a plane acoustic wave impinging on a two dimensional body. In the modal element method, the scattering body is represented by finite elements, which are coupled to an eigenfunction expansion representing the acoustic pressure in the infinite computational domain surrounding the body. The present paper extends the previous work by developing the algorithm necessary to calculate the acoustics scattering cross section by the modal element method. The scattering cross section is the acoustical equivalent to the Radar Cross Section (RCS) in electromagnetic theory. Since the scattering cross section is evaluated at infinite distance from the body, an asymptotic approximation is used in conjunction with the standard modal element method. For validation, the scattering cross section of the rigid circular cylinder is computed for the frequency range 0.1 is less than or equal to ka is less than or equal to 100. Results show excellent agreement with the analytic solution.

Battery element and method for making same

NASA Technical Reports Server (NTRS)

Clough, Thomas J. (Inventor); Pinsky, Naum (Inventor)

1989-01-01

In a method for producing a battery element useful as at least a positive plate in a lead-acid battery, the element comprising a fluid impervious, electrically conductive matrix having mutually opposing first and second surfaces and positive active electrode material associated with the first surface of the matrix, the improvement which comprises: conditioning the first surface to enhance the association of the positive active electrode material and the first surface; and applying and associating the positive active electrode material to the first surface.

A new method for measuring the rotational accuracy of rolling element bearings

NASA Astrophysics Data System (ADS)

Chen, Ye; Zhao, Xiangsong; Gao, Weiguo; Hu, Gaofeng; Zhang, Shizhen; Zhang, Dawei

2016-12-01

The rotational accuracy of a machine tool spindle has critical influence upon the geometric shape and surface roughness of finished workpiece. The rotational performance of the rolling element bearings is a main factor which affects the spindle accuracy, especially in the ultra-precision machining. In this paper, a new method is developed to measure the rotational accuracy of rolling element bearings of machine tool spindles. Variable and measurable axial preload is applied to seat the rolling elements in the bearing races, which is used to simulate the operating conditions. A high-precision (radial error is less than 300 nm) and high-stiffness (radial stiffness is 600 N/μm) hydrostatic reference spindle is adopted to rotate the inner race of the test bearing. To prevent the outer race from rotating, a 2-degrees of freedom flexure hinge mechanism (2-DOF FHM) is designed. Correction factors by using stiffness analysis are adopted to eliminate the influences of 2-DOF FHM in the radial direction. Two capacitive displacement sensors with nano-resolution (the highest resolution is 9 nm) are used to measure the radial error motion of the rolling element bearing, without separating the profile error as the traditional rotational accuracy metrology of the spindle. Finally, experimental measurements are performed at different spindle speeds (100-4000 rpm) and axial preloads (75-780 N). Synchronous and asynchronous error motion values are evaluated to demonstrate the feasibility and repeatability of the developed method and instrument.

On modelling three-dimensional piezoelectric smart structures with boundary spectral element method

NASA Astrophysics Data System (ADS)

Zou, Fangxin; Aliabadi, M. H.

2017-05-01

The computational efficiency of the boundary element method in elastodynamic analysis can be significantly improved by employing high-order spectral elements for boundary discretisation. In this work, for the first time, the so-called boundary spectral element method is utilised to formulate the piezoelectric smart structures that are widely used in structural health monitoring (SHM) applications. The resultant boundary spectral element formulation has been validated by the finite element method (FEM) and physical experiments. The new formulation has demonstrated a lower demand on computational resources and a higher numerical stability than commercial FEM packages. Comparing to the conventional boundary element formulation, a significant reduction in computational expenses has been achieved. In summary, the boundary spectral element formulation presented in this paper provides a highly efficient and stable mathematical tool for the development of SHM applications.

A New Expanded Mixed Element Method for Convection-Dominated Sobolev Equation

PubMed Central

Wang, Jinfeng; Li, Hong; Fang, Zhichao

2014-01-01

We propose and analyze a new expanded mixed element method, whose gradient belongs to the simple square integrable space instead of the classical H(div; Ω) space of Chen's expanded mixed element method. We study the new expanded mixed element method for convection-dominated Sobolev equation, prove the existence and uniqueness for finite element solution, and introduce a new expanded mixed projection. We derive the optimal a priori error estimates in L 2-norm for the scalar unknown u and a priori error estimates in (L 2)2-norm for its gradient λ and its flux σ. Moreover, we obtain the optimal a priori error estimates in H 1-norm for the scalar unknown u. Finally, we obtained some numerical results to illustrate efficiency of the new method. PMID:24701153

Developments in variational methods for high performance plate and shell elements

NASA Technical Reports Server (NTRS)

Felippa, Carlos A.; Militello, Carmelo

1991-01-01

High performance elements are simple finite elements constructed to deliver engineering accuracy with coarse arbitrary grids. This is part of a series on the variational foundations of high-performance elements, with emphasis on plate and shell elements constructed with the free formulation (FF) and assumed natural strain (ANS) methods. Parameterized variational principles are studied that provide a common foundation for the FF and ANS methods, as well as for a combination of both. From this unified formulation a variant of the ANS formulation, called the assumed natural deviatoric strain (ANDES) formulation, emerges as an important special case. The first ANDES element, a high-performance 9 degrees of freedom triangular Kirchhoff plate bending element, is briefly described to illustrate the use of the new formulation.

Criticality safety evaluation for the Advanced Test Reactor enhanced low enriched uranium fuel elements

DOE Office of Scientific and Technical Information (OSTI.GOV)

Montierth, Leland M.

2016-07-19

The Global Threat Reduction Initiative (GTRI) convert program is developing a high uranium density fuel based on a low enriched uranium (LEU) uranium-molybdenum alloy. Testing of prototypic GTRI fuel elements is necessary to demonstrate integrated fuel performance behavior and scale-up of fabrication techniques. GTRI Enhanced LEU Fuel (ELF) elements based on the ATR-Standard Size elements (all plates fueled) are to be fabricated for testing in the Advanced Test Reactor (ATR). While a specific ELF element design will eventually be provided for detailed analyses and in-core testing, this criticality safety evaluation (CSE) is intended to evaluate a hypothetical ELF element designmore » for criticality safety purposes. Existing criticality analyses have analyzed Standard (HEU) ATR elements from which controls have been derived. This CSE documents analysis that determines the reactivity of the hypothetical ELF fuel elements relative to HEU ATR elements and whether the existing HEU ATR element controls bound the ELF element. The initial calculations presented in this CSE analyzed the original ELF design, now referred to as Mod 0.1. In addition, as part of a fuel meat thickness optimization effort for reactor performance, other designs have been evaluated. As of early 2014 the most current conceptual designs are Mk1A and Mk1B, that were previously referred to as conceptual designs Mod 0.10 and Mod 0.11, respectively. Revision 1 evaluates the reactivity of the ATR HEU Mark IV elements for a comparison with the Mark VII elements.« less

A survey of mixed finite element methods

NASA Technical Reports Server (NTRS)

Brezzi, F.

1987-01-01

This paper is an introduction to and an overview of mixed finite element methods. It discusses the mixed formulation of certain basic problems in elasticity and hydrodynamics. It also discusses special techniques for solving the discrete problem.

High-order local maximum principle preserving (MPP) discontinuous Galerkin finite element method for the transport equation

NASA Astrophysics Data System (ADS)

Anderson, R.; Dobrev, V.; Kolev, Tz.; Kuzmin, D.; Quezada de Luna, M.; Rieben, R.; Tomov, V.

2017-04-01

In this work we present a FCT-like Maximum-Principle Preserving (MPP) method to solve the transport equation. We use high-order polynomial spaces; in particular, we consider up to 5th order spaces in two and three dimensions and 23rd order spaces in one dimension. The method combines the concepts of positive basis functions for discontinuous Galerkin finite element spatial discretization, locally defined solution bounds, element-based flux correction, and non-linear local mass redistribution. We consider a simple 1D problem with non-smooth initial data to explain and understand the behavior of different parts of the method. Convergence tests in space indicate that high-order accuracy is achieved. Numerical results from several benchmarks in two and three dimensions are also reported.

A multilevel correction adaptive finite element method for Kohn-Sham equation

NASA Astrophysics Data System (ADS)

Hu, Guanghui; Xie, Hehu; Xu, Fei

2018-02-01

In this paper, an adaptive finite element method is proposed for solving Kohn-Sham equation with the multilevel correction technique. In the method, the Kohn-Sham equation is solved on a fixed and appropriately coarse mesh with the finite element method in which the finite element space is kept improving by solving the derived boundary value problems on a series of adaptively and successively refined meshes. A main feature of the method is that solving large scale Kohn-Sham system is avoided effectively, and solving the derived boundary value problems can be handled efficiently by classical methods such as the multigrid method. Hence, the significant acceleration can be obtained on solving Kohn-Sham equation with the proposed multilevel correction technique. The performance of the method is examined by a variety of numerical experiments.

Static Buckling Model Tests and Elasto-plastic Finite Element Analysis of a Pile in Layers with Various Thicknesses

NASA Astrophysics Data System (ADS)

Okajima, Kenji; Imai, Junichi; Tanaka, Tadatsugu; Iida, Toshiaki

Damage to piles in the liquefied ground is frequently reported. Buckling by the excess vertical load could be one of the causes of the pile damage, as well as the lateral flow of the ground and the lateral load at the pile head. The buckling mechanism is described as a complicated interaction between the pile deformation by the vertical load and the earth pressure change cased by the pile deformation. In this study, series of static buckling model tests of a pile were carried out in dried sand ground with various thickness of the layer. Finite element analysis was applied to the test results to verify the effectiveness of the elasto-plastic finite element analysis combining the implicit-explicit mixed type dynamic relaxation method with the return mapping method to the pile buckling problems. The test results and the analysis indicated the possibility that the buckling load of a pile decreases greatly where the thickness of the layer increases.

Random element method for numerical modeling of diffusional processes

NASA Technical Reports Server (NTRS)

Ghoniem, A. F.; Oppenheim, A. K.

1982-01-01

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Analysis of Brick Masonry Wall using Applied Element Method

NASA Astrophysics Data System (ADS)

Lincy Christy, D.; Madhavan Pillai, T. M.; Nagarajan, Praveen

2018-03-01

The Applied Element Method (AEM) is a versatile tool for structural analysis. Analysis is done by discretising the structure as in the case of Finite Element Method (FEM). In AEM, elements are connected by a set of normal and shear springs instead of nodes. AEM is extensively used for the analysis of brittle materials. Brick masonry wall can be effectively analyzed in the frame of AEM. The composite nature of masonry wall can be easily modelled using springs. The brick springs and mortar springs are assumed to be connected in series. The brick masonry wall is analyzed and failure load is determined for different loading cases. The results were used to find the best aspect ratio of brick to strengthen brick masonry wall.

Structure-preserving spectral element method in attenuating seismic wave modeling

NASA Astrophysics Data System (ADS)

Cai, Wenjun; Zhang, Huai

2016-04-01

This work describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems which has superior behaviors in long-time stability and dissipation preservation. To construct the conformal symplectic method, we first reformulate the damped acoustic wave equation and the elastic wave equations in their equivalent conformal multi-symplectic structures, which naturally reveal the intrinsic properties of the original systems, especially, the dissipation laws. We thereafter separate each structures into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed numerical scheme, which is conformal symplectic and can therefore guarantee the numerical stability and dissipation preservation after a large time modeling. Additionally, a relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh-wave propagation. More comprehensive

The Spectral Element Method for Geophysical Flows

NASA Astrophysics Data System (ADS)

Taylor, Mark

1998-11-01

We will describe SEAM, a Spectral Element Atmospheric Model. SEAM solves the 3D primitive equations used in climate modeling and medium range forecasting. SEAM uses a spectral element discretization for the surface of the globe and finite differences in the vertical direction. The model is spectrally accurate, as demonstrated by a variety of test cases. It is well suited for modern distributed-shared memory computers, sustaining over 24 GFLOPS on a 240 processor HP Exemplar. This performance has allowed us to run several interesting simulations in full spherical geometry at high resolution (over 22 million grid points).

Improved accuracy for finite element structural analysis via an integrated force method

NASA Technical Reports Server (NTRS)

Patnaik, S. N.; Hopkins, D. A.; Aiello, R. A.; Berke, L.

1992-01-01

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

Discontinuous Galerkin finite element method for solving population density functions of cortical pyramidal and thalamic neuronal populations.

PubMed

Huang, Chih-Hsu; Lin, Chou-Ching K; Ju, Ming-Shaung

2015-02-01

Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons. Copyright © 2014 Elsevier Ltd. All rights reserved.

Finite-Element Methods for Real-Time Simulation of Surgery

NASA Technical Reports Server (NTRS)

Basdogan, Cagatay

2003-01-01

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Sensitivity of Particle Size in Discrete Element Method to Particle Gas Method (DEM_PGM) Coupling in Underbody Blast Simulations

DTIC Science & Technology

2016-06-12

Particle Size in Discrete Element Method to Particle Gas Method (DEM_PGM) Coupling in Underbody Blast Simulations Venkatesh Babu, Kumar Kulkarni, Sanjay...buried in soil viz., (1) coupled discrete element & particle gas methods (DEM-PGM) and (2) Arbitrary Lagrangian-Eulerian (ALE), are investigated. The...DEM_PGM and identify the limitations/strengths compared to the ALE method. Discrete Element Method (DEM) can model individual particle directly, and

Experimental Optimization Methods for Multi-Element Airfoils

NASA Technical Reports Server (NTRS)

Landman, Drew; Britcher, Colin P.

1996-01-01

A modern three element airfoil model with a remotely activated flap was used to investigate optimum flap testing position using an automated optimization algorithm in wind tunnel tests. Detailed results for lift coefficient versus flap vertical and horizontal position are presented for two angles of attack: 8 and 14 degrees. An on-line first order optimizer is demonstrated which automatically seeks the optimum lift as a function of flap position. Future work with off-line optimization techniques is introduced and aerodynamic hysteresis effects due to flap movement with flow on are discussed.

Design Methods for Load-bearing Elements from Crosslaminated Timber

NASA Astrophysics Data System (ADS)

Vilguts, A.; Serdjuks, D.; Goremikins, V.

2015-11-01

Cross-laminated timber is an environmentally friendly material, which possesses a decreased level of anisotropy in comparison with the solid and glued timber. Cross-laminated timber could be used for load-bearing walls and slabs of multi-storey timber buildings as well as decking structures of pedestrian and road bridges. Design methods of cross-laminated timber elements subjected to bending and compression with bending were considered. The presented methods were experimentally validated and verified by FEM. Two cross-laminated timber slabs were tested at the action of static load. Pine wood was chosen as a board's material. Freely supported beam with the span equal to 1.9 m, which was loaded by the uniformly distributed load, was a design scheme of the considered plates. The width of the plates was equal to 1 m. The considered cross-laminated timber plates were analysed by FEM method. The comparison of stresses acting in the edge fibres of the plate and the maximum vertical displacements shows that both considered methods can be used for engineering calculations. The difference between the results obtained experimentally and analytically is within the limits from 2 to 31%. The difference in results obtained by effective strength and stiffness and transformed sections methods was not significant.

Hybrid finite difference/finite element immersed boundary method.

PubMed

E Griffith, Boyce; Luo, Xiaoyu

2017-12-01

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach uses a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach. © 2017 The Authors International  Journal  for  Numerical  Methods  in  Biomedical  Engineering Published by John Wiley & Sons Ltd.

Three dimensional finite element methods: Their role in the design of DC accelerator systems

NASA Astrophysics Data System (ADS)

Podaru, Nicolae C.; Gottdang, A.; Mous, D. J. W.

2013-04-01

High Voltage Engineering has designed, built and tested a 2 MV dual irradiation system that will be applied for radiation damage studies and ion beam material modification. The system consists of two independent accelerators which support simultaneous proton and electron irradiation (energy range 100 keV - 2 MeV) of target sizes of up to 300 × 300 mm2. Three dimensional finite element methods were used in the design of various parts of the system. The electrostatic solver was used to quantify essential parameters of the solid-state power supply generating the DC high voltage. The magnetostatic solver and ray tracing were used to optimize the electron/ion beam transport. Close agreement between design and measurements of the accelerator characteristics as well as beam performance indicate the usefulness of three dimensional finite element methods during accelerator system design.

Modeling and Simulation of a Nuclear Fuel Element Test Section

NASA Technical Reports Server (NTRS)

Moran, Robert P.; Emrich, William

2011-01-01

"The Nuclear Thermal Rocket Element Environmental Simulator" test section closely simulates the internal operating conditions of a thermal nuclear rocket. The purpose of testing is to determine the ideal fuel rod characteristics for optimum thermal heat transfer to their hydrogen cooling/working fluid while still maintaining fuel rod structural integrity. Working fluid exhaust temperatures of up to 5,000 degrees Fahrenheit can be encountered. The exhaust gas is rendered inert and massively reduced in temperature for analysis using a combination of water cooling channels and cool N2 gas injectors in the H2-N2 mixer portion of the test section. An extensive thermal fluid analysis was performed in support of the engineering design of the H2-N2 mixer in order to determine the maximum "mass flow rate"-"operating temperature" curve of the fuel elements hydrogen exhaust gas based on the test facilities available cooling N2 mass flow rate as the limiting factor.

Elemental Water Impact Test: Phase 1 20-Inch Hemisphere

NASA Technical Reports Server (NTRS)

Vassilakos, Gregory J.

2015-01-01

Spacecraft are being designed based on LS-DYNA simulations of water landing impacts. The Elemental Water Impact Test (EWIT) series was undertaken to assess the accuracy of LS-DYNA water impact simulations. Phase 1 of the EWIT series featured water impact tests of a 20-inch hemisphere dropped from heights of 5 feet and 10 feet. The hemisphere was outfitted with an accelerometer and three pressure gages. The focus of this report is the correlation of analytical models against test data.

Verification of a non-hydrostatic dynamical core using horizontally spectral element vertically finite difference method: 2-D aspects

NASA Astrophysics Data System (ADS)

Choi, S.-J.; Giraldo, F. X.; Kim, J.; Shin, S.

2014-06-01

The non-hydrostatic (NH) compressible Euler equations of dry atmosphere are solved in a simplified two dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and a finite difference method (FDM) for the vertical discretization. The SEM uses high-order nodal basis functions associated with Lagrange polynomials based on Gauss-Lobatto-Legendre (GLL) quadrature points. The FDM employs a third-order upwind biased scheme for the vertical flux terms and a centered finite difference scheme for the vertical derivative terms and quadrature. The Euler equations used here are in a flux form based on the hydrostatic pressure vertical coordinate, which are the same as those used in the Weather Research and Forecasting (WRF) model, but a hybrid sigma-pressure vertical coordinate is implemented in this model. We verified the model by conducting widely used standard benchmark tests: the inertia-gravity wave, rising thermal bubble, density current wave, and linear hydrostatic mountain wave. The results from those tests demonstrate that the horizontally spectral element vertically finite difference model is accurate and robust. By using the 2-D slice model, we effectively show that the combined spatial discretization method of the spectral element and finite difference method in the horizontal and vertical directions, respectively, offers a viable method for the development of a NH dynamical core.

A class of hybrid finite element methods for electromagnetics: A review

NASA Technical Reports Server (NTRS)

Volakis, J. L.; Chatterjee, A.; Gong, J.

1993-01-01

Integral equation methods have generally been the workhorse for antenna and scattering computations. In the case of antennas, they continue to be the prominent computational approach, but for scattering applications the requirement for large-scale computations has turned researchers' attention to near neighbor methods such as the finite element method, which has low O(N) storage requirements and is readily adaptable in modeling complex geometrical features and material inhomogeneities. In this paper, we review three hybrid finite element methods for simulating composite scatterers, conformal microstrip antennas, and finite periodic arrays. Specifically, we discuss the finite element method and its application to electromagnetic problems when combined with the boundary integral, absorbing boundary conditions, and artificial absorbers for terminating the mesh. Particular attention is given to large-scale simulations, methods, and solvers for achieving low memory requirements and code performance on parallel computing architectures.

A boundary element alternating method for two-dimensional mixed-mode fracture problems

NASA Technical Reports Server (NTRS)

Raju, I. S.; Krishnamurthy, T.

1992-01-01

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Methods and devices for fabricating and assembling printable semiconductor elements

DOEpatents

Nuzzo, Ralph G; Rogers, John A; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao

2014-03-04

The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

Methods and devices for fabricating and assembling printable semiconductor elements

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nuzzo, Ralph G.; Rogers, John A.; Menard, Etienne

The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

Discrete Element Method Simulation of a Boulder Extraction From an Asteroid

NASA Technical Reports Server (NTRS)

Kulchitsky, Anton K.; Johnson, Jerome B.; Reeves, David M.; Wilkinson, Allen

2014-01-01

The force required to pull 7t and 40t polyhedral boulders from the surface of an asteroid is simulated using the discrete element method considering the effects of microgravity, regolith cohesion and boulder acceleration. The connection between particle surface energy and regolith cohesion is estimated by simulating a cohesion sample tearing test. An optimal constant acceleration is found where the peak net force from inertia and cohesion is a minimum. Peak pulling forces can be further reduced by using linear and quadratic acceleration functions with up to a 40% reduction in force for quadratic acceleration.

A Moving Discontinuous Galerkin Finite Element Method for Flows with Interfaces

DTIC Science & Technology

2017-12-07

Naval Research Laboratory Washington, DC 20375-5320 NRL/MR/6040--17-9765 A Moving Discontinuous Galerkin Finite Element Method for Flows with...guidance to revise the method to ensure such properties. Acknowledgements This work was sponsored by the Office of Naval Research through the Naval...18. NUMBER OF PAGES 17. LIMITATION OF ABSTRACT A Moving Discontinuous Galerkin Finite Element Method for Flows with Interfaces Andrew Corrigan, Andrew

METHOD OF TESTING FOR LEAKS

DOEpatents

Creutz, E.C.; McAdams, Wm.A.; Foss, M.H.

1958-07-22

A method is described for detecting minute holes In fuel element jackets. The method comprises submerging the jacketed body in an atmosphere of a radioactive gas under pressure, the radioactive emanations from said gas being sufficientiy penetratIng to penetrate the jacket of the jacketed body. After the jacketed body is removed from the radtoactive gas atmosphere, it is exannined for the presence of emanations from radioactive gas which entered the jacketed body through the minute holes. In this manner, the detectton of radioactive emanations is a positive indication that the fuel element is not perfectly sealed.

Parallelized Three-Dimensional Resistivity Inversion Using Finite Elements And Adjoint State Methods

NASA Astrophysics Data System (ADS)

Schaa, Ralf; Gross, Lutz; Du Plessis, Jaco

2015-04-01

resistivity. The Hessian of the regularization term is used as preconditioner which requires an additional PDE solution in each iteration step. As it turns out, the relevant PDEs are naturally formulated in the finite element framework. Using the domain decomposition method provided in Escript, the inversion scheme has been parallelized for distributed memory computers with multi-core shared memory nodes. We show numerical examples from simple layered models to complex 3D models and compare with the results from other methods. The inversion scheme is furthermore tested on a field data example to characterise localised freshwater discharge in a coastal environment.. References: L. Gross and C. Kemp (2013) Large Scale Joint Inversion of Geophysical Data using the Finite Element Method in escript. ASEG Extended Abstracts 2013, http://dx.doi.org/10.1071/ASEG2013ab306

Stabilization of time domain acoustic boundary element method for the exterior problem avoiding the nonuniqueness.

PubMed

Jang, Hae-Won; Ih, Jeong-Guon

2013-03-01

The time domain boundary element method (TBEM) to calculate the exterior sound field using the Kirchhoff integral has difficulties in non-uniqueness and exponential divergence. In this work, a method to stabilize TBEM calculation for the exterior problem is suggested. The time domain CHIEF (Combined Helmholtz Integral Equation Formulation) method is newly formulated to suppress low order fictitious internal modes. This method constrains the surface Kirchhoff integral by forcing the pressures at the additional interior points to be zero when the shortest retarded time between boundary nodes and an interior point elapses. However, even after using the CHIEF method, the TBEM calculation suffers the exponential divergence due to the remaining unstable high order fictitious modes at frequencies higher than the frequency limit of the boundary element model. For complete stabilization, such troublesome modes are selectively adjusted by projecting the time response onto the eigenspace. In a test example for a transiently pulsating sphere, the final average error norm of the stabilized response compared to the analytic solution is 2.5%.

Method of securing filter elements

DOE Office of Scientific and Technical Information (OSTI.GOV)

Brown, Erik P.; Haslam, Jeffery L.; Mitchell, Mark A.

2016-10-04

A filter securing system including a filter unit body housing; at least one tubular filter element positioned in the filter unit body housing, the tubular filter element having a closed top and an open bottom; a dimple in either the filter unit body housing or the top of the tubular filter element; and a socket in either the filter unit body housing or the top of the tubular filter element that receives the dimple in either the filter unit body housing or the top of the tubular filter element to secure the tubular filter element to the filter unit bodymore » housing.« less

Axisymmetric solid elements by a rational hybrid stress method

NASA Technical Reports Server (NTRS)

Tian, Z.; Pian, T. H. H.

1985-01-01

Four-node axisymmetric solid elements are derived by a new version of hybrid method for which the assumed stresses are expressed in complete polynomials in natural coordinates. The stress equilibrium conditions are introduced through the use of additional displacements as Lagrange multipliers. A rational procedure is to choose the displacement terms such that the resulting strains are also of complete polynomials of the same order. Example problems all indicate that elements obtained by this procedure lead to better results in displacements and stresses than that by other finite elements.

Development and Application of the p-version of the Finite Element Method.

DTIC Science & Technology

1985-11-21

this property hierarchic families of finite elements. The h-version of the finite element method has been the subject of inten- sive study since the...early 1950’s and perhaps even earlier. Study of the p-version of the finite element method, on the other hand, began at Washington University in St...Louis in the early 1970’s and led to a more recent study of * .the h-p version. Research in the p-version (formerly called The Constraint Method) has

Analysis of Tire Tractive Performance on Deformable Terrain by Finite Element-Discrete Element Method

NASA Astrophysics Data System (ADS)

Nakashima, Hiroshi; Takatsu, Yuzuru

The goal of this study is to develop a practical and fast simulation tool for soil-tire interaction analysis, where finite element method (FEM) and discrete element method (DEM) are coupled together, and which can be realized on a desktop PC. We have extended our formerly proposed dynamic FE-DE method (FE-DEM) to include practical soil-tire system interaction, where not only the vertical sinkage of a tire, but also the travel of a driven tire was considered. Numerical simulation by FE-DEM is stable, and the relationships between variables, such as load-sinkage and sinkage-travel distance, and the gross tractive effort and running resistance characteristics, are obtained. Moreover, the simulation result is accurate enough to predict the maximum drawbar pull for a given tire, once the appropriate parameter values are provided. Therefore, the developed FE-DEM program can be applied with sufficient accuracy to interaction problems in soil-tire systems.

A Least-Squares Finite Element Method for Electromagnetic Scattering Problems

NASA Technical Reports Server (NTRS)

Wu, Jie; Jiang, Bo-nan

1996-01-01

The least-squares finite element method (LSFEM) is applied to electromagnetic scattering and radar cross section (RCS) calculations. In contrast to most existing numerical approaches, in which divergence-free constraints are omitted, the LSFF-M directly incorporates two divergence equations in the discretization process. The importance of including the divergence equations is demonstrated by showing that otherwise spurious solutions with large divergence occur near the scatterers. The LSFEM is based on unstructured grids and possesses full flexibility in handling complex geometry and local refinement Moreover, the LSFEM does not require any special handling, such as upwinding, staggered grids, artificial dissipation, flux-differencing, etc. Implicit time discretization is used and the scheme is unconditionally stable. By using a matrix-free iterative method, the computational cost and memory requirement for the present scheme is competitive with other approaches. The accuracy of the LSFEM is verified by several benchmark test problems.

Application of the Spectral Element Method to Acoustic Radiation

NASA Technical Reports Server (NTRS)

Doyle, James F.; Rizzi, Stephen A. (Technical Monitor)

2000-01-01

This report summarizes research to develop a capability for analysis of interior noise in enclosed structures when acoustically excited by an external random source. Of particular interest was the application to the study of noise and vibration transmission in thin-walled structures as typified by aircraft fuselages. Three related topics are focused upon. The first concerns the development of a curved frame spectral element, the second shows how the spectral element method for wave propagation in folded plate structures is extended to problems involving curved segmented plates. These are of significance because by combining these curved spectral elements with previously presented flat spectral elements, the dynamic response of geometrically complex structures can be determined. The third topic shows how spectral elements, which incorporate the effect of fluid loading on the structure, are developed for analyzing acoustic radiation from dynamically loaded extended plates.

Tree stability under wind: simulating uprooting with root breakage using a finite element method

PubMed Central

Yang, Ming; Défossez, Pauline; Danjon, Frédéric; Fourcaud, Thierry

2014-01-01

Background and Aims Windstorms are the major natural hazard affecting European forests, causing tree damage and timber losses. Modelling tree anchorage mechanisms has progressed with advances in plant architectural modelling, but it is still limited in terms of estimation of anchorage strength. This paper aims to provide a new model for root anchorage, including the successive breakage of roots during uprooting. Methods The model was based on the finite element method. The breakage of individual roots was taken into account using a failure law derived from previous work carried out on fibre metal laminates. Soil mechanical plasticity was considered using the Mohr–Coulomb failure criterion. The mechanical model for roots was implemented in the numerical code ABAQUS using beam elements embedded in a soil block meshed with 3-D solid elements. The model was tested by simulating tree-pulling experiments previously carried out on a tree of Pinus pinaster (maritime pine). Soil mechanical parameters were obtained from laboratory tests. Root system architecture was digitized and imported into ABAQUS while root material properties were estimated from the literature. Key Results Numerical simulations of tree-pulling tests exhibited realistic successive root breakages during uprooting, which could be seen in the resulting response curves. Broken roots could be visually located within the root system at any stage of the simulations. The model allowed estimation of anchorage strength in terms of the critical turning moment and accumulated energy, which were in good agreement with in situ measurements. Conclusions This study provides the first model of tree anchorage strength for P. pinaster derived from the mechanical strength of individual roots. The generic nature of the model permits its further application to other tree species and soil conditions. PMID:25006178

Method of holding optical elements without deformation during their fabrication

DOEpatents

Hed, P. Paul

1997-01-01

An improved method for securing and removing an optical element to and from a blocking tool without causing deformation of the optical element. A lens tissue is placed on the top surface of the blocking tool. Dots of UV cement are applied to the lens tissue without any of the dots contacting each other. An optical element is placed on top of the blocking tool with the lens tissue sandwiched therebetween. The UV cement is then cured. After subsequent fabrication steps, the bonded blocking tool, lens tissue, and optical element are placed in a debonding solution to soften the UV cement. The optical element is then removed from the blocking tool.

Spectral element method for elastic and acoustic waves in frequency domain

DOE Office of Scientific and Technical Information (OSTI.GOV)

Shi, Linlin; Zhou, Yuanguo; Wang, Jia-Min

Numerical techniques in time domain are widespread in seismic and acoustic modeling. In some applications, however, frequency-domain techniques can be advantageous over the time-domain approach when narrow band results are desired, especially if multiple sources can be handled more conveniently in the frequency domain. Moreover, the medium attenuation effects can be more accurately and conveniently modeled in the frequency domain. In this paper, we present a spectral-element method (SEM) in frequency domain to simulate elastic and acoustic waves in anisotropic, heterogeneous, and lossy media. The SEM is based upon the finite-element framework and has exponential convergence because of the usemore » of GLL basis functions. The anisotropic perfectly matched layer is employed to truncate the boundary for unbounded problems. Compared with the conventional finite-element method, the number of unknowns in the SEM is significantly reduced, and higher order accuracy is obtained due to its spectral accuracy. To account for the acoustic-solid interaction, the domain decomposition method (DDM) based upon the discontinuous Galerkin spectral-element method is proposed. Numerical experiments show the proposed method can be an efficient alternative for accurate calculation of elastic and acoustic waves in frequency domain.« less

Prediction of shear critical behavior of high-strength reinforced concrete columns using finite element methods

NASA Astrophysics Data System (ADS)

Alrasyid, Harun; Safi, Fahrudin; Iranata, Data; Chen-Ou, Yu

2017-11-01

This research shows the prediction of shear behavior of High-Strength Reinforced Concrete Columns using Finite-Element Method. The experimental data of nine half scale high-strength reinforced concrete were selected. These columns using specified concrete compressive strength of 70 MPa, specified yield strength of longitudinal and transverse reinforcement of 685 and 785 MPa, respectively. The VecTor2 finite element software was used to simulate the shear critical behavior of these columns. The combination axial compression load and monotonic loading were applied at this prediction. It is demonstrated that VecTor2 finite element software provides accurate prediction of load-deflection up to peak at applied load, but provide similar behavior at post peak load. The shear strength prediction provide by VecTor 2 are slightly conservative compare to test result.

Generalization of mixed multiscale finite element methods with applications

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, C S

Many science and engineering problems exhibit scale disparity and high contrast. The small scale features cannot be omitted in the physical models because they can affect the macroscopic behavior of the problems. However, resolving all the scales in these problems can be prohibitively expensive. As a consequence, some types of model reduction techniques are required to design efficient solution algorithms. For practical purpose, we are interested in mixed finite element problems as they produce solutions with certain conservative properties. Existing multiscale methods for such problems include the mixed multiscale finite element methods. We show that for complicated problems, the mixedmore » multiscale finite element methods may not be able to produce reliable approximations. This motivates the need of enrichment for coarse spaces. Two enrichment approaches are proposed, one is based on generalized multiscale finte element metthods (GMsFEM), while the other is based on spectral element-based algebraic multigrid (rAMGe). The former one, which is called mixed GMsFEM, is developed for both Darcy’s flow and linear elasticity. Application of the algorithm in two-phase flow simulations are demonstrated. For linear elasticity, the algorithm is subtly modified due to the symmetry requirement of the stress tensor. The latter enrichment approach is based on rAMGe. The algorithm differs from GMsFEM in that both of the velocity and pressure spaces are coarsened. Due the multigrid nature of the algorithm, recursive application is available, which results in an efficient multilevel construction of the coarse spaces. Stability, convergence analysis, and exhaustive numerical experiments are carried out to validate the proposed enrichment approaches. iii« less

Modelling of single walled carbon nanotube cylindrical structures with finite element method simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Günay, E.

In this study, the modulus of elasticity and shear modulus values of single-walled carbon nanotubes SWCNTs were modelled by using both finite element method and the Matlab code. Initially, cylindrical armchair and zigzag single walled 3D space frames were demonstrated as carbon nanostructures. Thereafter, macro programs were written by the Matlab code producing the space truss for zigzag and armchair models. 3D space frames were introduced to the ANSYS software and then tension, compression and additionally torsion tests were performed on zigzag and armchair carbon nanotubes with BEAM4 element in obtaining the exact values of elastic and shear modulus values.more » In this study, two different boundary conditions were tested and especially used in torsion loading. The equivalent shear modulus data was found by averaging the corresponding values obtained from ten different nodal points on the nanotube path. Finally, in this study it was determined that the elastic constant values showed proportional changes by increasing the carbon nanotube diameters up to a certain level but beyond this level these values remained stable.« less

A new method for constructing analytic elements for groundwater flow.

NASA Astrophysics Data System (ADS)

Strack, O. D.

2007-12-01

The analytic element method is based upon the superposition of analytic functions that are defined throughout the infinite domain, and can be used to meet a variety of boundary conditions. Analytic elements have been use successfully for a number of problems, mainly dealing with the Poisson equation (see, e.g., Theory and Applications of the Analytic Element Method, Reviews of Geophysics, 41,2/1005 2003 by O.D.L. Strack). The majority of these analytic elements consists of functions that exhibit jumps along lines or curves. Such linear analytic elements have been developed also for other partial differential equations, e.g., the modified Helmholz equation and the heat equation, and were constructed by integrating elementary solutions, the point sink and the point doublet, along a line. This approach is limiting for two reasons. First, the existence is required of the elementary solutions, and, second, the integration tends to limit the range of solutions that can be obtained. We present a procedure for generating analytic elements that requires merely the existence of a harmonic function with the desired properties; such functions exist in abundance. The procedure to be presented is used to generalize this harmonic function in such a way that the resulting expression satisfies the applicable differential equation. The approach will be applied, along with numerical examples, for the modified Helmholz equation and for the heat equation, while it is noted that the method is in no way restricted to these equations. The procedure is carried out entirely in terms of complex variables, using Wirtinger calculus.

A blended continuous–discontinuous finite element method for solving the multi-fluid plasma model

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sousa, E.M., E-mail: sousae@uw.edu; Shumlak, U., E-mail: shumlak@uw.edu

The multi-fluid plasma model represents electrons, multiple ion species, and multiple neutral species as separate fluids that interact through short-range collisions and long-range electromagnetic fields. The model spans a large range of temporal and spatial scales, which renders the model stiff and presents numerical challenges. To address the large range of timescales, a blended continuous and discontinuous Galerkin method is proposed, where the massive ion and neutral species are modeled using an explicit discontinuous Galerkin method while the electrons and electromagnetic fields are modeled using an implicit continuous Galerkin method. This approach is able to capture large-gradient ion and neutralmore » physics like shock formation, while resolving high-frequency electron dynamics in a computationally efficient manner. The details of the Blended Finite Element Method (BFEM) are presented. The numerical method is benchmarked for accuracy and tested using two-fluid one-dimensional soliton problem and electromagnetic shock problem. The results are compared to conventional finite volume and finite element methods, and demonstrate that the BFEM is particularly effective in resolving physics in stiff problems involving realistic physical parameters, including realistic electron mass and speed of light. The benefit is illustrated by computing a three-fluid plasma application that demonstrates species separation in multi-component plasmas.« less

Scientific use of the finite element method in Orthodontics

PubMed Central

Knop, Luegya; Gandini, Luiz Gonzaga; Shintcovsk, Ricardo Lima; Gandini, Marcia Regina Elisa Aparecida Schiavon

2015-01-01

INTRODUCTION: The finite element method (FEM) is an engineering resource applied to calculate the stress and deformation of complex structures, and has been widely used in orthodontic research. With the advantage of being a non-invasive and accurate method that provides quantitative and detailed data on the physiological reactions possible to occur in tissues, applying the FEM can anticipate the visualization of these tissue responses through the observation of areas of stress created from applied orthodontic mechanics. OBJECTIVE: This article aims at reviewing and discussing the stages of the finite element method application and its applicability in Orthodontics. RESULTS: FEM is able to evaluate the stress distribution at the interface between periodontal ligament and alveolar bone, and the shifting trend in various types of tooth movement when using different types of orthodontic devices. Therefore, it is necessary to know specific software for this purpose. CONCLUSIONS: FEM is an important experimental method to answer questions about tooth movement, overcoming the disadvantages of other experimental methods. PMID:25992996

Simply scan--optical methods for elemental carbon measurement in diesel exhaust particulate.

PubMed

Forder, James A

2014-08-01

This article describes a performance assessment of three optical methods, a Magee Scientific OT21 Transmissometer, a Hach-Lange Microcolor II difference gloss meter, and a combination of an office scanner with Adobe Photoshop software. The optical methods measure filter staining as a proxy for elemental carbon in diesel exhaust particulate (DEP) exposure assessment and the suitability of each as a replacement for the existing Bosch meter optical method. Filters loaded with DEP were produced from air in a non-coal mine and the exhaust gases from a mobile crane. These were measured with each apparatus and then by combustion to obtain a reference elemental carbon value. The results from each apparatus were then plotted against both the Bosch number and reference elemental carbon values. The equations of the best fit lines for these plots were derived, and these gave functions for elemental carbon and Bosch number from the output of each new optical method. For each optical method, the range of DEP loadings which can be measured has been determined, and conversion equations for elemental carbon and Bosch number have been obtained. All three optical methods studied will effectively quantify blackness as a measure of elemental carbon. Of these the Magee Scientific OT21 transmissometer has the best performance. The Microcolor II and scanner/photoshop methods will in addition allow conversion to Bosch number which may be useful if historical Bosch data are available and functions for this are described. The scanner/photoshop method demonstrates a technique to obtain measurements of DEP exposure without the need to purchase specialized instrumentation. © The Author 2014. Published by Oxford University Press on behalf of the British Occupational Hygiene Society.

Base excitation testing system using spring elements to pivotally mount wind turbine blades

DOEpatents

Cotrell, Jason; Hughes, Scott; Butterfield, Sandy; Lambert, Scott

2013-12-10

A system (1100) for fatigue testing wind turbine blades (1102) through forced or resonant excitation of the base (1104) of a blade (1102). The system (1100) includes a test stand (1112) and a restoring spring assembly (1120) mounted on the test stand (1112). The restoring spring assembly (1120) includes a primary spring element (1124) that extends outward from the test stand (1112) to a blade mounting plate (1130) configured to receive a base (1104) of blade (1102). During fatigue testing, a supported base (1104) of a blad (1102) may be pivotally mounted to the test stand (1112) via the restoring spring assembly (1120). The system (1100) may include an excitation input assembly (1140) that is interconnected with the blade mouting plate (1130) to selectively apply flapwise, edgewise, and/or pitch excitation forces. The restoring spring assemply (1120) may include at least one tuning spring member (1127) positioned adjacent to the primary spring element (1124) used to tune the spring constant or stiffness of the primary spring element (1124) in one of the excitation directions.

Improved accuracy for finite element structural analysis via a new integrated force method

NASA Technical Reports Server (NTRS)

Patnaik, Surya N.; Hopkins, Dale A.; Aiello, Robert A.; Berke, Laszlo

1992-01-01

A comparative study was carried out to determine the accuracy of finite element analyses based on the stiffness method, a mixed method, and the new integrated force and dual integrated force methods. The numerical results were obtained with the following software: MSC/NASTRAN and ASKA for the stiffness method; an MHOST implementation method for the mixed method; and GIFT for the integrated force methods. The results indicate that on an overall basis, the stiffness and mixed methods present some limitations. The stiffness method generally requires a large number of elements in the model to achieve acceptable accuracy. The MHOST method tends to achieve a higher degree of accuracy for course models than does the stiffness method implemented by MSC/NASTRAN and ASKA. The two integrated force methods, which bestow simultaneous emphasis on stress equilibrium and strain compatibility, yield accurate solutions with fewer elements in a model. The full potential of these new integrated force methods remains largely unexploited, and they hold the promise of spawning new finite element structural analysis tools.

METHOD OF MAKING WIRE FUEL ELEMENTS

DOEpatents

Zambrow, J.L.

1960-08-01

A method is given for making a nuclear reactor fuel element in the form of a uranium-bearing wire clad with zirconium. A uranium bar is enclosed in a zirconium sheath which is coated with an oxide of magnesium, beryllium, or zirconium. The sheathed bar is then placed in a steel tube and reduced to the desired diameter by swaging at 800 to 900 deg C, after which the steel and oxide are removed.

Three-dimensional Stress Analysis Using the Boundary Element Method

NASA Technical Reports Server (NTRS)

Wilson, R. B.; Banerjee, P. K.

1984-01-01

The boundary element method is to be extended (as part of the NASA Inelastic Analysis Methods program) to the three-dimensional stress analysis of gas turbine engine hot section components. The analytical basis of the method (as developed in elasticity) is outlined, its numerical implementation is summarized, and the approaches to be followed in extending the method to include inelastic material response indicated.

A method for trace element determination of single Daphnia specimens using total reflection X-ray fluorescence spectrometry

NASA Astrophysics Data System (ADS)

Mages, Margarete; Woelfl, Stefan; v. Tümpling jun, Wolf

2001-11-01

Two new preparation techniques for total-reflection X-ray fluorescence (TXRF) element determination of single freshwater crustacean specimens (dry weight: 3-40 μg ind -1) have been developed and tested using Daphnia pulex from a deep, oligotrophic freshwater lake located in southern Chile. Dry method: Specimens were washed with 0.2 μm filtered lake water and frozen in liquid nitrogen. The freeze-dried Daphnia specimens were weighed using an ultra-fine microbalance and placed on quartz glass carriers for TXRF analysis. Wet method: Specimens were washed with 0.2 μm filtered lake water and placed on quartz glass carriers for TXRF analysis and dried in air. The dry weight was determined using the previously established body length-dry weight relationship. Method validation for both the dry and the wet preparation method in combination with TXRF spectrometry for the element determination in small single freshwater crustaceans showed that both methods can be used for routine investigations. There were no significant differences between the dry and the wet methods concerning the elements Ca, K, Fe, Zn, Br, P, Cu, but the determination of Mn, S and Sr revealed significant differences between the two methods. It seems that the dry method yields more precise results, but the wet method is easier to handle in the field when samples cannot be fixed with liquid nitrogen.

Bending behaviors of fully covered biodegradable polydioxanone biliary stent for human body by finite element method.

PubMed

Liu, Yanhui; Zhu, Guoqing; Yang, Huazhe; Wang, Conger; Zhang, Peihua; Han, Guangting

2018-01-01

This paper presents a study of the bending flexibility of fully covered biodegradable polydioxanone biliary stents (FCBPBs) developed for human body. To investigate the relationship between the bending load and structure parameter (monofilament diameter and braid-pin number), biodegradable polydioxanone biliary stents derived from braiding method were covered with membrane prepared via electrospinning method, and nine FCBPBSs were then obtained for bending test to evaluate the bending flexibility. In addition, by the finite element method, nine numerical models based on actual biliary stent were established and the bending load was calculated through the finite element method. Results demonstrate that the simulation and experimental results are in good agreement with each other, indicating that the simulation results can be provided a useful reference to the investigation of biliary stents. Furthermore, the stress distribution on FCBPBSs was studied, and the plastic dissipation analysis and plastic strain of FCBPBSs were obtained via the bending simulation. Copyright © 2017 Elsevier Ltd. All rights reserved.

Method of holding optical elements without deformation during their fabrication

DOEpatents

Hed, P.P.

1997-04-29

An improved method for securing and removing an optical element to and from a blocking tool without causing deformation of the optical element is disclosed. A lens tissue is placed on the top surface of the blocking tool. Dots of UV cement are applied to the lens tissue without any of the dots contacting each other. An optical element is placed on top of the blocking tool with the lens tissue sandwiched therebetween. The UV cement is then cured. After subsequent fabrication steps, the bonded blocking tool, lens tissue, and optical element are placed in a debonding solution to soften the UV cement. The optical element is then removed from the blocking tool. 16 figs.

Finite-Element Analysis of a Mach-8 Flight Test Article Using Nonlinear Contact Elements

NASA Technical Reports Server (NTRS)

Richards, W. Lance

1997-01-01

A flight test article, called a glove, is required for a Mach-8 boundary-layer experiment to be conducted on a flight mission of the air-launched Pegasus(reg) space booster. The glove is required to provide a smooth, three-dimensional, structurally stable, aerodynamic surface and includes instrumentation to determine when and where boundary-layer transition occurs during the hypersonic flight trajectory. A restraint mechanism has been invented to attach the glove to the wing of the space booster. The restraint mechanism securely attaches the glove to the wing in directions normal to the wing/glove interface surface, but allows the glove to thermally expand and contract to alleviate stresses in directions parallel to the interface surface. A finite-element analysis has been performed using nonlinear contact elements to model the complex behavior of the sliding restraint mechanism. This paper provides an overview of the glove design and presents details of the analysis that were essential to demonstrate the flight worthiness of the wing-glove test article. Results show that all glove components are well within the allowable stress and deformation requirements to satisfy the objectives of the flight research experiment.

Flow tests of a single fuel element coolant channel for a compact fast reactor for space power

NASA Technical Reports Server (NTRS)

Springborn, R. H.

1971-01-01

Water flow tests were conducted on a single-fuel-element cooling channel for a nuclear concept to be used for space power. The tests established a method for measuring coolant flow rate which is applicable to water flow testing of a complete mockup of the reference reactor. The inlet plenum-to-outlet plenum pressure drop, which approximates the overall core pressure drop, was measured and correlated with flow rate. This information can be used for reactor coolant flow and heat transfer calculations. An analytical study of the flow characteristics was also conducted.

Method and apparatus for diagnosing breached fuel elements

DOEpatents

Gross, K.C.; Lambert, J.D.B.; Nomura, S.

1987-03-02

The invention provides an apparatus and method for diagnosing breached fuel elements in a nuclear reactor. A detection system measures the activity of isotopes from the cover gas in the reactor. A data acquisition and processing system monitors the detection system and corrects for the effects of the cover-gas clean up system on the measured activity and further calculates the derivative curve of the corrected activity as a function of time. A plotting system graphs the derivative curve, which represents the instantaneous release rate of fission gas from a breached fuel element. 8 figs.

Creating a Test Validated Structural Dynamic Finite Element Model of the Multi-Utility Technology Test Bed Aircraft

NASA Technical Reports Server (NTRS)

Pak, Chan-Gi; Truong, Samson S.

2014-01-01

Small modeling errors in the finite element model will eventually induce errors in the structural flexibility and mass, thus propagating into unpredictable errors in the unsteady aerodynamics and the control law design. One of the primary objectives of Multi Utility Technology Test Bed, X-56A, aircraft is the flight demonstration of active flutter suppression, and therefore in this study, the identification of the primary and secondary modes for the structural model tuning based on the flutter analysis of X-56A. The ground vibration test validated structural dynamic finite element model of the X-56A is created in this study. The structural dynamic finite element model of the X-56A is improved using a model tuning tool. In this study, two different weight configurations of the X-56A have been improved in a single optimization run.

An adaptive finite element method for the inequality-constrained Reynolds equation

NASA Astrophysics Data System (ADS)

Gustafsson, Tom; Rajagopal, Kumbakonam R.; Stenberg, Rolf; Videman, Juha

2018-07-01

We present a stabilized finite element method for the numerical solution of cavitation in lubrication, modeled as an inequality-constrained Reynolds equation. The cavitation model is written as a variable coefficient saddle-point problem and approximated by a residual-based stabilized method. Based on our recent results on the classical obstacle problem, we present optimal a priori estimates and derive novel a posteriori error estimators. The method is implemented as a Nitsche-type finite element technique and shown in numerical computations to be superior to the usually applied penalty methods.

The NonConforming Virtual Element Method for the Stokes Equations

DOE PAGES

Cangiani, Andrea; Gyrya, Vitaliy; Manzini, Gianmarco

2016-01-01

In this paper, we present the nonconforming virtual element method (VEM) for the numerical approximation of velocity and pressure in the steady Stokes problem. The pressure is approximated using discontinuous piecewise polynomials, while each component of the velocity is approximated using the nonconforming virtual element space. On each mesh element the local virtual space contains the space of polynomials of up to a given degree, plus suitable nonpolynomial functions. The virtual element functions are implicitly defined as the solution of local Poisson problems with polynomial Neumann boundary conditions. As typical in VEM approaches, the explicit evaluation of the non-polynomial functionsmore » is not required. This approach makes it possible to construct nonconforming (virtual) spaces for any polynomial degree regardless of the parity, for two- and three-dimensional problems, and for meshes with very general polygonal and polyhedral elements. We show that the nonconforming VEM is inf-sup stable and establish optimal a priori error estimates for the velocity and pressure approximations. Finally, numerical examples confirm the convergence analysis and the effectiveness of the method in providing high-order accurate approximations.« less

A stochastic method for computing hadronic matrix elements

DOE PAGES

Alexandrou, Constantia; Constantinou, Martha; Dinter, Simon; ...

2014-01-24

In this study, we present a stochastic method for the calculation of baryon 3-point functions which is an alternative to the typically used sequential method offering more versatility. We analyze the scaling of the error of the stochastically evaluated 3-point function with the lattice volume and find a favorable signal to noise ratio suggesting that the stochastic method can be extended to large volumes providing an efficient approach to compute hadronic matrix elements and form factors.

Testing for change in structural elements of forest inventories

Treesearch

Melinda Vokoun; David Wear; Robert Abt

2009-01-01

In this article we develop a methodology to test for changes in the underlying relationships between measures of forest productivity (structural elements) and site characteristics, herein referred to as structural changes, using standard forest inventories. Changes in measures of forest growing stock volume and number of trees for both...

Updating the Finite Element Model of the Aerostructures Test Wing Using Ground Vibration Test Data

NASA Technical Reports Server (NTRS)

Lung, Shun-Fat; Pak, Chan-Gi

2009-01-01

Improved and/or accelerated decision making is a crucial step during flutter certification processes. Unfortunately, most finite element structural dynamics models have uncertainties associated with model validity. Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. The model tuning process requires not only satisfactory correlations between analytical and experimental results, but also the retention of the mass and stiffness properties of the structures. Minimizing the difference between analytical and experimental results is a type of optimization problem. By utilizing the multidisciplinary design, analysis, and optimization (MDAO) tool in order to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes can be matched to the target data to retain the mass matrix orthogonality. This approach has been applied to minimize the model uncertainties for the structural dynamics model of the aerostructures test wing (ATW), which was designed and tested at the National Aeronautics and Space Administration Dryden Flight Research Center (Edwards, California). This study has shown that natural frequencies and corresponding mode shapes from the updated finite element model have excellent agreement with corresponding measured data.

Updating the Finite Element Model of the Aerostructures Test Wing using Ground Vibration Test Data

NASA Technical Reports Server (NTRS)

Lung, Shun-fat; Pak, Chan-gi

2009-01-01

Improved and/or accelerated decision making is a crucial step during flutter certification processes. Unfortunately, most finite element structural dynamics models have uncertainties associated with model validity. Tuning the finite element model using measured data to minimize the model uncertainties is a challenging task in the area of structural dynamics. The model tuning process requires not only satisfactory correlations between analytical and experimental results, but also the retention of the mass and stiffness properties of the structures. Minimizing the difference between analytical and experimental results is a type of optimization problem. By utilizing the multidisciplinary design, analysis, and optimization (MDAO) tool in order to optimize the objective function and constraints; the mass properties, the natural frequencies, and the mode shapes can be matched to the target data to retain the mass matrix orthogonality. This approach has been applied to minimize the model uncertainties for the structural dynamics model of the Aerostructures Test Wing (ATW), which was designed and tested at the National Aeronautics and Space Administration (NASA) Dryden Flight Research Center (DFRC) (Edwards, California). This study has shown that natural frequencies and corresponding mode shapes from the updated finite element model have excellent agreement with corresponding measured data.

Discontinuous Galerkin Finite Element Method for Parabolic Problems

NASA Technical Reports Server (NTRS)

Kaneko, Hideaki; Bey, Kim S.; Hou, Gene J. W.

2004-01-01

In this paper, we develop a time and its corresponding spatial discretization scheme, based upon the assumption of a certain weak singularity of parallel ut(t) parallel Lz(omega) = parallel ut parallel2, for the discontinuous Galerkin finite element method for one-dimensional parabolic problems. Optimal convergence rates in both time and spatial variables are obtained. A discussion of automatic time-step control method is also included.

Finite element method for optimal guidance of an advanced launch vehicle

NASA Technical Reports Server (NTRS)

Hodges, Dewey H.; Bless, Robert R.; Calise, Anthony J.; Leung, Martin

1992-01-01

A temporal finite element based on a mixed form of Hamilton's weak principle is summarized for optimal control problems. The resulting weak Hamiltonian finite element method is extended to allow for discontinuities in the states and/or discontinuities in the system equations. An extension of the formulation to allow for control inequality constraints is also presented. The formulation does not require element quadrature, and it produces a sparse system of nonlinear algebraic equations. To evaluate its feasibility for real-time guidance applications, this approach is applied to the trajectory optimization of a four-state, two-stage model with inequality constraints for an advanced launch vehicle. Numerical results for this model are presented and compared to results from a multiple-shooting code. The results show the accuracy and computational efficiency of the finite element method.

A Runge-Kutta discontinuous finite element method for high speed flows

NASA Technical Reports Server (NTRS)

Bey, Kim S.; Oden, J. T.

1991-01-01

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

A finite element conjugate gradient FFT method for scattering

NASA Technical Reports Server (NTRS)

Collins, Jeffery D.; Ross, Dan; Jin, J.-M.; Chatterjee, A.; Volakis, John L.

1991-01-01

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

A proposed method for enhanced eigen-pair extraction using finite element methods: Theory and application

NASA Technical Reports Server (NTRS)

Jara-Almonte, J.; Mitchell, L. D.

1988-01-01

The paper covers two distinct parts: theory and application. The goal of this work was the reduction of model size with an increase in eigenvalue/vector accuracy. This method is ideal for the condensation of large truss- or beam-type structures. The theoretical approach involves the conversion of a continuum transfer matrix beam element into an 'Exact' dynamic stiffness element. This formulation is implemented in a finite element environment. This results in the need to solve a transcendental eigenvalue problem. Once the eigenvalue is determined the eigenvectors can be reconstructed with any desired spatial precision. No discretization limitations are imposed on the reconstruction. The results of such a combined finite element and transfer matrix formulation is a much smaller FEM eigenvalue problem. This formulation has the ability to extract higher eigenvalues as easily and as accurately as lower eigenvalues. Moreover, one can extract many more eigenvalues/vectors from the model than the number of degrees of freedom in the FEM formulation. Typically, the number of eigenvalues accurately extractable via the 'Exact' element method are at least 8 times the number of degrees of freedom. In contrast, the FEM usually extracts one accurate (within 5 percent) eigenvalue for each 3-4 degrees of freedom. The 'Exact' element results in a 20-30 improvement in the number of accurately extractable eigenvalues and eigenvectors.

Validity of Combining History Elements and Physical Examination Tests to Diagnose Patellofemoral Pain.

PubMed

Décary, Simon; Frémont, Pierre; Pelletier, Bruno; Fallaha, Michel; Belzile, Sylvain; Martel-Pelletier, Johanne; Pelletier, Jean-Pierre; Feldman, Debbie; Sylvestre, Marie-Pierre; Vendittoli, Pascal-André; Desmeules, François

2018-04-01

To assess the validity of diagnostic clusters combining history elements and physical examination tests to diagnose or exclude patellofemoral pain (PFP). Prospective diagnostic study. Orthopedic outpatient clinics, family medicine clinics, and community-dwelling. Consecutive patients (N=279) consulting one of the participating orthopedic surgeons (n=3) or sport medicine physicians (n=2) for any knee complaint. Not applicable. History elements and physical examination tests were obtained by a trained physiotherapist blinded to the reference standard: a composite diagnosis including both physical examination tests and imaging results interpretation performed by an expert physician. Penalized logistic regression (least absolute shrinkage and selection operator) was used to identify history elements and physical examination tests associated with the diagnosis of PFP, and recursive partitioning was used to develop diagnostic clusters. Diagnostic accuracy measures including sensitivity, specificity, positive and negative predictive values, and positive and negative likelihood ratios with associated 95% confidence intervals (CIs) were calculated. Two hundred seventy-nine participants were evaluated, and 75 had a diagnosis of PFP (26.9%). Different combinations of history elements and physical examination tests including the age of participants, knee pain location, difficulty descending stairs, patellar facet palpation, and passive knee extension range of motion were associated with a diagnosis of PFP and used in clusters to accurately discriminate between individuals with PFP and individuals without PFP. Two diagnostic clusters developed to confirm the presence of PFP yielded a positive likelihood ratio of 8.7 (95% CI, 5.2-14.6) and 3 clusters to exclude PFP yielded a negative likelihood ratio of .12 (95% CI, .06-.27). Diagnostic clusters combining common history elements and physical examination tests that can accurately diagnose or exclude PFP compared to various knee

A partially penalty immersed Crouzeix-Raviart finite element method for interface problems.

PubMed

An, Na; Yu, Xijun; Chen, Huanzhen; Huang, Chaobao; Liu, Zhongyan

2017-01-01

The elliptic equations with discontinuous coefficients are often used to describe the problems of the multiple materials or fluids with different densities or conductivities or diffusivities. In this paper we develop a partially penalty immersed finite element (PIFE) method on triangular grids for anisotropic flow models, in which the diffusion coefficient is a piecewise definite-positive matrix. The standard linear Crouzeix-Raviart type finite element space is used on non-interface elements and the piecewise linear Crouzeix-Raviart type immersed finite element (IFE) space is constructed on interface elements. The piecewise linear functions satisfying the interface jump conditions are uniquely determined by the integral averages on the edges as degrees of freedom. The PIFE scheme is given based on the symmetric, nonsymmetric or incomplete interior penalty discontinuous Galerkin formulation. The solvability of the method is proved and the optimal error estimates in the energy norm are obtained. Numerical experiments are presented to confirm our theoretical analysis and show that the newly developed PIFE method has optimal-order convergence in the [Formula: see text] norm as well. In addition, numerical examples also indicate that this method is valid for both the isotropic and the anisotropic elliptic interface problems.

METHOD OF MAKING FUEL ELEMENTS

DOEpatents

Bean, C.H.; Macherey, R.E.

1959-12-01

A method is described for fabricating fuel elements, particularly for enclosing a plate of metal with a second metal by inserting the plate into an aperture of a frame of a second plate, placing a sheet of the second metal on each of opposite faces of the assembled plate and frame, purging with an inert gas the air from the space within the frame and the sheets while sealing the seams between the frame and the sheets, exhausting the space, purging the space with air, re-exhausting the spaces, sealing the second aperture, and applying heat and pressure to bond the sheets, the plate, and the frame to one another.

Finite cover method with mortar elements for elastoplasticity problems

NASA Astrophysics Data System (ADS)

Kurumatani, M.; Terada, K.

2005-06-01

Finite cover method (FCM) is extended to elastoplasticity problems. The FCM, which was originally developed under the name of manifold method, has recently been recognized as one of the generalized versions of finite element methods (FEM). Since the mesh for the FCM can be regular and squared regardless of the geometry of structures to be analyzed, structural analysts are released from a burdensome task of generating meshes conforming to physical boundaries. Numerical experiments are carried out to assess the performance of the FCM with such discretization in elastoplasticity problems. Particularly to achieve this accurately, the so-called mortar elements are introduced to impose displacement boundary conditions on the essential boundaries, and displacement compatibility conditions on material interfaces of two-phase materials or on joint surfaces between mutually incompatible meshes. The validity of the mortar approximation is also demonstrated in the elastic-plastic FCM.

Test methods for determining the suitability of metal alloys for use in oxygen-enriched environments

NASA Technical Reports Server (NTRS)

Stoltzfus, Joel M.; Gunaji, Mohan V.

1991-01-01

Materials are more flammable in oxygen rich environments than in air. When the structural elements of a system containing oxygen ignite and burn, the results are often catastrophic, causing loss of equipment and perhaps even human lives. Therefore, selection of the proper metallic and non-metallic materials for use in oxygen systems is extremely important. While test methods for the selection of non-metallic materials have been available for years, test methods for the selection of alloys have not been available until recently. Presented here are several test methods that were developed recently at NASA's White Sands Test Facility (WSTF) to study the ignition and combustion of alloys, including the supersonic and subsonic speed particle impact tests, the frictional heating and coefficient of friction tests, and the promoted combustion test. These test methods are available for commercial use.

Nitsche Extended Finite Element Methods for Earthquake Simulation

NASA Astrophysics Data System (ADS)

Coon, Ethan T.

Modeling earthquakes and geologically short-time-scale events on fault networks is a difficult problem with important implications for human safety and design. These problems demonstrate a. rich physical behavior, in which distributed loading localizes both spatially and temporally into earthquakes on fault systems. This localization is governed by two aspects: friction and fault geometry. Computationally, these problems provide a stern challenge for modelers --- static and dynamic equations must be solved on domains with discontinuities on complex fault systems, and frictional boundary conditions must be applied on these discontinuities. The most difficult aspect of modeling physics on complicated domains is the mesh. Most numerical methods involve meshing the geometry; nodes are placed on the discontinuities, and edges are chosen to coincide with faults. The resulting mesh is highly unstructured, making the derivation of finite difference discretizations difficult. Therefore, most models use the finite element method. Standard finite element methods place requirements on the mesh for the sake of stability, accuracy, and efficiency. The formation of a mesh which both conforms to fault geometry and satisfies these requirements is an open problem, especially for three dimensional, physically realistic fault. geometries. In addition, if the fault system evolves over the course of a dynamic simulation (i.e. in the case of growing cracks or breaking new faults), the geometry must he re-meshed at each time step. This can be expensive computationally. The fault-conforming approach is undesirable when complicated meshes are required, and impossible to implement when the geometry is evolving. Therefore, meshless and hybrid finite element methods that handle discontinuities without placing them on element boundaries are a desirable and natural way to discretize these problems. Several such methods are being actively developed for use in engineering mechanics involving crack

[Research Progress and Prospect of Applications of Finite Element Method in Lumbar Spine Biomechanics].

PubMed

Zhang, Zhenjun; Li, Yang; Liao, Zhenhua; Liu, Weiqiang

2016-12-01

Based on the application of finite element analysis in spine biomechanics,the research progress of finite element method applied in lumbar spine mechanics is reviewed and the prospect is forecasted.The related works,including lumbar ontology modeling,clinical application research,and occupational injury and protection,are summarized.The main research areas of finite element method are as follows:new accurate modeling process,the optimized simulation method,diversified clinical effect evaluation,and the clinical application of artificial lumbar disc.According to the recent research progress,the application prospects of finite element method,such as automation and individuation of modeling process,evaluation and analysis of new operation methods and simulation of mechanical damage and dynamic response,are discussed.The purpose of this paper is to provide the theoretical reference and practical guidance for the clinical lumbar problems by reviewing the application of finite element method in the field of the lumbar spine biomechanics.

Robust Hybrid Finite Element Methods for Antennas and Microwave Circuits

NASA Technical Reports Server (NTRS)

Gong, J.; Volakis, John L.

1996-01-01

One of the primary goals in this dissertation is concerned with the development of robust hybrid finite element-boundary integral (FE-BI) techniques for modeling and design of conformal antennas of arbitrary shape. Both the finite element and integral equation methods will be first overviewed in this chapter with an emphasis on recently developed hybrid FE-BI methodologies for antennas, microwave and millimeter wave applications. The structure of the dissertation is then outlined. We conclude the chapter with discussions of certain fundamental concepts and methods in electromagnetics, which are important to this study.

Injector element characterization methodology

NASA Technical Reports Server (NTRS)

Cox, George B., Jr.

1988-01-01

Characterization of liquid rocket engine injector elements is an important part of the development process for rocket engine combustion devices. Modern nonintrusive instrumentation for flow velocity and spray droplet size measurement, and automated, computer-controlled test facilities allow rapid, low-cost evaluation of injector element performance and behavior. Application of these methods in rocket engine development, paralleling their use in gas turbine engine development, will reduce rocket engine development cost and risk. The Alternate Turbopump (ATP) Hot Gas Systems (HGS) preburner injector elements were characterized using such methods, and the methodology and some of the results obtained will be shown.

Error analysis of finite element method for Poisson–Nernst–Planck equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sun, Yuzhou; Sun, Pengtao; Zheng, Bin

A priori error estimates of finite element method for time-dependent Poisson-Nernst-Planck equations are studied in this work. We obtain the optimal error estimates in L∞(H1) and L2(H1) norms, and suboptimal error estimates in L∞(L2) norm, with linear element, and optimal error estimates in L∞(L2) norm with quadratic or higher-order element, for both semi- and fully discrete finite element approximations. Numerical experiments are also given to validate the theoretical results.

Three-dimensional photoacoustic tomography based on graphics-processing-unit-accelerated finite element method.

PubMed

Peng, Kuan; He, Ling; Zhu, Ziqiang; Tang, Jingtian; Xiao, Jiaying

2013-12-01

Compared with commonly used analytical reconstruction methods, the frequency-domain finite element method (FEM) based approach has proven to be an accurate and flexible algorithm for photoacoustic tomography. However, the FEM-based algorithm is computationally demanding, especially for three-dimensional cases. To enhance the algorithm's efficiency, in this work a parallel computational strategy is implemented in the framework of the FEM-based reconstruction algorithm using a graphic-processing-unit parallel frame named the "compute unified device architecture." A series of simulation experiments is carried out to test the accuracy and accelerating effect of the improved method. The results obtained indicate that the parallel calculation does not change the accuracy of the reconstruction algorithm, while its computational cost is significantly reduced by a factor of 38.9 with a GTX 580 graphics card using the improved method.

A Priori Method of Using Photon Activation Analysis to Determine Unknown Trace Element Concentrations in NIST Standards

DOE Office of Scientific and Technical Information (OSTI.GOV)

Green, Jaromy; Sun Zaijing; Wells, Doug

2009-03-10

Photon activation analysis detected elements in two NIST standards that did not have reported concentration values. A method is currently being developed to infer these concentrations by using scaling parameters and the appropriate known quantities within the NIST standard itself. Scaling parameters include: threshold, peak and endpoint energies; photo-nuclear cross sections for specific isotopes; Bremstrahlung spectrum; target thickness; and photon flux. Photo-nuclear cross sections and energies from the unknown elements must also be known. With these quantities, the same integral was performed for both the known and unknown elements resulting in an inference of the concentration of the un-reported elementmore » based on the reported value. Since Rb and Mn were elements that were reported in the standards, and because they had well-identified peaks, they were used as the standards of inference to determine concentrations of the unreported elements of As, I, Nb, Y, and Zr. This method was tested by choosing other known elements within the standards and inferring a value based on the stated procedure. The reported value of Mn in the first NIST standard was 403{+-}15 ppm and the reported value of Ca in the second NIST standard was 87000 ppm (no reported uncertainty). The inferred concentrations were 370{+-}23 ppm and 80200{+-}8700 ppm respectively.« less

High-Resolution Genuinely Multidimensional Solution of Conservation Laws by the Space-Time Conservation Element and Solution Element Method

NASA Technical Reports Server (NTRS)

Himansu, Ananda; Chang, Sin-Chung; Yu, Sheng-Tao; Wang, Xiao-Yen; Loh, Ching-Yuen; Jorgenson, Philip C. E.

1999-01-01

In this overview paper, we review the basic principles of the method of space-time conservation element and solution element for solving the conservation laws in one and two spatial dimensions. The present method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. In contrast to the modern upwind schemes, the approach here does not use the Riemann solver and the reconstruction procedure as the building blocks. The drawbacks of the upwind approach, such as the difficulty of rationally extending the 1D scalar approach to systems of equations and particularly to multiple dimensions is here contrasted with the uniformity and ease of generalization of the Conservation Element and Solution Element (CE/SE) 1D scalar schemes to systems of equations and to multiple spatial dimensions. The assured compatibility with the simplest type of unstructured meshes, and the uniquely simple nonreflecting boundary conditions of the present method are also discussed. The present approach has yielded high-resolution shocks, rarefaction waves, acoustic waves, vortices, ZND detonation waves, and shock/acoustic waves/vortices interactions. Moreover, since no directional splitting is employed, numerical resolution of two-dimensional calculations is comparable to that of the one-dimensional calculations. Some sample applications displaying the strengths and broad applicability of the CE/SE method are reviewed.

Adaptive finite element methods for two-dimensional problems in computational fracture mechanics

NASA Technical Reports Server (NTRS)

Min, J. B.; Bass, J. M.; Spradley, L. W.

1994-01-01

Some recent results obtained using solution-adaptive finite element methods in two-dimensional problems in linear elastic fracture mechanics are presented. The focus is on the basic issue of adaptive finite element methods for validating the new methodology by computing demonstration problems and comparing the stress intensity factors to analytical results.

A hybridized formulation for the weak Galerkin mixed finite element method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Wang, Junping; Ye, Xiu

This paper presents a hybridized formulation for the weak Galerkin mixed finite element method (WG-MFEM) which was introduced and analyzed in Wang and Ye (2014) for second order elliptic equations. The WG-MFEM method was designed by using discontinuous piecewise polynomials on finite element partitions consisting of polygonal or polyhedral elements of arbitrary shape. The key to WG-MFEM is the use of a discrete weak divergence operator which is defined and computed by solving inexpensive problems locally on each element. The hybridized formulation of this paper leads to a significantly reduced system of linear equations involving only the unknowns arising frommore » the Lagrange multiplier in hybridization. Optimal-order error estimates are derived for the hybridized WG-MFEM approximations. In conclusion, some numerical results are reported to confirm the theory and a superconvergence for the Lagrange multiplier.« less

A hybridized formulation for the weak Galerkin mixed finite element method

DOE PAGES

Mu, Lin; Wang, Junping; Ye, Xiu

2016-01-14

This paper presents a hybridized formulation for the weak Galerkin mixed finite element method (WG-MFEM) which was introduced and analyzed in Wang and Ye (2014) for second order elliptic equations. The WG-MFEM method was designed by using discontinuous piecewise polynomials on finite element partitions consisting of polygonal or polyhedral elements of arbitrary shape. The key to WG-MFEM is the use of a discrete weak divergence operator which is defined and computed by solving inexpensive problems locally on each element. The hybridized formulation of this paper leads to a significantly reduced system of linear equations involving only the unknowns arising frommore » the Lagrange multiplier in hybridization. Optimal-order error estimates are derived for the hybridized WG-MFEM approximations. In conclusion, some numerical results are reported to confirm the theory and a superconvergence for the Lagrange multiplier.« less

Rolling-element fatigue life of silicon nitride balls: Preliminary test results

NASA Technical Reports Server (NTRS)

Parker, R. J.; Zaretsky, E. V.

1972-01-01

Hot pressed silicon nitride was evaluated as a rolling element bearing material. The five-ball fatigue tester was used to test 12.7 mm (0.500 in.) diameter balls at a maximum Hertz stress of 800,000 psi at a race temperature of 130 F. The fatigue spalls in the silicon nitride resembled those in typical bearing steels. The ten-percent fatigue life of the silicon nitride balls was approximately one-eighth to one-fifth that of typical bearing steels (52100 and M-50). The load capacity of the silicon nitride was approximately one-third that of typical bearing steels. The load capacity of the silicon nitride was significantly higher than previously tested ceramic materials for rolling element bearings.

An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics

NASA Astrophysics Data System (ADS)

Hahn, Oliver; Angulo, Raul E.

2016-01-01

N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new approach to simulate the gravitational evolution of cold collisionless fluids by solving the Vlasov-Poisson equations in terms of adaptively refineable `Lagrangian phase-space elements'. These geometrical elements are piecewise smooth maps between Lagrangian space and Eulerian phase-space and approximate the continuum structure of the distribution function. They allow for dynamical adaptive splitting to accurately follow the evolution even in regions of very strong mixing. We discuss in detail various one-, two- and three-dimensional test problems to demonstrate the performance of our method. Its advantages compared to N-body algorithms are: (I) explicit tracking of the fine-grained distribution function, (II) natural representation of caustics, (III) intrinsically smooth gravitational potential fields, thus (IV) eliminating the need for any type of ad hoc force softening. We show the potential of our method by simulating structure formation in a warm dark matter scenario. We discuss how spurious collisionality and large-scale discreteness noise of N-body methods are both strongly suppressed, which eliminates the artificial fragmentation of filaments. Therefore, we argue that our new approach improves on the N-body method when simulating self-gravitating cold and collisionless fluids, and is the first method that allows us to explicitly follow the fine-grained evolution in six-dimensional phase-space.

Satellite Vibration Testing: Angle optimisation method to Reduce Overtesting

NASA Astrophysics Data System (ADS)

Knight, Charly; Remedia, Marcello; Aglietti, Guglielmo S.; Richardson, Guy

2018-06-01

Spacecraft overtesting is a long running problem, and the main focus of most attempts to reduce it has been to adjust the base vibration input (i.e. notching). Instead this paper examines testing alternatives for secondary structures (equipment) coupled to the main structure (satellite) when they are tested separately. Even if the vibration source is applied along one of the orthogonal axes at the base of the coupled system (satellite plus equipment), the dynamics of the system and potentially the interface configuration mean the vibration at the interface may not occur all along one axis much less the corresponding orthogonal axis of the base excitation. This paper proposes an alternative testing methodology in which the testing of a piece of equipment occurs at an offset angle. This Angle Optimisation method may have multiple tests but each with an altered input direction allowing for the best match between all specified equipment system responses with coupled system tests. An optimisation process that compares the calculated equipment RMS values for a range of inputs with the maximum coupled system RMS values, and is used to find the optimal testing configuration for the given parameters. A case study was performed to find the best testing angles to match the acceleration responses of the centre of mass and sum of interface forces for all three axes, as well as the von Mises stress for an element by a fastening point. The angle optimisation method resulted in RMS values and PSD responses that were much closer to the coupled system when compared with traditional testing. The optimum testing configuration resulted in an overall average error significantly smaller than the traditional method. Crucially, this case study shows that the optimum test campaign could be a single equipment level test opposed to the traditional three orthogonal direction tests.

The Method of Space-time Conservation Element and Solution Element: Development of a New Implicit Solver

NASA Technical Reports Server (NTRS)

Chang, S. C.; Wang, X. Y.; Chow, C. Y.; Himansu, A.

1995-01-01

The method of space-time conservation element and solution element is a nontraditional numerical method designed from a physicist's perspective, i.e., its development is based more on physics than numerics. It uses only the simplest approximation techniques and yet is capable of generating nearly perfect solutions for a 2-D shock reflection problem used by Helen Yee and others. In addition to providing an overall view of the new method, we introduce a new concept in the design of implicit schemes, and use it to construct a highly accurate solver for a convection-diffusion equation. It is shown that, in the inviscid case, this new scheme becomes explicit and its amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, its principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme.

Finite Element Analysis and Test Results Comparison for the Hybrid Wing Body Center Section Test Article

NASA Technical Reports Server (NTRS)

Przekop, Adam; Jegley, Dawn C.; Rouse, Marshall; Lovejoy, Andrew E.

2016-01-01

This report documents the comparison of test measurements and predictive finite element analysis results for a hybrid wing body center section test article. The testing and analysis efforts were part of the Airframe Technology subproject within the NASA Environmentally Responsible Aviation project. Test results include full field displacement measurements obtained from digital image correlation systems and discrete strain measurements obtained using both unidirectional and rosette resistive gauges. Most significant results are presented for the critical five load cases exercised during the test. Final test to failure after inflicting severe damage to the test article is also documented. Overall, good comparison between predicted and actual behavior of the test article is found.

A new weak Galerkin finite element method for elliptic interface problems

DOE PAGES

Mu, Lin; Wang, Junping; Ye, Xiu; ...

2016-08-26

We introduce and analyze a new weak Galerkin (WG) finite element method in this paper for solving second order elliptic equations with discontinuous coefficients and interfaces. Comparing with the existing WG algorithm for solving the same type problems, the present WG method has a simpler variational formulation and fewer unknowns. Moreover, the new WG algorithm allows the use of finite element partitions consisting of general polytopal meshes and can be easily generalized to high orders. Optimal order error estimates in both H1 and L2 norms are established for the present WG finite element solutions. We conducted extensive numerical experiments inmore » order to examine the accuracy, flexibility, and robustness of the proposed WG interface approach. In solving regular elliptic interface problems, high order convergences are numerically confirmed by using piecewise polynomial basis functions of high degrees. Moreover, the WG method is shown to be able to accommodate very complicated interfaces, due to its flexibility in choosing finite element partitions. Finally, in dealing with challenging problems with low regularities, the piecewise linear WG method is capable of delivering a second order of accuracy in L∞ norm for both C1 and H2 continuous solutions.« less

A new weak Galerkin finite element method for elliptic interface problems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Wang, Junping; Ye, Xiu

We introduce and analyze a new weak Galerkin (WG) finite element method in this paper for solving second order elliptic equations with discontinuous coefficients and interfaces. Comparing with the existing WG algorithm for solving the same type problems, the present WG method has a simpler variational formulation and fewer unknowns. Moreover, the new WG algorithm allows the use of finite element partitions consisting of general polytopal meshes and can be easily generalized to high orders. Optimal order error estimates in both H1 and L2 norms are established for the present WG finite element solutions. We conducted extensive numerical experiments inmore » order to examine the accuracy, flexibility, and robustness of the proposed WG interface approach. In solving regular elliptic interface problems, high order convergences are numerically confirmed by using piecewise polynomial basis functions of high degrees. Moreover, the WG method is shown to be able to accommodate very complicated interfaces, due to its flexibility in choosing finite element partitions. Finally, in dealing with challenging problems with low regularities, the piecewise linear WG method is capable of delivering a second order of accuracy in L∞ norm for both C1 and H2 continuous solutions.« less

Updating finite element dynamic models using an element-by-element sensitivity methodology

NASA Technical Reports Server (NTRS)

Farhat, Charbel; Hemez, Francois M.

1993-01-01

A sensitivity-based methodology for improving the finite element model of a given structure using test modal data and a few sensors is presented. The proposed method searches for both the location and sources of the mass and stiffness errors and does not interfere with the theory behind the finite element model while correcting these errors. The updating algorithm is derived from the unconstrained minimization of the squared L sub 2 norms of the modal dynamic residuals via an iterative two-step staggered procedure. At each iteration, the measured mode shapes are first expanded assuming that the model is error free, then the model parameters are corrected assuming that the expanded mode shapes are exact. The numerical algorithm is implemented in an element-by-element fashion and is capable of 'zooming' on the detected error locations. Several simulation examples which demonstate the potential of the proposed methodology are discussed.

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

NASA Technical Reports Server (NTRS)

Feng, Huiyu; Mavriplis, Catherine; VanderWijngaart, Rob; Biswas, Rupak

2004-01-01

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel

A weak Galerkin least-squares finite element method for div-curl systems

NASA Astrophysics Data System (ADS)

Li, Jichun; Ye, Xiu; Zhang, Shangyou

2018-06-01

In this paper, we introduce a weak Galerkin least-squares method for solving div-curl problem. This finite element method leads to a symmetric positive definite system and has the flexibility to work with general meshes such as hybrid mesh, polytopal mesh and mesh with hanging nodes. Error estimates of the finite element solution are derived. The numerical examples demonstrate the robustness and flexibility of the proposed method.

Development of quadrilateral spline thin plate elements using the B-net method

NASA Astrophysics Data System (ADS)

Chen, Juan; Li, Chong-Jun

2013-08-01

The quadrilateral discrete Kirchhoff thin plate bending element DKQ is based on the isoparametric element Q8, however, the accuracy of the isoparametric quadrilateral elements will drop significantly due to mesh distortions. In a previouswork, we constructed an 8-node quadrilateral spline element L8 using the triangular area coordinates and the B-net method, which can be insensitive to mesh distortions and possess the second order completeness in the Cartesian coordinates. In this paper, a thin plate spline element is developed based on the spline element L8 and the refined technique. Numerical examples show that the present element indeed possesses higher accuracy than the DKQ element for distorted meshes.

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

NASA Technical Reports Server (NTRS)

Kaljevic, Igor; Patnaik, Surya N.; Hopkins, Dale A.

1996-01-01

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

An Automated Method for Landmark Identification and Finite-Element Modeling of the Lumbar Spine.

PubMed

Campbell, Julius Quinn; Petrella, Anthony J

2015-11-01

The purpose of this study was to develop a method for the automated creation of finite-element models of the lumbar spine. Custom scripts were written to extract bone landmarks of lumbar vertebrae and assemble L1-L5 finite-element models. End-plate borders, ligament attachment points, and facet surfaces were identified. Landmarks were identified to maintain mesh correspondence between meshes for later use in statistical shape modeling. 90 lumbar vertebrae were processed creating 18 subject-specific finite-element models. Finite-element model surfaces and ligament attachment points were reproduced within 1e-5 mm of the bone surface, including the critical contact surfaces of the facets. Element quality exceeded specifications in 97% of elements for the 18 models created. The current method is capable of producing subject-specific finite-element models of the lumbar spine with good accuracy, quality, and robustness. The automated methods developed represent advancement in the state of the art of subject-specific lumbar spine modeling to a scale not possible with prior manual and semiautomated methods.

Finite element analysis of CFRP reinforced silo structure design method

NASA Astrophysics Data System (ADS)

Yuan, Long; Xu, Xinsheng

2017-11-01

Because of poor construction, there is a serious problem of concrete quality in the silo project, which seriously affects the safe use of the structure. Concrete quality problems are mainly seen in three aspects: concrete strength cannot meet the design requirements, concrete cracking phenomenon is serious, and the unreasonable concrete vibration leads to a lot of honeycombs and surface voids. Silos are usually reinforced by carbon fiber cloth in order to ensure the safe use of silos. By the example of an alumina silo in a fly ash plant in Binzhou, Shandong Province, the alumina silo project was tested and examined on site. According to filed test results, the actual concrete strength was determined, and the damage causes of the silo was analysed. Then, a finite element analysis model of this silo was established, the CFRP cloth reinforcement method was adopted to strengthen the silo, and other technology like additional reinforcement, rebar planting, carbon fiber bonding technology was also expounded. The research of this paper is of great significance to the design and construction of silo structure.

Using Finite Element Method to Estimate the Material Properties of a Bearing Cage

DTIC Science & Technology

2018-02-01

UNCLASSIFIED UNCLASSIFIED AD-E403 988 Technical Report ARMET-TR-17035 USING FINITE ELEMENT METHOD TO ESTIMATE THE MATERIAL...TITLE AND SUBTITLE USING FINITE ELEMENT METHOD TO ESTIMATE THE MATERIAL PROPERTIES OF A BEARING CAGE 5a. CONTRACT NUMBER 5b. GRANT...specifications of non-metallic bearing cages are typically not supplied by the manufacturer. In order to setup a finite element analysis of a

77 FR 16868 - Quality Verification for Plate-Type Uranium-Aluminum Fuel Elements for Use in Research and Test...

Federal Register 2010, 2011, 2012, 2013, 2014

2012-03-22

... Fuel Elements for Use in Research and Test Reactors AGENCY: Nuclear Regulatory Commission. ACTION... Plate-Type Uranium-Aluminum Fuel Elements for Use in Research and Test Reactors.'' This guide describes... plate-type uranium-aluminum fuel elements used in research and test reactors (RTRs). DATES: Submit...

Numerical simulation of self-sustained oscillation of a voice-producing element based on Navier-Stokes equations and the finite element method.

PubMed

de Vries, Martinus P; Hamburg, Marc C; Schutte, Harm K; Verkerke, Gijsbertus J; Veldman, Arthur E P

2003-04-01

Surgical removal of the larynx results in radically reduced production of voice and speech. To improve voice quality a voice-producing element (VPE) is developed, based on the lip principle, called after the lips of a musician while playing a brass instrument. To optimize the VPE, a numerical model is developed. In this model, the finite element method is used to describe the mechanical behavior of the VPE. The flow is described by two-dimensional incompressible Navier-Stokes equations. The interaction between VPE and airflow is modeled by placing the grid of the VPE model in the grid of the aerodynamical model, and requiring continuity of forces and velocities. By applying and increasing pressure to the numerical model, pulses comparable to glottal volume velocity waveforms are obtained. By variation of geometric parameters their influence can be determined. To validate this numerical model, an in vitro test with a prototype of the VPE is performed. Experimental and numerical results show an acceptable agreement.

A new compound control method for sine-on-random mixed vibration test

NASA Astrophysics Data System (ADS)

Zhang, Buyun; Wang, Ruochen; Zeng, Falin

2017-09-01

Vibration environmental test (VET) is one of the important and effective methods to provide supports for the strength design, reliability and durability test of mechanical products. A new separation control strategy was proposed to apply in multiple-input multiple-output (MIMO) sine on random (SOR) mixed mode vibration test, which is the advanced and intensive test type of VET. As the key problem of the strategy, correlation integral method was applied to separate the mixed signals which included random and sinusoidal components. The feedback control formula of MIMO linear random vibration system was systematically deduced in frequency domain, and Jacobi control algorithm was proposed in view of the elements, such as self-spectrum, coherence, and phase of power spectral density (PSD) matrix. Based on the excessive correction of excitation in sine vibration test, compression factor was introduced to reduce the excitation correction, avoiding the destruction to vibration table or other devices. The two methods were synthesized to be applied in MIMO SOR vibration test system. In the final, verification test system with the vibration of a cantilever beam as the control object was established to verify the reliability and effectiveness of the methods proposed in the paper. The test results show that the exceeding values can be controlled in the tolerance range of references accurately, and the method can supply theory and application supports for mechanical engineering.

An Error Analysis for the Finite Element Method Applied to Convection Diffusion Problems.

DTIC Science & Technology

1981-03-01

D TFhG-]NOLOGY k 4b 00 \\" ) ’b Technical Note BN-962 AN ERROR ANALYSIS FOR THE FINITE ELEMENT METHOD APPLIED TO CONVECTION DIFFUSION PROBLEM by I...Babu~ka and W. G. Szym’czak March 1981 V.. UNVI I Of- ’i -S AN ERROR ANALYSIS FOR THE FINITE ELEMENT METHOD P. - 0 w APPLIED TO CONVECTION DIFFUSION ...AOAO98 895 MARYLAND UNIVYCOLLEGE PARK INST FOR PHYSICAL SCIENCE--ETC F/G 12/I AN ERROR ANALYIS FOR THE FINITE ELEMENT METHOD APPLIED TO CONV..ETC (U

A variational numerical method based on finite elements for the nonlinear solution characteristics of the periodically forced Chen system

NASA Astrophysics Data System (ADS)

Khan, Sabeel M.; Sunny, D. A.; Aqeel, M.

2017-09-01

Nonlinear dynamical systems and their solutions are very sensitive to initial conditions and therefore need to be approximated carefully. In this article, we present and analyze nonlinear solution characteristics of the periodically forced Chen system with the application of a variational method based on the concept of finite time-elements. Our approach is based on the discretization of physical time space into finite elements where each time-element is mapped to a natural time space. The solution of the system is then determined in natural time space using a set of suitable basis functions. The numerical algorithm is presented and implemented to compute and analyze nonlinear behavior at different time-step sizes. The obtained results show an excellent agreement with the classical RK-4 and RK-5 methods. The accuracy and convergence of the method is shown by comparing numerically computed results with the exact solution for a test problem. The presented method has shown a great potential in dealing with the solutions of nonlinear dynamical systems and thus can be utilized in delineating different features and characteristics of their solutions.

Application of the Boundary Element Method to Fatigue Crack Growth Analysis

DTIC Science & Technology

1988-09-01

III, and Noetic PROBE in Section IV. Correlation of the boundary element method and modeling techniques employed in this study were shown with the...distribution unlimited I I I Preface! 3 The purpose of this study was to apply the boundary element method (BEM) to two dimensional fracture mechanics...problems, and to use the BEM to analyze the interference effects of holes on cracks through a parametric study of a two hole 3 tension strip. The study

Data processing device test apparatus and method therefor

DOEpatents

Wilcox, Richard Jacob; Mulig, Jason D.; Eppes, David; Bruce, Michael R.; Bruce, Victoria J.; Ring, Rosalinda M.; Cole, Jr., Edward I.; Tangyunyong, Paiboon; Hawkins, Charles F.; Louie, Arnold Y.

2003-04-08

A method and apparatus mechanism for testing data processing devices are implemented. The test mechanism isolates critical paths by correlating a scanning microscope image with a selected speed path failure. A trigger signal having a preselected value is generated at the start of each pattern vector. The sweep of the scanning microscope is controlled by a computer, which also receives and processes the image signals returned from the microscope. The value of the trigger signal is correlated with a set of pattern lines being driven on the DUT. The trigger is either asserted or negated depending the detection of a pattern line failure and the particular line that failed. In response to the detection of the particular speed path failure being characterized, and the trigger signal, the control computer overlays a mask on the image of the device under test (DUT). The overlaid image provides a visual correlation of the failure with the structural elements of the DUT at the level of resolution of the microscope itself.

Flow Applications of the Least Squares Finite Element Method

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan

1998-01-01

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

The constraint method: A new finite element technique. [applied to static and dynamic loads on plates

NASA Technical Reports Server (NTRS)

Tsai, C.; Szabo, B. A.

1973-01-01

An approch to the finite element method which utilizes families of conforming finite elements based on complete polynomials is presented. Finite element approximations based on this method converge with respect to progressively reduced element sizes as well as with respect to progressively increasing orders of approximation. Numerical results of static and dynamic applications of plates are presented to demonstrate the efficiency of the method. Comparisons are made with plate elements in NASTRAN and the high-precision plate element developed by Cowper and his co-workers. Some considerations are given to implementation of the constraint method into general purpose computer programs such as NASTRAN.

Boundary element method for 2D materials and thin films.

PubMed

Hrtoň, M; Křápek, V; Šikola, T

2017-10-02

2D materials emerge as a viable platform for the control of light at the nanoscale. In this context the need has arisen for a fast and reliable tool capable of capturing their strictly 2D nature in 3D light scattering simulations. So far, 2D materials and their patterned structures (ribbons, discs, etc.) have been mostly treated as very thin films of subnanometer thickness with an effective dielectric function derived from their 2D optical conductivity. In this study an extension to the existing framework of the boundary element method (BEM) with 2D materials treated as a conductive interface between two media is presented. The testing of our enhanced method on problems with known analytical solutions reveals that for certain types of tasks the new modification is faster than the original BEM algorithm. Furthermore, the representation of 2D materials as an interface allows us to simulate problems in which their optical properties depend on spatial coordinates. Such spatial dependence can occur naturally or can be tailored artificially to attain new functional properties.

A DFFD simulation method combined with the spectral element method for solid-fluid-interaction problems

NASA Astrophysics Data System (ADS)

Chen, Li-Chieh; Huang, Mei-Jiau

2017-02-01

A 2D simulation method for a rigid body moving in an incompressible viscous fluid is proposed. It combines one of the immersed-boundary methods, the DFFD (direct forcing fictitious domain) method with the spectral element method; the former is employed for efficiently capturing the two-way FSI (fluid-structure interaction) and the geometric flexibility of the latter is utilized for any possibly co-existing stationary and complicated solid or flow boundary. A pseudo body force is imposed within the solid domain to enforce the rigid body motion and a Lagrangian mesh composed of triangular elements is employed for tracing the rigid body. In particular, a so called sub-cell scheme is proposed to smooth the discontinuity at the fluid-solid interface and to execute integrations involving Eulerian variables over the moving-solid domain. The accuracy of the proposed method is verified through an observed agreement of the simulation results of some typical flows with analytical solutions or existing literatures.

Hybrid Residual Flexibility/Mass-Additive Method for Structural Dynamic Testing

NASA Technical Reports Server (NTRS)

Tinker, M. L.

2003-01-01

A large fixture was designed and constructed for modal vibration testing of International Space Station elements. This fixed-base test fixture, which weighs thousands of pounds and is anchored to a massive concrete floor, initially utilized spherical bearings and pendulum mechanisms to simulate Shuttle orbiter boundary constraints for launch of the hardware. Many difficulties were encountered during a checkout test of the common module prototype structure, mainly due to undesirable friction and excessive clearances in the test-article-to-fixture interface bearings. Measured mode shapes and frequencies were not representative of orbiter-constrained modes due to the friction and clearance effects in the bearings. As a result, a major redesign effort for the interface mechanisms was undertaken. The total cost of the fixture design, construction and checkout, and redesign was over $2 million. Because of the problems experienced with fixed-base testing, alternative free-suspension methods were studied, including the residual flexibility and mass-additive approaches. Free-suspension structural dynamics test methods utilize soft elastic bungee cords and overhead frame suspension systems that are less complex and much less expensive than fixed-base systems. The cost of free-suspension fixturing is on the order of tens of thousands of dollars as opposed to millions, for large fixed-base fixturing. In addition, free-suspension test configurations are portable, allowing modal tests to be done at sites without modal test facilities. For example, a mass-additive modal test of the ASTRO-1 Shuttle payload was done at the Kennedy Space Center launch site. In this Technical Memorandum, the mass-additive and residual flexibility test methods are described in detail. A discussion of a hybrid approach that combines the best characteristics of each method follows and is the focus of the study.

A higher-order conservation element solution element method for solving hyperbolic differential equations on unstructured meshes

NASA Astrophysics Data System (ADS)

Bilyeu, David

This dissertation presents an extension of the Conservation Element Solution Element (CESE) method from second- to higher-order accuracy. The new method retains the favorable characteristics of the original second-order CESE scheme, including (i) the use of the space-time integral equation for conservation laws, (ii) a compact mesh stencil, (iii) the scheme will remain stable up to a CFL number of unity, (iv) a fully explicit, time-marching integration scheme, (v) true multidimensionality without using directional splitting, and (vi) the ability to handle two- and three-dimensional geometries by using unstructured meshes. This algorithm has been thoroughly tested in one, two and three spatial dimensions and has been shown to obtain the desired order of accuracy for solving both linear and non-linear hyperbolic partial differential equations. The scheme has also shown its ability to accurately resolve discontinuities in the solutions. Higher order unstructured methods such as the Discontinuous Galerkin (DG) method and the Spectral Volume (SV) methods have been developed for one-, two- and three-dimensional application. Although these schemes have seen extensive development and use, certain drawbacks of these methods have been well documented. For example, the explicit versions of these two methods have very stringent stability criteria. This stability criteria requires that the time step be reduced as the order of the solver increases, for a given simulation on a given mesh. The research presented in this dissertation builds upon the work of Chang, who developed a fourth-order CESE scheme to solve a scalar one-dimensional hyperbolic partial differential equation. The completed research has resulted in two key deliverables. The first is a detailed derivation of a high-order CESE methods on unstructured meshes for solving the conservation laws in two- and three-dimensional spaces. The second is the code implementation of these numerical methods in a computer code. For

Highly Accurate Beam Torsion Solutions Using the p-Version Finite Element Method

NASA Technical Reports Server (NTRS)

Smith, James P.

1996-01-01

A new treatment of the classical beam torsion boundary value problem is applied. Using the p-version finite element method with shape functions based on Legendre polynomials, torsion solutions for generic cross-sections comprised of isotropic materials are developed. Element shape functions for quadrilateral and triangular elements are discussed, and numerical examples are provided.

Self-supporting method; an alternative method for steel truss bridge element replacement

NASA Astrophysics Data System (ADS)

Arsyad, Muhammad; Sangadji, Senot; As'ad, Sholihin

2017-11-01

Steel truss bridge often requires replacement of its element due to serious damage caused by traffic accidents. This replacement is carried out using temporary supporting structure. It would be difficult when the available space for the temporary structure is quite limited and or the position of work is at a high elevation. The self-supporting method is proposed instead of temporary supporting structure. This paper will discuss an innovative method of bridge rehabilitation by utilizing the existing bridge structure. It requires such temporary connecting structure that installed on the existing bridge element, therefore, the forces during replacement process could be transferred to the bridge foundation directly. By taking the case on a steel truss bridge Jetis Salatiga which requires element replacement due to its damages on two main diagonals, a modeling is carried out to get a proper repair method. Structural analysis is conducted for three temporary connecting structure models: “I,” “V,” and triangular model. Stresses and translations that occur in the structure are used as constraints. Bridge bearings are modeled in two different modes: fixed-fixed system and fixed-free one. Temperature load is given in each condition to obtain the appropriate time for execution. The triangular model is chosen as the best one. In the fixed-fixed mode, this method can be carried out in a temperature range 27-28.8° C, while in fixed-free one, the temperature it is allowed between 27-43.4 °C. The D4 is dismantled first by cutting the D4 leaving an area of 1140.2 mm2 or 127 mm web length to enable plastic condition until the D4 collapses. At the beginning of elongation occurs, immediately performed a slowly jacking on a temporary connecting structure so that the force on D4 is gradually transferred to the temporary connecting structure then the D4 and D5 are set in their place.

Semi-automatic sparse preconditioners for high-order finite element methods on non-uniform meshes

NASA Astrophysics Data System (ADS)

Austin, Travis M.; Brezina, Marian; Jamroz, Ben; Jhurani, Chetan; Manteuffel, Thomas A.; Ruge, John

2012-05-01

High-order finite elements often have a higher accuracy per degree of freedom than the classical low-order finite elements. However, in the context of implicit time-stepping methods, high-order finite elements present challenges to the construction of efficient simulations due to the high cost of inverting the denser finite element matrix. There are many cases where simulations are limited by the memory required to store the matrix and/or the algorithmic components of the linear solver. We are particularly interested in preconditioned Krylov methods for linear systems generated by discretization of elliptic partial differential equations with high-order finite elements. Using a preconditioner like Algebraic Multigrid can be costly in terms of memory due to the need to store matrix information at the various levels. We present a novel method for defining a preconditioner for systems generated by high-order finite elements that is based on a much sparser system than the original high-order finite element system. We investigate the performance for non-uniform meshes on a cube and a cubed sphere mesh, showing that the sparser preconditioner is more efficient and uses significantly less memory. Finally, we explore new methods to construct the sparse preconditioner and examine their effectiveness for non-uniform meshes. We compare results to a direct use of Algebraic Multigrid as a preconditioner and to a two-level additive Schwarz method.

JWST center of curvature test method and results

NASA Astrophysics Data System (ADS)

Saif, Babak; Chaney, David; Greenfield, Perry; Van Gorkom, Kyle; Brooks, Keira; Hack, Warren; Bluth, Marcel; Bluth, Josh; Sanders, James; Smith, Koby; Carey, Larkin; Chaung, Sze; Keski-Kuha, Ritva; Feinberg, Lee; Tournois, Severine; Smith, W. Scott; Kradinov, Vladimir

2017-09-01

The James Webb Space Telescope (JWST) recently saw the completion of the assembly process for the Optical Telescope Element and Integrated Science Instrument Module (OTIS). This integration effort was performed at Goddard Space Flight Center (GSFC) in Greenbelt, Maryland. In conjunction with this assembly process a series of vibration and acoustic tests were performed. To help assure the telescope's primary mirror was not adversely impacted by this environmental testing an optical center of curvature (CoC) test was performed to measure changes in the mirror's optical performance. The primary is a 6.5 meter diameter mirror consisting of 18 individual hexagonal segments. Each segment is an off-axis asphere. There are a total of three prescriptions repeated six times each. As part of the CoC test each segment was individually measured using a high-speed interferometer (HSI) designed and built specifically for this test. This interferometer is capable of characterizing both static and dynamic characteristics of the mirrors. The latter capability was used, with the aid of a vibration stinger applying a low-level input force, to measure the dynamic characteristic changes of the PM backplane structure. This paper describes the CoC test setup, an innovative alignment method, and both static and dynamic test results.

Ignitability test method

NASA Technical Reports Server (NTRS)

Bement, Laurence J.; Schimmel, Morry L.

1989-01-01

To overcome serious weaknesses in determining the performance of initiating devices, a novel 'ignitability test method', representing actual design interfaces and ignition materials, has been developed. Ignition device output consists of heat, light, gas an burning particles. Past research methods have evaluated these parameters individually. This paper describes the development and demonstration of an ignitability test method combining all these parameters, and the quantitative assessment of the ignition performance of two widely used percussion primers, the M42C1-PA101 and the M42C2-793. The ignition materials used for this evaluation were several powder, granule and pellet sizes of black powder and boron-potassium nitrate. This test method should be useful for performance evaluation of all initiator types, quality assurance, evaluation of ignition interfaces, and service life studies of initiators and ignition materials.

Evaluation of a transfinite element numerical solution method for nonlinear heat transfer problems

NASA Technical Reports Server (NTRS)

Cerro, J. A.; Scotti, S. J.

1991-01-01

Laplace transform techniques have been widely used to solve linear, transient field problems. A transform-based algorithm enables calculation of the response at selected times of interest without the need for stepping in time as required by conventional time integration schemes. The elimination of time stepping can substantially reduce computer time when transform techniques are implemented in a numerical finite element program. The coupling of transform techniques with spatial discretization techniques such as the finite element method has resulted in what are known as transfinite element methods. Recently attempts have been made to extend the transfinite element method to solve nonlinear, transient field problems. This paper examines the theoretical basis and numerical implementation of one such algorithm, applied to nonlinear heat transfer problems. The problem is linearized and solved by requiring a numerical iteration at selected times of interest. While shown to be acceptable for weakly nonlinear problems, this algorithm is ineffective as a general nonlinear solution method.

A chemical-spectrochemical method for the determination of rare earth elements and thorium in cerium minerals

USGS Publications Warehouse

Rose, H.J.; Murata, K.J.; Carron, M.K.

1954-01-01

In a combined chemical-spectrochemical procedure for quantitatively determining rare earth elements in cerium minerals, cerium is determined volumetrically, a total rare earths plus thoria precipitate is separated chemically, the ceria content of the precipitate is raised to 80??0 percent by adding pure ceria, and the resulting mixture is analyzed for lanthanum, praseodymium, neodymium, samarium, gadolinium, yttrium, and thorium spectrochemically by means of the d.c. carbon arc. Spectral lines of singly ionized cerium are used as internal standard lines in the spectrochemical determination which is patterned after Fassel's procedure [1]. Results of testing the method with synthetic mixtures of rare earths and with samples of chemically analyzed cerium minerals show that the coefficient of variation for a quadruplicate determination of any element does not exceed 5??0 (excepting yttrium at concentrations less than 1 percent) and that the method is free of serious systematic error. ?? 1954.

Tree stability under wind: simulating uprooting with root breakage using a finite element method.

PubMed

Yang, Ming; Défossez, Pauline; Danjon, Frédéric; Fourcaud, Thierry

2014-09-01

Windstorms are the major natural hazard affecting European forests, causing tree damage and timber losses. Modelling tree anchorage mechanisms has progressed with advances in plant architectural modelling, but it is still limited in terms of estimation of anchorage strength. This paper aims to provide a new model for root anchorage, including the successive breakage of roots during uprooting. The model was based on the finite element method. The breakage of individual roots was taken into account using a failure law derived from previous work carried out on fibre metal laminates. Soil mechanical plasticity was considered using the Mohr-Coulomb failure criterion. The mechanical model for roots was implemented in the numerical code ABAQUS using beam elements embedded in a soil block meshed with 3-D solid elements. The model was tested by simulating tree-pulling experiments previously carried out on a tree of Pinus pinaster (maritime pine). Soil mechanical parameters were obtained from laboratory tests. Root system architecture was digitized and imported into ABAQUS while root material properties were estimated from the literature. Numerical simulations of tree-pulling tests exhibited realistic successive root breakages during uprooting, which could be seen in the resulting response curves. Broken roots could be visually located within the root system at any stage of the simulations. The model allowed estimation of anchorage strength in terms of the critical turning moment and accumulated energy, which were in good agreement with in situ measurements. This study provides the first model of tree anchorage strength for P. pinaster derived from the mechanical strength of individual roots. The generic nature of the model permits its further application to other tree species and soil conditions.

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Wang, Junping; Ye, Xiu

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for bothmore » the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.« less

A Least-Squares-Based Weak Galerkin Finite Element Method for Second Order Elliptic Equations

DOE PAGES

Mu, Lin; Wang, Junping; Ye, Xiu

2017-08-17

Here, in this article, we introduce a least-squares-based weak Galerkin finite element method for the second order elliptic equation. This new method is shown to provide very accurate numerical approximations for both the primal and the flux variables. In contrast to other existing least-squares finite element methods, this new method allows us to use discontinuous approximating functions on finite element partitions consisting of arbitrary polygon/polyhedron shapes. We also develop a Schur complement algorithm for the resulting discretization problem by eliminating all the unknowns that represent the solution information in the interior of each element. Optimal order error estimates for bothmore » the primal and the flux variables are established. An extensive set of numerical experiments are conducted to demonstrate the robustness, reliability, flexibility, and accuracy of the least-squares-based weak Galerkin finite element method. Finally, the numerical examples cover a wide range of applied problems, including singularly perturbed reaction-diffusion equations and the flow of fluid in porous media with strong anisotropy and heterogeneity.« less

Rolling element fatigue testing of gear materials

NASA Technical Reports Server (NTRS)

Nahm, A. H.

1978-01-01

Rolling element fatigue lives of nine alloys were evaluated in Rolling Contact (RC) rigs. Test conditions included a Hertzian stress at 4,826 MPa (700 ksi), a rolling speed of 6.23 m/sec (245 in/sec.). Tests were run with a Type I oil (MIL-L-7808G) at room temperature. B-10 lives (10% failure rate) of alloys were compared versus reference alloys, VIM-VAR AISI M-50 and VAR AISI 9310. Six case carburizing alloys (AISI 9310, CBS600, CBS1000M, EX00014, Vasco X-2 and EX00053) and three through-hardening alloys (AISI M-50, VascoMax 350 and Vasco Matrix 2 evaluated, showed RCF performance inferior or equivalent to that of AISI 9310 and AISI M-50. It was also found that the effects of vacuum melting processes, different tempering temperatures, freezing cycle during heat treating, shot peening, gold plating and chrome plating employed in the present investigation did not significantly affect RCF life.

System-Integrated Finite Element Analysis of a Full-Scale Helicopter Crash Test with Deployable Energy Absorbers

NASA Technical Reports Server (NTRS)

Annett, Martin S.; Polanco, Michael A.

2010-01-01

A full-scale crash test of an MD-500 helicopter was conducted in December 2009 at NASA Langley's Landing and Impact Research facility (LandIR). The MD-500 helicopter was fitted with a composite honeycomb Deployable Energy Absorber (DEA) and tested under vertical and horizontal impact velocities of 26-ft/sec and 40-ft/sec, respectively. The objectives of the test were to evaluate the performance of the DEA concept under realistic crash conditions and to generate test data for validation of a system integrated finite element model. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests was conducted to evaluate the impact performances of various components, including a new crush tube and the DEA blocks. Parameters defined within the system integrated finite element model were determined from these tests. The objective of this paper is to summarize the finite element models developed and analyses performed, beginning with pre-test predictions and continuing through post-test validation.

High-precision solution to the moving load problem using an improved spectral element method

NASA Astrophysics Data System (ADS)

Wen, Shu-Rui; Wu, Zhi-Jing; Lu, Nian-Li

2018-02-01

In this paper, the spectral element method (SEM) is improved to solve the moving load problem. In this method, a structure with uniform geometry and material properties is considered as a spectral element, which means that the element number and the degree of freedom can be reduced significantly. Based on the variational method and the Laplace transform theory, the spectral stiffness matrix and the equivalent nodal force of the beam-column element are established. The static Green function is employed to deduce the improved function. The proposed method is applied to two typical engineering practices—the one-span bridge and the horizontal jib of the tower crane. The results have revealed the following. First, the new method can yield extremely high-precision results of the dynamic deflection, the bending moment and the shear force in the moving load problem. In most cases, the relative errors are smaller than 1%. Second, by comparing with the finite element method, one can obtain the highly accurate results using the improved SEM with smaller element numbers. Moreover, the method can be widely used for statically determinate as well as statically indeterminate structures. Third, the dynamic deflection of the twin-lift jib decreases with the increase in the moving load speed, whereas the curvature of the deflection increases. Finally, the dynamic deflection, the bending moment and the shear force of the jib will all increase as the magnitude of the moving load increases.

The Mixed Finite Element Multigrid Method for Stokes Equations

PubMed Central

Muzhinji, K.; Shateyi, S.; Motsa, S. S.

2015-01-01

The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q 2-Q 1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results. PMID:25945361

Lichen elements as pollution indicators: evaluation of methods for large monitoring programmes

Treesearch

Susan Will-Wolf; Sarah Jovan; Michael C. Amacher

2017-01-01

Lichen element content is a reliable indicator for relative air pollution load in research and monitoring programmes requiring both efficiency and representation of many sites. We tested the value of costly rigorous field and handling protocols for sample element analysis using five lichen species. No relaxation of rigour was supported; four relaxed protocols generated...

Method for the removal of elemental mercury from a gas stream

DOEpatents

Mendelsohn, Marshall H.; Huang, Hann-Sheng

1999-01-01

A method is provided to remove elemental mercury from a gas stream by reacting the gas stream with an oxidizing solution to convert the elemental mercury to soluble mercury compounds. Other constituents are also oxidized. The gas stream is then passed through a wet scrubber to remove the mercuric compounds and oxidized constituents.

Finite Element Simulation of Three Full-Scale Crash Tests for Cessna 172 Aircraft

NASA Technical Reports Server (NTRS)

Mason, Brian H.; Warren, Jerry E., Jr.

2017-01-01

The NASA Emergency Locator Transmitter Survivability and Reliability (ELT-SAR) project was initiated in 2013 to assess the crash performance standards for the next generation of emergency locator transmitter (ELT) systems. Three Cessna 172 aircraft were acquired to perform crash testing at NASA Langley Research Center's Landing and Impact Research Facility. Full-scale crash tests were conducted in the summer of 2015 and each test article was subjected to severe, but survivable, impact conditions including a flare-to-stall during emergency landing, and two controlled-flight-into-terrain scenarios. Full-scale finite element analyses were performed using a commercial explicit solver, ABAQUS. The first test simulated impacting a concrete surface represented analytically by a rigid plane. Tests 2 and 3 simulated impacting a dirt surface represented analytically by an Eulerian grid of brick elements using a Mohr-Coulomb material model. The objective of this paper is to summarize the test and analysis results for the three full-scale crash tests. Simulation models of the airframe which correlate well with the tests are needed for future studies of alternate ELT mounting configurations.

Traditional lecture versus jigsaw learning method for teaching Medication Therapy Management (MTM) core elements.

PubMed

Wilson, Jennifer A; Pegram, Angela H; Battise, Dawn M; Robinson, April M

2017-11-01

To determine if traditional didactic lecture or the jigsaw learning method is more effective to teach the medication therapy management (MTM) core elements in a first year pharmacy course. Traditional didactic lecture and a pre-class reading assignment were used in the fall semester cohort, and the jigsaw method was used in the spring semester cohort. Jigsaw is a cooperative learning strategy requiring students to assume responsibility for learning, and subsequently teaching peers. The students were responsible for reading specific sections of the pre-class reading, and then teaching other students in small groups about their specific reading assignments. To assess potential differences, identical pre- and post-tests were administered before and after the MTM section. Additionally, grade performance on an in-class project and final exam questions were compared, and students were surveyed on perceptions of teaching method used. A total of 45 and 43 students completed both the pre- and post-test in the fall and spring (96% and 93% response rate), respectively. Improvement in post-test scores favored the traditional method (p = 0.001). No statistical differences were noted between groups with grade performance on the in-class project and final exam questions. However, students favored the jigsaw method over traditional lecture and perceived improvements in problem solving skills, listening/communication skills and encouragement of cooperative learning (p = 0.018, 0.025 and 0.031). Although students favored the jigsaw learning method, traditional didactic lecture was more effective for the pre- and post-knowledge test performance. This may indicate that traditional didactic lecture is more effective for more foundational content. Copyright © 2017 Elsevier Inc. All rights reserved.

Cognitive-graphic method for constructing of hierarchical forms of basic functions of biquadratic finite element

NASA Astrophysics Data System (ADS)

Astionenko, I. O.; Litvinenko, O. I.; Osipova, N. V.; Tuluchenko, G. Ya.; Khomchenko, A. N.

2016-10-01

Recently the interpolation bases of the hierarchical type have been used for the problem solving of the approximation of multiple arguments functions (such as in the finite-element method). In this work the cognitive graphical method of constructing of the hierarchical form bases on the serendipity finite elements is suggested, which allowed to get the alternative bases on a biquadratic finite element from the serendipity family without internal knots' inclusion. The cognitive-graphic method allowed to improve the known interpolation procedure of Taylor and to get the modified elements with irregular arrangement of knots. The proposed procedures are universal and are spread in the area of finite-elements.

On conforming mixed finite element methods for incompressible viscous flow problems

NASA Technical Reports Server (NTRS)

Gunzburger, M. D; Nicolaides, R. A.; Peterson, J. S.

1982-01-01

The application of conforming mixed finite element methods to obtain approximate solutions of linearized Navier-Stokes equations is examined. Attention is given to the convergence rates of various finite element approximations of the pressure and the velocity field. The optimality of the convergence rates are addressed in terms of comparisons of the approximation convergence to a smooth solution in relation to the best approximation available for the finite element space used. Consideration is also devoted to techniques for efficient use of a Gaussian elimination algorithm to obtain a solution to a system of linear algebraic equations derived by finite element discretizations of linear partial differential equations.

Stabilized finite element methods to simulate the conductances of ion channels

NASA Astrophysics Data System (ADS)

Tu, Bin; Xie, Yan; Zhang, Linbo; Lu, Benzhuo

2015-03-01

We have previously developed a finite element simulator, ichannel, to simulate ion transport through three-dimensional ion channel systems via solving the Poisson-Nernst-Planck equations (PNP) and Size-modified Poisson-Nernst-Planck equations (SMPNP), and succeeded in simulating some ion channel systems. However, the iterative solution between the coupled Poisson equation and the Nernst-Planck equations has difficulty converging for some large systems. One reason we found is that the NP equations are advection-dominated diffusion equations, which causes troubles in the usual FE solution. The stabilized schemes have been applied to compute fluids flow in various research fields. However, they have not been studied in the simulation of ion transport through three-dimensional models based on experimentally determined ion channel structures. In this paper, two stabilized techniques, the SUPG and the Pseudo Residual-Free Bubble function (PRFB) are introduced to enhance the numerical robustness and convergence performance of the finite element algorithm in ichannel. The conductances of the voltage dependent anion channel (VDAC) and the anthrax toxin protective antigen pore (PA) are simulated to validate the stabilization techniques. Those two stabilized schemes give reasonable results for the two proteins, with decent agreement with both experimental data and Brownian dynamics (BD) simulations. For a variety of numerical tests, it is found that the simulator effectively avoids previous numerical instability after introducing the stabilization methods. Comparison based on our test data set between the two stabilized schemes indicates both SUPG and PRFB have similar performance (the latter is slightly more accurate and stable), while SUPG is relatively more convenient to implement.

Fatigue assessment for selected connections of structural steel bridge components using the finite elements method

NASA Astrophysics Data System (ADS)

Śledziewski, Krzysztof

2018-01-01

Material fatigue it is one of the most frequent causes of steel bridge failures, particularly the bridges already existing. Thus, the procedure of fatigue life assessment is one of the most relevant procedures in a comprehensive assessment of load-carrying capacity and service life of the structure. A reliable assessment of the fatigue life is predominantly decisive for estimation of the remaining service life. Hitherto, calculation methods of welded joints took into account only stresses occurring in cross sections of whole elements and did not take into account stress concentration occurring in the vicinity of the weld, caused by geometrical aspects of the detail. At present, use of the Finite Element Analysis, makes possible looking for more accurate approach to the fatigue design of steel structures. The method of geometrical stresses is just such approach which is based on definition of stresses which take into account geometry of the detail. The study presents fatigue assessment of a representative type of welded joint in welded bridge structures. The testing covered longitudinal attachments. The main analyses were carried out on the basis of FEM and the method of local stresses, so-called "hot-spot" stresses. The obtained values of stresses were compared with the values obtained in accordance with the method of nominal stress.

Simulation and Flight Test Capability for Testing Prototype Sense and Avoid System Elements

NASA Technical Reports Server (NTRS)

Howell, Charles T.; Stock, Todd M.; Verstynen, Harry A.; Wehner, Paul J.

2012-01-01

NASA Langley Research Center (LaRC) and The MITRE Corporation (MITRE) have developed, and successfully demonstrated, an integrated simulation-to-flight capability for evaluating sense and avoid (SAA) system elements. This integrated capability consists of a MITRE developed fast-time computer simulation for evaluating SAA algorithms, and a NASA LaRC surrogate unmanned aircraft system (UAS) equipped to support hardware and software in-the-loop evaluation of SAA system elements (e.g., algorithms, sensors, architecture, communications, autonomous systems), concepts, and procedures. The fast-time computer simulation subjects algorithms to simulated flight encounters/ conditions and generates a fitness report that records strengths, weaknesses, and overall performance. Reviewed algorithms (and their fitness report) are then transferred to NASA LaRC where additional (joint) airworthiness evaluations are performed on the candidate SAA system-element configurations, concepts, and/or procedures of interest; software and hardware components are integrated into the Surrogate UAS research systems; and flight safety and mission planning activities are completed. Onboard the Surrogate UAS, candidate SAA system element configurations, concepts, and/or procedures are subjected to flight evaluations and in-flight performance is monitored. The Surrogate UAS, which can be controlled remotely via generic Ground Station uplink or automatically via onboard systems, operates with a NASA Safety Pilot/Pilot in Command onboard to permit safe operations in mixed airspace with manned aircraft. An end-to-end demonstration of a typical application of the capability was performed in non-exclusionary airspace in October 2011; additional research, development, flight testing, and evaluation efforts using this integrated capability are planned throughout fiscal year 2012 and 2013.

Method for the removal of elemental mercury from a gas stream

DOEpatents

Mendelsohn, M.H.; Huang, H.S.

1999-05-04

A method is provided to remove elemental mercury from a gas stream by reacting the gas stream with an oxidizing solution to convert the elemental mercury to soluble mercury compounds. Other constituents are also oxidized. The gas stream is then passed through a wet scrubber to remove the mercuric compounds and oxidized constituents. 7 figs.

Elemental X-ray mapping of agglutinated foraminifer tests: A non- destructive technique for determining compositional characteristics.

USGS Publications Warehouse

Commeau, R.F.; Reynolds, Leslie A.; Poag, C.W.

1985-01-01

The composition of agglutinated foraminiferal tests vary remarkably in response to local substrate characteristics, physiochemical properties of the water column and species- dependant selectivity of test components. We have employed a technique that combines a scanning electron microscope with an energy dispersive X-ray spectrometer system to identify major and minor elemental constituents of agglutinated foraminiferal walls. As a sample is bombarded with a beam of high energy electrons, X-rays are generated that are characteristic of the elements present. As a result, X- ray density maps can be produced for each of several elements present in the tests of agglutinated foraminifers. 

Coupling finite element and spectral methods: First results

NASA Technical Reports Server (NTRS)

Bernardi, Christine; Debit, Naima; Maday, Yvon

1987-01-01

A Poisson equation on a rectangular domain is solved by coupling two methods: the domain is divided in two squares, a finite element approximation is used on the first square and a spectral discretization is used on the second one. Two kinds of matching conditions on the interface are presented and compared. In both cases, error estimates are proved.

An Ellipsoidal Particle-Finite Element Method for Hypervelocity Impact Simulation. Chapter 1

NASA Technical Reports Server (NTRS)

Shivarama, Ravishankar; Fahrenthold, Eric P.

2004-01-01

A number of coupled particle-element and hybrid particle-element methods have been developed for the simulation of hypervelocity impact problems, to avoid certain disadvantages associated with the use of pure continuum based or pure particle based methods. To date these methods have employed spherical particles. In recent work a hybrid formulation has been extended to the ellipsoidal particle case. A model formulation approach based on Lagrange's equations, with particles entropies serving as generalized coordinates, avoids the angular momentum conservation problems which have been reported with ellipsoidal smooth particle hydrodynamics models.

Discontinuous dual-primal mixed finite elements for elliptic problems

NASA Technical Reports Server (NTRS)

Bottasso, Carlo L.; Micheletti, Stefano; Sacco, Riccardo

2000-01-01

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Naturally occurring levels of elements in fishes as determined by PIXE and XRF methods

NASA Astrophysics Data System (ADS)

Tallandini, L.; Giacobini, F.; Turchetto, M.; Galassini, S.; Liu, Q. X.; Shao, H. R.; Moschini, G.; Moro, R.; Gialanella, G.; Ghermandi, G.; Cecchi, R.; Injuk, J.; Valković, V.

1989-04-01

Naturally occurring levels of S, Cl, K, Ca, Cr, Mn, Fe, Ni, Cu, Zn, As, Se, Br, Sb, Sr and Pb were measured in the gills, liver and muscles of fishes ( Zosterisessor ophiocephalus Pall) in the northwestern region of the Adriatic Sea. The overall performance of PIXE and XRF methods was tested by the analysis of standard reference materials. The mean concentration values for elements were calculated from the distribution of experimentally determined concentration values. The obtained data are discussed in the framework of metal metabolism and toxicology.

Iterative methods for mixed finite element equations

NASA Technical Reports Server (NTRS)

Nakazawa, S.; Nagtegaal, J. C.; Zienkiewicz, O. C.

1985-01-01

Iterative strategies for the solution of indefinite system of equations arising from the mixed finite element method are investigated in this paper with application to linear and nonlinear problems in solid and structural mechanics. The augmented Hu-Washizu form is derived, which is then utilized to construct a family of iterative algorithms using the displacement method as the preconditioner. Two types of iterative algorithms are implemented. Those are: constant metric iterations which does not involve the update of preconditioner; variable metric iterations, in which the inverse of the preconditioning matrix is updated. A series of numerical experiments is conducted to evaluate the numerical performance with application to linear and nonlinear model problems.

Use of adjoint methods in the probabilistic finite element approach to fracture mechanics

NASA Technical Reports Server (NTRS)

Liu, Wing Kam; Besterfield, Glen; Lawrence, Mark; Belytschko, Ted

1988-01-01

The adjoint method approach to probabilistic finite element methods (PFEM) is presented. When the number of objective functions is small compared to the number of random variables, the adjoint method is far superior to the direct method in evaluating the objective function derivatives with respect to the random variables. The PFEM is extended to probabilistic fracture mechanics (PFM) using an element which has the near crack-tip singular strain field embedded. Since only two objective functions (i.e., mode I and II stress intensity factors) are needed for PFM, the adjoint method is well suited.

Adaptive mixed finite element methods for Darcy flow in fractured porous media

NASA Astrophysics Data System (ADS)

Chen, Huangxin; Salama, Amgad; Sun, Shuyu

2016-10-01

In this paper, we propose adaptive mixed finite element methods for simulating the single-phase Darcy flow in two-dimensional fractured porous media. The reduced model that we use for the simulation is a discrete fracture model coupling Darcy flows in the matrix and the fractures, and the fractures are modeled by one-dimensional entities. The Raviart-Thomas mixed finite element methods are utilized for the solution of the coupled Darcy flows in the matrix and the fractures. In order to improve the efficiency of the simulation, we use adaptive mixed finite element methods based on novel residual-based a posteriori error estimators. In addition, we develop an efficient upscaling algorithm to compute the effective permeability of the fractured porous media. Several interesting examples of Darcy flow in the fractured porous media are presented to demonstrate the robustness of the algorithm.

New formulation of the discrete element method

NASA Astrophysics Data System (ADS)

Rojek, Jerzy; Zubelewicz, Aleksander; Madan, Nikhil; Nosewicz, Szymon

2018-01-01

A new original formulation of the discrete element method based on the soft contact approach is presented in this work. The standard DEM has heen enhanced by the introduction of the additional (global) deformation mode caused by the stresses in the particles induced by the contact forces. Uniform stresses and strains are assumed for each particle. The stresses are calculated from the contact forces. The strains are obtained using an inverse constitutive relationship. The strains allow us to obtain deformed particle shapes. The deformed shapes (ellipses) are taken into account in contact detection and evaluation of the contact forces. A simple example of a uniaxial compression of a rectangular specimen, discreti.zed with equal sized particles is simulated to verify the DDEM algorithm. The numerical example shows that a particle deformation changes the particle interaction and the distribution of forces in the discrete element assembly. A quantitative study of micro-macro elastic properties proves the enhanced capabilities of the DDEM as compared to standard DEM.

A boundary element method for steady incompressible thermoviscous flow

NASA Technical Reports Server (NTRS)

Dargush, G. F.; Banerjee, P. K.

1991-01-01

A boundary element formulation is presented for moderate Reynolds number, steady, incompressible, thermoviscous flows. The governing integral equations are written exclusively in terms of velocities and temperatures, thus eliminating the need for the computation of any gradients. Furthermore, with the introduction of reference velocities and temperatures, volume modeling can often be confined to only a small portion of the problem domain, typically near obstacles or walls. The numerical implementation includes higher order elements, adaptive integration and multiregion capability. Both the integral formulation and implementation are discussed in detail. Several examples illustrate the high level of accuracy that is obtainable with the current method.

Comparison of the convolution quadrature method and enhanced inverse FFT with application in elastodynamic boundary element method

NASA Astrophysics Data System (ADS)

Schanz, Martin; Ye, Wenjing; Xiao, Jinyou

2016-04-01

Transient problems can often be solved with transformation methods, where the inverse transformation is usually performed numerically. Here, the discrete Fourier transform in combination with the exponential window method is compared with the convolution quadrature method formulated as inverse transformation. Both are inverse Laplace transforms, which are formally identical but use different complex frequencies. A numerical study is performed, first with simple convolution integrals and, second, with a boundary element method (BEM) for elastodynamics. Essentially, when combined with the BEM, the discrete Fourier transform needs less frequency calculations, but finer mesh compared to the convolution quadrature method to obtain the same level of accuracy. If further fast methods like the fast multipole method are used to accelerate the boundary element method the convolution quadrature method is better, because the iterative solver needs much less iterations to converge. This is caused by the larger real part of the complex frequencies necessary for the calculation, which improves the conditions of system matrix.

Characterization of human passive muscles for impact loads using genetic algorithm and inverse finite element methods.

PubMed

Chawla, A; Mukherjee, S; Karthikeyan, B

2009-02-01

The objective of this study is to identify the dynamic material properties of human passive muscle tissues for the strain rates relevant to automobile crashes. A novel methodology involving genetic algorithm (GA) and finite element method is implemented to estimate the material parameters by inverse mapping the impact test data. Isolated unconfined impact tests for average strain rates ranging from 136 s(-1) to 262 s(-1) are performed on muscle tissues. Passive muscle tissues are modelled as isotropic, linear and viscoelastic material using three-element Zener model available in PAMCRASH(TM) explicit finite element software. In the GA based identification process, fitness values are calculated by comparing the estimated finite element forces with the measured experimental forces. Linear viscoelastic material parameters (bulk modulus, short term shear modulus and long term shear modulus) are thus identified at strain rates 136 s(-1), 183 s(-1) and 262 s(-1) for modelling muscles. Extracted optimal parameters from this study are comparable with reported parameters in literature. Bulk modulus and short term shear modulus are found to be more influential in predicting the stress-strain response than long term shear modulus for the considered strain rates. Variations within the set of parameters identified at different strain rates indicate the need for new or improved material model, which is capable of capturing the strain rate dependency of passive muscle response with single set of material parameters for wide range of strain rates.

Simulation of Reacting Flow with a Discontinuous Spectral Element Method

NASA Astrophysics Data System (ADS)

Ghiasi, Zia; Mashayek, Farzad; Komperda, Jonathan

2013-11-01

While using high order methods is desirable in order to accurately capture the small scale mixing effects in reacting flows, the challenge is to develop and implement such methods for complex geometries. In this work, a high-order Discontinuous Spectral Element Method (DSEM) code, which solves for the Navier-Stokes equations, has been modified by adding the appropriate components to solve for scalar transport equations in order to simulate the chemical reaction. Dealing with discontinuous solution at element interfaces is a challenge that is met by patching the fluxes at mortars thus making them continuous on interfaces. The patching is performed using the Lax-Fredrichs numerical flux for scalars, whereas a generalized Riemann solver is used for the Navier-Stokes equations. Direct numerical simulation is conducted in a temporally developing mixing layer to validate the method for a single step reaction (F + rO --> [ 1 + r ] P). Next, the method is implemented to simulate a subsonic reacting flow in a slanted cavity combustor with gaseous fuel injectors to demonstrate the capability of the method to handle complex geometries. The results will be used for physical understanding of mixing and reaction in this type of combustors.

Research on Finite Element Model Generating Method of General Gear Based on Parametric Modelling

NASA Astrophysics Data System (ADS)

Lei, Yulong; Yan, Bo; Fu, Yao; Chen, Wei; Hou, Liguo

2017-06-01

Aiming at the problems of low efficiency and poor quality of gear meshing in the current mainstream finite element software, through the establishment of universal gear three-dimensional model, and explore the rules of unit and node arrangement. In this paper, a finite element model generation method of universal gear based on parameterization is proposed. Visual Basic program is used to realize the finite element meshing, give the material properties, and set the boundary / load conditions and other pre-processing work. The dynamic meshing analysis of the gears is carried out with the method proposed in this pape, and compared with the calculated values to verify the correctness of the method. The method greatly shortens the workload of gear finite element pre-processing, improves the quality of gear mesh, and provides a new idea for the FEM pre-processing.

Efficient Preconditioning for the p-Version Finite Element Method in Two Dimensions

DTIC Science & Technology

1989-10-01

paper, we study fast parallel preconditioners for systems of equations arising from the p-version finite element method. The p-version finite element...computations and the solution of a relatively small global auxiliary problem. We study two different methods. In the first (Section 3), the global...20], will be studied in the next section. Problem (3.12) is obviously much more easily solved than the original problem ,nd the procedure is highly

Analysis of Within-Test Variability of Non-Destructive Test Methods to Evaluate Compressive Strength of Normal Vibrated and Self-Compacting Concretes

NASA Astrophysics Data System (ADS)

Nepomuceno, Miguel C. S.; Lopes, Sérgio M. R.

2017-10-01

Non-destructive tests (NDT) have been used in the last decades for the assessment of in-situ quality and integrity of concrete elements. An important step in the application of NDT methods concerns to the interpretation and validation of the test results. In general, interpretation of NDT results should involve three distinct phases leading to the development of conclusions: processing of collected data, analysis of within-test variability and quantitative evaluation of property under investigation. The analysis of within-test variability can provide valuable information, since this can be compared with that of within-test variability associated with the NDT method in use, either to provide a measure of the quality control or to detect the presence of abnormal circumstances during the in-situ application. This paper reports the analysis of the experimental results of within-test variability of NDT obtained for normal vibrated concrete and self-compacting concrete. The NDT reported includes the surface hardness test, ultrasonic pulse velocity test, penetration resistance test, pull-off test, pull-out test and maturity test. The obtained results are discussed and conclusions are presented.

Verification of a non-hydrostatic dynamical core using the horizontal spectral element method and vertical finite difference method: 2-D aspects

NASA Astrophysics Data System (ADS)

Choi, S.-J.; Giraldo, F. X.; Kim, J.; Shin, S.

2014-11-01

The non-hydrostatic (NH) compressible Euler equations for dry atmosphere were solved in a simplified two-dimensional (2-D) slice framework employing a spectral element method (SEM) for the horizontal discretization and a finite difference method (FDM) for the vertical discretization. By using horizontal SEM, which decomposes the physical domain into smaller pieces with a small communication stencil, a high level of scalability can be achieved. By using vertical FDM, an easy method for coupling the dynamics and existing physics packages can be provided. The SEM uses high-order nodal basis functions associated with Lagrange polynomials based on Gauss-Lobatto-Legendre (GLL) quadrature points. The FDM employs a third-order upwind-biased scheme for the vertical flux terms and a centered finite difference scheme for the vertical derivative and integral terms. For temporal integration, a time-split, third-order Runge-Kutta (RK3) integration technique was applied. The Euler equations that were used here are in flux form based on the hydrostatic pressure vertical coordinate. The equations are the same as those used in the Weather Research and Forecasting (WRF) model, but a hybrid sigma-pressure vertical coordinate was implemented in this model. We validated the model by conducting the widely used standard tests: linear hydrostatic mountain wave, tracer advection, and gravity wave over the Schär-type mountain, as well as density current, inertia-gravity wave, and rising thermal bubble. The results from these tests demonstrated that the model using the horizontal SEM and the vertical FDM is accurate and robust provided sufficient diffusion is applied. The results with various horizontal resolutions also showed convergence of second-order accuracy due to the accuracy of the time integration scheme and that of the vertical direction, although high-order basis functions were used in the horizontal. By using the 2-D slice model, we effectively showed that the combined spatial

Direct tests of a pixelated microchannel plate as the active element of a shower maximum detector

DOE PAGES

Apresyan, A.; Los, S.; Pena, C.; ...

2016-05-07

One possibility to make a fast and radiation resistant shower maximum detector is to use a secondary emitter as an active element. We report our studies of microchannel plate photomultipliers (MCPs) as the active element of a shower-maximum detector. We present test beam results obtained using Photonis XP85011 to detect secondary particles of an electromagnetic shower. We focus on the use of the multiple pixels on the Photonis MCP in order to find a transverse two-dimensional shower distribution. A spatial resolution of 0.8 mm was obtained with an 8 GeV electron beam. As a result, a method for measuring themore » arrival time resolution for electromagnetic showers is presented, and we show that time resolution better than 40 ps can be achieved.« less

Direct tests of a pixelated microchannel plate as the active element of a shower maximum detector

DOE Office of Scientific and Technical Information (OSTI.GOV)

Apresyan, A.; Los, S.; Pena, C.

One possibility to make a fast and radiation resistant shower maximum detector is to use a secondary emitter as an active element. We report our studies of microchannel plate photomultipliers (MCPs) as the active element of a shower-maximum detector. We present test beam results obtained using Photonis XP85011 to detect secondary particles of an electromagnetic shower. We focus on the use of the multiple pixels on the Photonis MCP in order to find a transverse two-dimensional shower distribution. A spatial resolution of 0.8 mm was obtained with an 8 GeV electron beam. As a result, a method for measuring themore » arrival time resolution for electromagnetic showers is presented, and we show that time resolution better than 40 ps can be achieved.« less

Finite Element Simulation for Analysing the Design and Testing of an Energy Absorption System

PubMed Central

Segade, Abraham; López-Campos, José A.; Fernández, José R.; Casarejos, Enrique; Vilán, José A.

2016-01-01

It is not uncommon to use profiles to act as energy absorption parts in vehicle safety systems. This work analyses an impact attenuator based on a simple design and discusses the use of a thermoplastic material. We present the design of the impact attenuator and a mechanical test for the prototype. We develop a simulation model using the finite element method and explicit dynamics, and we evaluate the most appropriate mesh size and integration for describing the test results. Finally, we consider the performance of different materials, metallic ones (steel AISI 4310, Aluminium 5083-O) and a thermoplastic foam (IMPAXX500™). This reflects the car industry’s interest in using new materials to make high-performance, low-mass energy absorbers. We show the strength of the models when it comes to providing reliable results for large deformations and strong non-linearities, and how they are highly correlated with respect to the test results both in value and behaviour. PMID:28773778

Augmented Method to Improve Thermal Data for the Figure Drift Thermal Distortion Predictions of the JWST OTIS Cryogenic Vacuum Test

NASA Technical Reports Server (NTRS)

Park, Sang C.; Carnahan, Timothy M.; Cohen, Lester M.; Congedo, Cherie B.; Eisenhower, Michael J.; Ousley, Wes; Weaver, Andrew; Yang, Kan

2017-01-01

The JWST Optical Telescope Element (OTE) assembly is the largest optically stable infrared-optimized telescope currently being manufactured and assembled, and is scheduled for launch in 2018. The JWST OTE, including the 18 segment primary mirror, secondary mirror, and the Aft Optics Subsystem (AOS) are designed to be passively cooled and operate near 45K. These optical elements are supported by a complex composite backplane structure. As a part of the structural distortion model validation efforts, a series of tests are planned during the cryogenic vacuum test of the fully integrated flight hardware at NASA JSC Chamber A. The successful ends to the thermal-distortion phases are heavily dependent on the accurate temperature knowledge of the OTE structural members. However, the current temperature sensor allocations during the cryo-vac test may not have sufficient fidelity to provide accurate knowledge of the temperature distributions within the composite structure. A method based on an inverse distance relationship among the sensors and thermal model nodes was developed to improve the thermal data provided for the nanometer scale WaveFront Error (WFE) predictions. The Linear Distance Weighted Interpolation (LDWI) method was developed to augment the thermal model predictions based on the sparse sensor information. This paper will encompass the development of the LDWI method using the test data from the earlier pathfinder cryo-vac tests, and the results of the notional and as tested WFE predictions from the structural finite element model cases to characterize the accuracies of this LDWI method.

A Galleria Boundary Element Method for two-dimensional nonlinear magnetostatics

NASA Astrophysics Data System (ADS)

Brovont, Aaron D.

The Boundary Element Method (BEM) is a numerical technique for solving partial differential equations that is used broadly among the engineering disciplines. The main advantage of this method is that one needs only to mesh the boundary of a solution domain. A key drawback is the myriad of integrals that must be evaluated to populate the full system matrix. To this day these integrals have been evaluated using numerical quadrature. In this research, a Galerkin formulation of the BEM is derived and implemented to solve two-dimensional magnetostatic problems with a focus on accurate, rapid computation. To this end, exact, closed-form solutions have been derived for all the integrals comprising the system matrix as well as those required to compute fields in post-processing; the need for numerical integration has been eliminated. It is shown that calculation of the system matrix elements using analytical solutions is 15-20 times faster than with numerical integration of similar accuracy. Furthermore, through the example analysis of a c-core inductor, it is demonstrated that the present BEM formulation is a competitive alternative to the Finite Element Method (FEM) for linear magnetostatic analysis. Finally, the BEM formulation is extended to analyze nonlinear magnetostatic problems via the Dual Reciprocity Method (DRBEM). It is shown that a coarse, meshless analysis using the DRBEM is able to achieve RMS error of 3-6% compared to a commercial FEM package in lightly saturated conditions.

Mathematical aspects of finite element methods for incompressible viscous flows

NASA Technical Reports Server (NTRS)

Gunzburger, M. D.

1986-01-01

Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.

Discontinuous Galerkin finite element method for the nonlinear hyperbolic problems with entropy-based artificial viscosity stabilization

NASA Astrophysics Data System (ADS)

Zingan, Valentin Nikolaevich

This work develops a discontinuous Galerkin finite element discretization of non- linear hyperbolic conservation equations with efficient and robust high order stabilization built on an entropy-based artificial viscosity approximation. The solutions of equations are represented by elementwise polynomials of an arbitrary degree p > 0 which are continuous within each element but discontinuous on the boundaries. The discretization of equations in time is done by means of high order explicit Runge-Kutta methods identified with respective Butcher tableaux. To stabilize a numerical solution in the vicinity of shock waves and simultaneously preserve the smooth parts from smearing, we add some reasonable amount of artificial viscosity in accordance with the physical principle of entropy production in the interior of shock waves. The viscosity coefficient is proportional to the local size of the residual of an entropy equation and is bounded from above by the first-order artificial viscosity defined by a local wave speed. Since the residual of an entropy equation is supposed to be vanishingly small in smooth regions (of the order of the Local Truncation Error) and arbitrarily large in shocks, the entropy viscosity is almost zero everywhere except the shocks, where it reaches the first-order upper bound. One- and two-dimensional benchmark test cases are presented for nonlinear hyperbolic scalar conservation laws and the system of compressible Euler equations. These tests demonstrate the satisfactory stability properties of the method and optimal convergence rates as well. All numerical solutions to the test problems agree well with the reference solutions found in the literature. We conclude that the new method developed in the present work is a valuable alternative to currently existing techniques of viscous stabilization.

Correlation Results for a Mass Loaded Vehicle Panel Test Article Finite Element Models and Modal Survey Tests

NASA Technical Reports Server (NTRS)

Maasha, Rumaasha; Towner, Robert L.

2012-01-01

High-fidelity Finite Element Models (FEMs) were developed to support a recent test program at Marshall Space Flight Center (MSFC). The FEMs correspond to test articles used for a series of acoustic tests. Modal survey tests were used to validate the FEMs for five acoustic tests (a bare panel and four different mass-loaded panel configurations). An additional modal survey test was performed on the empty test fixture (orthogrid panel mounting fixture, between the reverb and anechoic chambers). Modal survey tests were used to test-validate the dynamic characteristics of FEMs used for acoustic test excitation. Modal survey testing and subsequent model correlation has validated the natural frequencies and mode shapes of the FEMs. The modal survey test results provide a basis for the analysis models used for acoustic loading response test and analysis comparisons

Evaluation of Strip Footing Bearing Capacity Built on the Anthropogenic Embankment by Random Finite Element Method

NASA Astrophysics Data System (ADS)

Pieczynska-Kozlowska, Joanna

2014-05-01

One of a geotechnical problem in the area of Wroclaw is an anthropogenic embankment layer delaying to the depth of 4-5m, arising as a result of historical incidents. In such a case an assumption of bearing capacity of strip footing might be difficult. The standard solution is to use a deep foundation or foundation soil replacement. However both methods generate significant costs. In the present paper the authors focused their attention on the influence of anthropogenic embankment variability on bearing capacity. Soil parameters were defined on the basis of CPT test and modeled as 2D anisotropic random fields and the assumption of bearing capacity were made according deterministic finite element methods. Many repeated of the different realizations of random fields lead to stable expected value of bearing capacity. The algorithm used to estimate the bearing capacity of strip footing was the random finite element method (e.g. [1]). In traditional approach of bearing capacity the formula proposed by [2] is taken into account. qf = c'Nc + qNq + 0.5γBN- γ (1) where: qf is the ultimate bearing stress, cis the cohesion, qis the overburden load due to foundation embedment, γ is the soil unit weight, Bis the footing width, and Nc, Nq and Nγ are the bearing capacity factors. The method of evaluation the bearing capacity of strip footing based on finite element method incorporate five parameters: Young's modulus (E), Poisson's ratio (ν), dilation angle (ψ), cohesion (c), and friction angle (φ). In the present study E, ν and ψ are held constant while c and φ are randomized. Although the Young's modulus does not affect the bearing capacity it governs the initial elastic response of the soil. Plastic stress redistribution is accomplished using a viscoplastic algorithm merge with an elastic perfectly plastic (Mohr - Coulomb) failure criterion. In this paper a typical finite element mesh was assumed with 8-node elements consist in 50 columns and 20 rows. Footings width B

Steady and Unsteady Nozzle Simulations Using the Conservation Element and Solution Element Method

NASA Technical Reports Server (NTRS)

Friedlander, David Joshua; Wang, Xiao-Yen J.

2014-01-01

This paper presents results from computational fluid dynamic (CFD) simulations of a three-stream plug nozzle. Time-accurate, Euler, quasi-1D and 2D-axisymmetric simulations were performed as part of an effort to provide a CFD-based approach to modeling nozzle dynamics. The CFD code used for the simulations is based on the space-time Conservation Element and Solution Element (CESE) method. Steady-state results were validated using the Wind-US code and a code utilizing the MacCormack method while the unsteady results were partially validated via an aeroacoustic benchmark problem. The CESE steady-state flow field solutions showed excellent agreement with solutions derived from the other methods and codes while preliminary unsteady results for the three-stream plug nozzle are also shown. Additionally, a study was performed to explore the sensitivity of gross thrust computations to the control surface definition. The results showed that most of the sensitivity while computing the gross thrust is attributed to the control surface stencil resolution and choice of stencil end points and not to the control surface definition itself.Finally, comparisons between the quasi-1D and 2D-axisymetric solutions were performed in order to gain insight on whether a quasi-1D solution can capture the steady and unsteady nozzle phenomena without the cost of a 2D-axisymmetric simulation. Initial results show that while the quasi-1D solutions are similar to the 2D-axisymmetric solutions, the inability of the quasi-1D simulations to predict two dimensional phenomena limits its accuracy.

Test of tree core sampling for screening of toxic elements in soils from a Norwegian site.

PubMed

Algreen, Mette; Rein, Arno; Legind, Charlotte N; Amundsen, Carl Einar; Karlson, Ulrich Gosewinkel; Trapp, Stefan

2012-04-01

Tree core samples have been used to delineate organic subsurface plumes. In 2009 and 2010, samples were taken at trees growing on a former dump site in Norway and analyzed for arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), nickel (Ni), and zinc (Zn). Concentrations in wood were in averages (dw) 30 mg/kg for Zn, 2 mg/kg for Cu, and < 1 mg/kg for Cd, Cr, As and Ni. The concentrations in wood samples from the polluted test site were compared to those derived from a reference site. For all except one case, mean concentrations from the test site were higher than those from the reference site, but the difference was small and not always significant. Differences between tree species were usually higher than differences between reference and test site. Furthermore, all these elements occur naturally, and Cu, Ni, and Zn are essential minerals. Thus, all trees will have a natural background of these elements, and the occurrence alone does not indicate soil pollution. For the interpretation of the results, a comparison to wood samples from an unpolluted reference site with same species and similar soil conditions is required. This makes the tree core screening method less reliable for heavy metals than, e.g., for chlorinated solvents.

IMS Infrasound Monitoring Test Site at Trafelberg, Austria: Preliminary Results from In-Situ Response Measurements of Infrasound Elements

DTIC Science & Technology

2010-09-01

IMS INFRASOUND MONITORING TEST SITE AT TRAFELBERG, AUSTRIA: PRELIMINARY RESULTS FROM IN-SITU RESPONSE MEASUREMENTS OF INFRASOUND ELEMENTS Thomas...The International Monitoring System (IMS) of the Comprehensive Nuclear Test-Ban Treaty Organization (CTBTO) has a unique infrasound test site...all four infrasound elements at I99 were made from 0.008 Hz to several Hz. For all four elements, from 0.01 to 0.1 Hz, the measured magnitude and

ICP-MS: Analytical Method for Identification and Detection of Elemental Impurities.

PubMed

Mittal, Mohini; Kumar, Kapil; Anghore, Durgadas; Rawal, Ravindra K

2017-01-01

Aim of this article is to review and discuss the currently used quantitative analytical method ICP-MS, which is used for quality control of pharmaceutical products. ICP-MS technique has several applications such as determination of single elements, multi element analysis in synthetic drugs, heavy metals in environmental water, trace element content of selected fertilizers and dairy manures. ICP-MS is also used for determination of toxic and essential elements in different varieties of food samples and metal pollutant present in the environment. The pharmaceuticals may generate impurities at various stages of development, transportation and storage which make them risky to be administered. Thus, it is essential that these impurities must be detected and quantified. ICP-MS plays an important function in the recognition and revealing of elemental impurities. Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.org.

Mixed-Mode Bending Method for Delamination Testing

NASA Technical Reports Server (NTRS)

Reeder, James R.; Crews, John R., Jr.

1990-01-01

A mixed mode delamination test procedure was developed combining double cantilever beam (DCB) mode I loading and end-notch fixture (ENF) mode II loading on a split unidirectional laminate. By loading with a lever, a single applied load simultaneously produces mode I and mode II bending loads on the specimen. This mixed-mode bending (MMB) test was analyzed using both finite-element procedures and beam theory to calculate the mode I and mode II components of strain-energy release rate G(sub I) and G(sub II), respectively. A wide range of G(sub I)/G(sub II) ratios can be produced by varying the load position on the lever. As the delamination extended, the G(sub I)/G(sub II) ratios varied by less than 5%. Beam theory equations agreed closely with the finite-element results and provide a basis for selection of G(sub I)/G(sub II) test ratios and a basis for computing the mode I and mode II components of measured delamination toughness. The MMB test was demonstrated using AS4/PEEK (APC2) unidirectional laminates. The MMB test introduced in this paper is rather simple and is believed to offer several advantages over most current mixed-mode test.

On mathematical modelling of aeroelastic problems with finite element method

NASA Astrophysics Data System (ADS)

Sváček, Petr

2018-06-01

This paper is interested in solution of two-dimensional aeroelastic problems. Two mathematical models are compared for a benchmark problem. First, the classical approach of linearized aerodynamical forces is described to determine the aeroelastic instability and the aeroelastic response in terms of frequency and damping coefficient. This approach is compared to the coupled fluid-structure model solved with the aid of finite element method used for approximation of the incompressible Navier-Stokes equations. The finite element approximations are coupled to the non-linear motion equations of a flexibly supported airfoil. Both methods are first compared for the case of small displacement, where the linearized approach can be well adopted. The influence of nonlinearities for the case of post-critical regime is discussed.

LEAKAGE TESTING METHOD

DOEpatents

McAdams, Wm.A.; Foss, M.H.

1958-08-12

A method of testing containers for leaks is described, particularly the testing of containers or cans in which the uranium slugs for nuelear reactors are jacketed. This method involves the immersion of the can in water under l50 pounds of pressure, then removing, drying, and coating the can with anhydrous copper sulfate. Amy water absorbed by the can under pressure will exude and discolor the copper sulfate in the area about the leak.

Spectral-element Method for 3D Marine Controlled-source EM Modeling

NASA Astrophysics Data System (ADS)

Liu, L.; Yin, C.; Zhang, B., Sr.; Liu, Y.; Qiu, C.; Huang, X.; Zhu, J.

2017-12-01

As one of the predrill reservoir appraisal methods, marine controlled-source EM (MCSEM) has been widely used in mapping oil reservoirs to reduce risk of deep water exploration. With the technical development of MCSEM, the need for improved forward modeling tools has become evident. We introduce in this paper spectral element method (SEM) for 3D MCSEM modeling. It combines the flexibility of finite-element and high accuracy of spectral method. We use Galerkin weighted residual method to discretize the vector Helmholtz equation, where the curl-conforming Gauss-Lobatto-Chebyshev (GLC) polynomials are chosen as vector basis functions. As a kind of high-order complete orthogonal polynomials, the GLC have the characteristic of exponential convergence. This helps derive the matrix elements analytically and improves the modeling accuracy. Numerical 1D models using SEM with different orders show that SEM method delivers accurate results. With increasing SEM orders, the modeling accuracy improves largely. Further we compare our SEM with finite-difference (FD) method for a 3D reservoir model (Figure 1). The results show that SEM method is more effective than FD method. Only when the mesh is fine enough, can FD achieve the same accuracy of SEM. Therefore, to obtain the same precision, SEM greatly reduces the degrees of freedom and cost. Numerical experiments with different models (not shown here) demonstrate that SEM is an efficient and effective tool for MSCEM modeling that has significant advantages over traditional numerical methods.This research is supported by Key Program of National Natural Science Foundation of China (41530320), China Natural Science Foundation for Young Scientists (41404093), and Key National Research Project of China (2016YFC0303100, 2017YFC0601900).

The use of Galerkin finite-element methods to solve mass-transport equations

USGS Publications Warehouse

Grove, David B.

1977-01-01

The partial differential equation that describes the transport and reaction of chemical solutes in porous media was solved using the Galerkin finite-element technique. These finite elements were superimposed over finite-difference cells used to solve the flow equation. Both convection and flow due to hydraulic dispersion were considered. Linear and Hermite cubic approximations (basis functions) provided satisfactory results: however, the linear functions were computationally more efficient for two-dimensional problems. Successive over relaxation (SOR) and iteration techniques using Tchebyschef polynomials were used to solve the sparce matrices generated using the linear and Hermite cubic functions, respectively. Comparisons of the finite-element methods to the finite-difference methods, and to analytical results, indicated that a high degree of accuracy may be obtained using the method outlined. The technique was applied to a field problem involving an aquifer contaminated with chloride, tritium, and strontium-90. (Woodard-USGS)

METHOD OF FORMING A FUEL ELEMENT FOR A NUCLEAR REACTOR

DOEpatents

Layer, E.H. Jr.; Peet, C.S.

1962-01-23

A method is given for preparing a fuel element for a nuclear reactor. The method includes the steps of sandblasting a body of uranium dioxide to roughen the surface thereof, depositing a thin layer of carbon thereon by thermal decomposition of methane, and cladding the uranium dioxide body with zirconium by gas pressure bonding. (AEC)

Modeling the mechanics of axonal fiber tracts using the embedded finite element method.

PubMed

Garimella, Harsha T; Kraft, Reuben H

2017-05-01

A subject-specific human head finite element model with embedded axonal fiber tractography obtained from diffusion tensor imaging was developed. The axonal fiber tractography finite element model was coupled with the volumetric elements in the head model using the embedded element method. This technique enables the calculation of axonal strains and real-time tracking of the mechanical response of the axonal fiber tracts. The coupled model was then verified using pressure and relative displacement-based (between skull and brain) experimental studies and was employed to analyze a head impact, demonstrating the applicability of this method in studying axonal injury. Following this, a comparison study of different injury criteria was performed. This model was used to determine the influence of impact direction on the extent of the axonal injury. The results suggested that the lateral impact loading is more dangerous compared to loading in the sagittal plane, a finding in agreement with previous studies. Through this analysis, we demonstrated the viability of the embedded element method as an alternative numerical approach for studying axonal injury in patient-specific human head models. Copyright © 2016 John Wiley & Sons, Ltd.

A simple finite element method for the Stokes equations

DOE PAGES

Mu, Lin; Ye, Xiu

2017-03-21

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.

A simple finite element method for the Stokes equations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mu, Lin; Ye, Xiu

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.

Wheat mill stream properties for discrete element method modeling

USDA-ARS?s Scientific Manuscript database

A discrete phase approach based on individual wheat kernel characteristics is needed to overcome the limitations of previous statistical models and accurately predict the milling behavior of wheat. As a first step to develop a discrete element method (DEM) model for the wheat milling process, this s...

A mesoscopic bridging scale method for fluids and coupling dissipative particle dynamics with continuum finite element method

PubMed Central

Kojic, Milos; Filipovic, Nenad; Tsuda, Akira

2012-01-01

A multiscale procedure to couple a mesoscale discrete particle model and a macroscale continuum model of incompressible fluid flow is proposed in this study. We call this procedure the mesoscopic bridging scale (MBS) method since it is developed on the basis of the bridging scale method for coupling molecular dynamics and finite element models [G.J. Wagner, W.K. Liu, Coupling of atomistic and continuum simulations using a bridging scale decomposition, J. Comput. Phys. 190 (2003) 249–274]. We derive the governing equations of the MBS method and show that the differential equations of motion of the mesoscale discrete particle model and finite element (FE) model are only coupled through the force terms. Based on this coupling, we express the finite element equations which rely on the Navier–Stokes and continuity equations, in a way that the internal nodal FE forces are evaluated using viscous stresses from the mesoscale model. The dissipative particle dynamics (DPD) method for the discrete particle mesoscale model is employed. The entire fluid domain is divided into a local domain and a global domain. Fluid flow in the local domain is modeled with both DPD and FE method, while fluid flow in the global domain is modeled by the FE method only. The MBS method is suitable for modeling complex (colloidal) fluid flows, where continuum methods are sufficiently accurate only in the large fluid domain, while small, local regions of particular interest require detailed modeling by mesoscopic discrete particles. Solved examples – simple Poiseuille and driven cavity flows illustrate the applicability of the proposed MBS method. PMID:23814322

Verification and application of the Iosipescu shear test method

NASA Technical Reports Server (NTRS)

Walrath, D. E.; Adams, D. F.

1984-01-01

Finite element models were used to study the effects of notch angle variations on the stress state within an Iosipescu shear test speciment. These analytical results were also studied to determine the feasibility of using strain gage rosettes and a modified extensometer to measure shear strains in this test specimen. Analytical results indicate that notch angle variations produced only small differences in simulated shear properties. Both strain gage rosettes and the modified extensometer were shown to be feasible shear strain transducers for the test method. The Iosipoescu shear test fixture was redesigned to incorporate several improvements. These improvements include accommodation of a 50 percent larger specimen for easier measurement of shear train, a clamping mechanism to relax strict tolerances on specimen width, and a self contained alignment tool for use during specimen installation. A set of in-plane and interlaminar shear properties were measured for three graphite fabric/epoxy composites of T300/934 composite material. The three weave patterns were Oxford, 5-harness satin, and 8-harness satin.

Validation of finite element and boundary element methods for predicting structural vibration and radiated noise

NASA Technical Reports Server (NTRS)

Seybert, A. F.; Wu, X. F.; Oswald, Fred B.

1992-01-01

Analytical and experimental validation of methods to predict structural vibration and radiated noise are presented. A rectangular box excited by a mechanical shaker was used as a vibrating structure. Combined finite element method (FEM) and boundary element method (BEM) models of the apparatus were used to predict the noise radiated from the box. The FEM was used to predict the vibration, and the surface vibration was used as input to the BEM to predict the sound intensity and sound power. Vibration predicted by the FEM model was validated by experimental modal analysis. Noise predicted by the BEM was validated by sound intensity measurements. Three types of results are presented for the total radiated sound power: (1) sound power predicted by the BEM modeling using vibration data measured on the surface of the box; (2) sound power predicted by the FEM/BEM model; and (3) sound power measured by a sound intensity scan. The sound power predicted from the BEM model using measured vibration data yields an excellent prediction of radiated noise. The sound power predicted by the combined FEM/BEM model also gives a good prediction of radiated noise except for a shift of the natural frequencies that are due to limitations in the FEM model.

Passive interior noise reduction analysis of King Air 350 turboprop aircraft using boundary element method/finite element method (BEM/FEM)

NASA Astrophysics Data System (ADS)

Dandaroy, Indranil; Vondracek, Joseph; Hund, Ron; Hartley, Dayton

2005-09-01

The objective of this study was to develop a vibro-acoustic computational model of the Raytheon King Air 350 turboprop aircraft with an intent to reduce propfan noise in the cabin. To develop the baseline analysis, an acoustic cavity model of the aircraft interior and a structural dynamics model of the aircraft fuselage were created. The acoustic model was an indirect boundary element method representation using SYSNOISE, while the structural model was a finite-element method normal modes representation in NASTRAN and subsequently imported to SYSNOISE. In the acoustic model, the fan excitation sources were represented employing the Ffowcs Williams-Hawkings equation. The acoustic and the structural models were fully coupled in SYSNOISE and solved to yield the baseline response of acoustic pressure in the aircraft interior and vibration on the aircraft structure due to fan noise. Various vibration absorbers, tuned to fundamental blade passage tone (100 Hz) and its first harmonic (200 Hz), were applied to the structural model to study their effect on cabin noise reduction. Parametric studies were performed to optimize the number and location of these passive devices. Effects of synchrophasing and absorptive noise treatments applied to the aircraft interior were also investigated for noise reduction.

Enhanced Low-Enriched Uranium Fuel Element for the Advanced Test Reactor

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pope, M. A.; DeHart, M. D.; Morrell, S. R.

2015-03-01

Under the current US Department of Energy (DOE) policy and planning scenario, the Advanced Test Reactor (ATR) and its associated critical facility (ATRC) will be reconfigured to operate on low-enriched uranium (LEU) fuel. This effort has produced a conceptual design for an Enhanced LEU Fuel (ELF) element. This fuel features monolithic U-10Mo fuel foils and aluminum cladding separated by a thin zirconium barrier. As with previous iterations of the ELF design, radial power peaking is managed using different U-10Mo foil thicknesses in different plates of the element. The lead fuel element design, ELF Mk1A, features only three fuel meat thicknesses,more » a reduction from the previous iterations meant to simplify manufacturing. Evaluation of the ELF Mk1A fuel design against reactor performance requirements is ongoing, as are investigations of the impact of manufacturing uncertainty on safety margins. The element design has been evaluated in what are expected to be the most demanding design basis accident scenarios and has met all initial thermal-hydraulic criteria.« less

Strength Analysis on Ship Ladder Using Finite Element Method

NASA Astrophysics Data System (ADS)

Budianto; Wahyudi, M. T.; Dinata, U.; Ruddianto; Eko P., M. M.

2018-01-01

In designing the ship’s structure, it should refer to the rules in accordance with applicable classification standards. In this case, designing Ladder (Staircase) on a Ferry Ship which is set up, it must be reviewed based on the loads during ship operations, either during sailing or at port operations. The classification rules in ship design refer to the calculation of the structure components described in Classification calculation method and can be analysed using the Finite Element Method. Classification Regulations used in the design of Ferry Ships used BKI (Bureau of Classification Indonesia). So the rules for the provision of material composition in the mechanical properties of the material should refer to the classification of the used vessel. The analysis in this structure used program structure packages based on Finite Element Method. By using structural analysis on Ladder (Ladder), it obtained strength and simulation structure that can withstand load 140 kg both in static condition, dynamic, and impact. Therefore, the result of the analysis included values of safety factors in the ship is to keep the structure safe but the strength of the structure is not excessive.

International Space Station Bacteria Filter Element Post-Flight Testing and Service Life Prediction

NASA Technical Reports Server (NTRS)

Perry, J. L.; von Jouanne, R. G.; Turner, E. H.

2003-01-01

The International Space Station uses high efficiency particulate air (HEPA) filters to remove particulate matter from the cabin atmosphere. Known as Bacteria Filter Elements (BFEs), there are 13 elements deployed on board the ISS's U.S. Segment. The pre-flight service life prediction of 1 year for the BFEs is based upon performance engineering analysis of data collected during developmental testing that used a synthetic dust challenge. While this challenge is considered reasonable and conservative from a design perspective, an understanding of the actual filter loading is required to best manage the critical ISS Program resources. Thus testing was conducted on BFEs returned from the ISS to refine the service life prediction. Results from this testing and implications to ISS resource management are discussed. Recommendations for realizing significant savings to the ISS Program are presented.

Combining existing numerical models with data assimilation using weighted least-squares finite element methods.

PubMed

Rajaraman, Prathish K; Manteuffel, T A; Belohlavek, M; Heys, Jeffrey J

2017-01-01

A new approach has been developed for combining and enhancing the results from an existing computational fluid dynamics model with experimental data using the weighted least-squares finite element method (WLSFEM). Development of the approach was motivated by the existence of both limited experimental blood velocity in the left ventricle and inexact numerical models of the same flow. Limitations of the experimental data include measurement noise and having data only along a two-dimensional plane. Most numerical modeling approaches do not provide the flexibility to assimilate noisy experimental data. We previously developed an approach that could assimilate experimental data into the process of numerically solving the Navier-Stokes equations, but the approach was limited because it required the use of specific finite element methods for solving all model equations and did not support alternative numerical approximation methods. The new approach presented here allows virtually any numerical method to be used for approximately solving the Navier-Stokes equations, and then the WLSFEM is used to combine the experimental data with the numerical solution of the model equations in a final step. The approach dynamically adjusts the influence of the experimental data on the numerical solution so that more accurate data are more closely matched by the final solution and less accurate data are not closely matched. The new approach is demonstrated on different test problems and provides significantly reduced computational costs compared with many previous methods for data assimilation. Copyright © 2016 John Wiley & Sons, Ltd. Copyright © 2016 John Wiley & Sons, Ltd.

Dissolution Behaviour of Metal Elements from Several Types of E-waste Using Leaching Test

NASA Astrophysics Data System (ADS)

Nor, Nik Hisyamudin Muhd; Amira Nordin, Nurul; Mohamad, Fariza; Jaibee, Shafizan; Ismail, Al Emran; Omar, Badrul; Fauzi Ahmad, Mohd; Rahim, Abd Khalil Abd; Kamaruddin, Muhamad Khalif Ikhwan Mohd; Turan, Faiz Mohd; Abu Bakar, Elmi; Yokoyama, Seiji

2017-08-01

Rapid development of the electrical and electronic was increasing annually due to the demand by the human being. Increasing production of electrical and electronic product led to the increasing of electric and electronic waste or can be called as the e-waste. The UN Environment Programme estimates that the world generates 20-50 million tons of the e-waste each year and the amount is raising three times faster than other forms of municipal waste. This study is focusing on the investigation of the dissolution behaviour of metal element from several types of e-waste by hydrometallurgical process. Leaching test was conducted on the e-waste by using acid as the reagent solution. Prior to the leaching test, manual dismantling, separation, and crushing process were carried out to the e-waste. The e-waste were characterized by Scanning Electron Microcopy (SEM) and the Energy Dispersive X-ray Spectroscopy (EDX) to define the elements inside the sample of e-waste. While the liquid residue from leaching test was analyzed by using Inductively Couple Plasma-Mass Spectrometer (ICP-MS) to define the dissolution behaviour of the metal element that contain in the e-waste. It was found that the longest time for dismantling process was the dismantling of laptop. The dissolution behaviour of Fe, Al, Zn and Pb elements in the e-waste has affected to the increase of pH. The increasing pH led to the reduction of the metals element during leaching process.

A Floating Node Method for the Modelling of Discontinuities Within a Finite Element

NASA Technical Reports Server (NTRS)

Pinho, Silvestre T.; Chen, B. Y.; DeCarvalho, Nelson V.; Baiz, P. M.; Tay, T. E.

2013-01-01

This paper focuses on the accurate numerical representation of complex networks of evolving discontinuities in solids, with particular emphasis on cracks. The limitation of the standard finite element method (FEM) in approximating discontinuous solutions has motivated the development of re-meshing, smeared crack models, the eXtended Finite Element Method (XFEM) and the Phantom Node Method (PNM). We propose a new method which has some similarities to the PNM, but crucially: (i) does not introduce an error on the crack geometry when mapping to natural coordinates; (ii) does not require numerical integration over only part of a domain; (iii) can incorporate weak discontinuities and cohesive cracks more readily; (iv) is ideally suited for the representation of multiple and complex networks of (weak, strong and cohesive) discontinuities; (v) leads to the same solution as a finite element mesh where the discontinuity is represented explicitly; and (vi) is conceptually simpler than the PNM.

A least-squares finite element method for 3D incompressible Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Lin, T. L.; Hou, Lin-Jun; Povinelli, Louis A.

1993-01-01

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.

Method for measuring recovery of catalytic elements from fuel cells

DOEpatents

Shore, Lawrence [Edison, NJ; Matlin, Ramail [Berkeley, NJ

2011-03-08

A method is provided for measuring the concentration of a catalytic clement in a fuel cell powder. The method includes depositing on a porous substrate at least one layer of a powder mixture comprising the fuel cell powder and an internal standard material, ablating a sample of the powder mixture using a laser, and vaporizing the sample using an inductively coupled plasma. A normalized concentration of catalytic element in the sample is determined by quantifying the intensity of a first signal correlated to the amount of catalytic element in the sample, quantifying the intensity of a second signal correlated to the amount of internal standard material in the sample, and using a ratio of the first signal intensity to the second signal intensity to cancel out the effects of sample size.

A pedagogical derivation of the matrix element method in particle physics data analysis

NASA Astrophysics Data System (ADS)

Sumowidagdo, Suharyo

2018-03-01

The matrix element method provides a direct connection between the underlying theory of particle physics processes and detector-level physical observables. I am presenting a pedagogically-oriented derivation of the matrix element method, drawing from elementary concepts in probability theory, statistics, and the process of experimental measurements. The level of treatment should be suitable for beginning research student in phenomenology and experimental high energy physics.

POSTPROCESSING MIXED FINITE ELEMENT METHODS FOR SOLVING CAHN-HILLIARD EQUATION: METHODS AND ERROR ANALYSIS

PubMed Central

Wang, Wansheng; Chen, Long; Zhou, Jie

2015-01-01

A postprocessing technique for mixed finite element methods for the Cahn-Hilliard equation is developed and analyzed. Once the mixed finite element approximations have been computed at a fixed time on the coarser mesh, the approximations are postprocessed by solving two decoupled Poisson equations in an enriched finite element space (either on a finer grid or a higher-order space) for which many fast Poisson solvers can be applied. The nonlinear iteration is only applied to a much smaller size problem and the computational cost using Newton and direct solvers is negligible compared with the cost of the linear problem. The analysis presented here shows that this technique remains the optimal rate of convergence for both the concentration and the chemical potential approximations. The corresponding error estimate obtained in our paper, especially the negative norm error estimates, are non-trivial and different with the existing results in the literatures. PMID:27110063

A projection hybrid high order finite volume/finite element method for incompressible turbulent flows

NASA Astrophysics Data System (ADS)

Busto, S.; Ferrín, J. L.; Toro, E. F.; Vázquez-Cendón, M. E.

2018-01-01

In this paper the projection hybrid FV/FE method presented in [1] is extended to account for species transport equations. Furthermore, turbulent regimes are also considered thanks to the k-ε model. Regarding the transport diffusion stage new schemes of high order of accuracy are developed. The CVC Kolgan-type scheme and ADER methodology are extended to 3D. The latter is modified in order to profit from the dual mesh employed by the projection algorithm and the derivatives involved in the diffusion term are discretized using a Galerkin approach. The accuracy and stability analysis of the new method are carried out for the advection-diffusion-reaction equation. Within the projection stage the pressure correction is computed by a piecewise linear finite element method. Numerical results are presented, aimed at verifying the formal order of accuracy of the scheme and to assess the performance of the method on several realistic test problems.

Finite element modeling of ROPS in static testing and rear overturns.

PubMed

Harris, J R; Mucino, V H; Etherton, J R; Snyder, K A; Means, K H

2000-08-01

Even with the technological advances of the last several decades, agricultural production remains one of the most hazardous occupations in the United States. Death due to tractor rollover is a prime contributor to this hazard. Standards for rollover protective structures (ROPS) performance and certification have been developed by groups such as the Society of Automotive Engineers (SAE) and the American Society of Agricultural Engineers (ASAE) to combat these problems. The current ROPS certification standard, SAE J2194, requires either a dynamic or static testing sequence or both. Although some ROPS manufacturers perform both the dynamic and static phases of SAE J2194 testing, it is possible for a ROPS to be certified for field operation using static testing alone. This research compared ROPS deformation response from a simulated SAE J2194 static loading sequence to ROPS deformation response as a result of a simulated rearward tractor rollover. Finite element analysis techniques for plastic deformation were used to simulate both the static and dynamic rear rollover scenarios. Stress results from the rear rollover model were compared to results from simulated static testing per SAE J2194. Maximum stress values from simulated rear rollovers exceeded maximum stress values recorded during simulated static testing for half of the elements comprising the uprights. In the worst case, the static model underpredicts dynamic model results by approximately 7%. In the best case, the static model overpredicts dynamic model results by approximately 32%. These results suggest the need for additional experimental work to characterize ROPS stress levels during staged overturns and during testing according to the SAE standard.

Simulation of thin slot spirals and dual circular patch antennas using the finite element method with mixed elements

NASA Technical Reports Server (NTRS)

Gong, Jian; Volakis, John L.; Nurnberger, Michael W.

1995-01-01

This semi-annual report describes progress up to mid-January 1995. The report contains five sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. Of significance is the development of decomposition schemes which separate the different regions of the antenna volume. Substantial effort was devoted to improving the feed model in the context of the finite element method (FEM). Finally, an innovative scheme for truncating finite element meshes is presented.

Development and Application of the p-Version of the Finite Element Method.

DTIC Science & Technology

1987-12-30

element method has been the subject of intensive study since the early 1950’s and perhaps even earlier. Study of the p-version of the finite element...method, on the other hand, began at *Washington University in St. Louis in the early 1970’s and led to a more recent study of the h-p version. Research...infinite strip to a bounded domain. 3.3 A Numerical Argument Principle In order to assure that all roots have indeed been obtained, we have studied the

The Matrix Element Method: Past, Present, and Future

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gainer, James S.; Lykken, Joseph; Matchev, Konstantin T.

2013-07-12

The increasing use of multivariate methods, and in particular the Matrix Element Method (MEM), represents a revolution in experimental particle physics. With continued exponential growth in computing capabilities, the use of sophisticated multivariate methods-- already common-- will soon become ubiquitous and ultimately almost compulsory. While the existence of sophisticated algorithms for disentangling signal and background might naively suggest a diminished role for theorists, the use of the MEM, with its inherent connection to the calculation of differential cross sections will benefit from collaboration between theorists and experimentalists. In this white paper, we will briefly describe the MEM and some ofmore » its recent uses, note some current issues and potential resolutions, and speculate about exciting future opportunities.« less

An optimized computational method for determining the beta dose distribution using a multiple-element thermoluminescent dosimeter system.

PubMed

Shen, L; Levine, S H; Catchen, G L

1987-07-01

This paper describes an optimization method for determining the beta dose distribution in tissue, and it describes the associated testing and verification. The method uses electron transport theory and optimization techniques to analyze the responses of a three-element thermoluminescent dosimeter (TLD) system. Specifically, the method determines the effective beta energy distribution incident on the dosimeter system, and thus the system performs as a beta spectrometer. Electron transport theory provides the mathematical model for performing the optimization calculation. In this calculation, parameters are determined that produce calculated doses for each of the chip/absorber components in the three-element TLD system. The resulting optimized parameters describe an effective incident beta distribution. This method can be used to determine the beta dose specifically at 7 mg X cm-2 or at any depth of interest. The doses at 7 mg X cm-2 in tissue determined by this method are compared to those experimentally determined using an extrapolation chamber. For a great variety of pure beta sources having different incident beta energy distributions, good agreement is found. The results are also compared to those produced by a commonly used empirical algorithm. Although the optimization method produces somewhat better results, the advantage of the optimization method is that its performance is not sensitive to the specific method of calibration.

METHOD OF PREPARING A CERAMIC FUEL ELEMENT

DOEpatents

Ross, W.T.; Bloomster, C.H.; Bardsley, R.E.

1963-09-01

A method is described for preparing a fuel element from -325 mesh PuO/ sub 2/ and -20 mesh UO/sub 2/, and the steps of screening --325 mesh UO/sub 2/ from the -20 mesh UO/sub 2/, mixing PuO/sub 2/ with the --325 mesh UO/sub 2/, blending this mixture with sufficient --20 mesh UO/sub 2/ to obtain the desired composition, introducing the blend into a metal tube, repeating the procedure until the tube is full, and vibrating the tube to compact the powder are included. (AEC)

Mixed-RKDG Finite Element Methods for the 2-D Hydrodynamic Model for Semiconductor Device Simulation

DOE PAGES

Chen, Zhangxin; Cockburn, Bernardo; Jerome, Joseph W.; ...

1995-01-01

In this paper we introduce a new method for numerically solving the equations of the hydrodynamic model for semiconductor devices in two space dimensions. The method combines a standard mixed finite element method, used to obtain directly an approximation to the electric field, with the so-called Runge-Kutta Discontinuous Galerkin (RKDG) method, originally devised for numerically solving multi-dimensional hyperbolic systems of conservation laws, which is applied here to the convective part of the equations. Numerical simulations showing the performance of the new method are displayed, and the results compared with those obtained by using Essentially Nonoscillatory (ENO) finite difference schemes. Frommore » the perspective of device modeling, these methods are robust, since they are capable of encompassing broad parameter ranges, including those for which shock formation is possible. The simulations presented here are for Gallium Arsenide at room temperature, but we have tested them much more generally with considerable success.« less

Architecting the Finite Element Method Pipeline for the GPU.

PubMed

Fu, Zhisong; Lewis, T James; Kirby, Robert M; Whitaker, Ross T

2014-02-01

The finite element method (FEM) is a widely employed numerical technique for approximating the solution of partial differential equations (PDEs) in various science and engineering applications. Many of these applications benefit from fast execution of the FEM pipeline. One way to accelerate the FEM pipeline is by exploiting advances in modern computational hardware, such as the many-core streaming processors like the graphical processing unit (GPU). In this paper, we present the algorithms and data-structures necessary to move the entire FEM pipeline to the GPU. First we propose an efficient GPU-based algorithm to generate local element information and to assemble the global linear system associated with the FEM discretization of an elliptic PDE. To solve the corresponding linear system efficiently on the GPU, we implement a conjugate gradient method preconditioned with a geometry-informed algebraic multi-grid (AMG) method preconditioner. We propose a new fine-grained parallelism strategy, a corresponding multigrid cycling stage and efficient data mapping to the many-core architecture of GPU. Comparison of our on-GPU assembly versus a traditional serial implementation on the CPU achieves up to an 87 × speedup. Focusing on the linear system solver alone, we achieve a speedup of up to 51 × versus use of a comparable state-of-the-art serial CPU linear system solver. Furthermore, the method compares favorably with other GPU-based, sparse, linear solvers.

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

NASA Technical Reports Server (NTRS)

Chung, T. J. (Editor); Karr, Gerald R. (Editor)

1989-01-01

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Finite Element Method Applied to Fuse Protection Design

NASA Astrophysics Data System (ADS)

Li, Sen; Song, Zhiquan; Zhang, Ming; Xu, Liuwei; Li, Jinchao; Fu, Peng; Wang, Min; Dong, Lin

2014-03-01

In a poloidal field (PF) converter module, fuse protection is of great importance to ensure the safety of the thyristors. The fuse is pre-selected in a traditional way and then verified by finite element analysis. A 3D physical model is built by ANSYS software to solve the thermal-electric coupled problem of transient process in case of external fault. The result shows that this method is feasible.

CFD Methods and Tools for Multi-Element Airfoil Analysis

NASA Technical Reports Server (NTRS)

Rogers, Stuart E.; George, Michael W. (Technical Monitor)

1995-01-01

This lecture will discuss the computational tools currently available for high-lift multi-element airfoil analysis. It will present an overview of a number of different numerical approaches, their current capabilities, short-comings, and computational costs. The lecture will be limited to viscous methods, including inviscid/boundary layer coupling methods, and incompressible and compressible Reynolds-averaged Navier-Stokes methods. Both structured and unstructured grid generation approaches will be presented. Two different structured grid procedures are outlined, one which uses multi-block patched grids, the other uses overset chimera grids. Turbulence and transition modeling will be discussed.

Unidirectional Fabric Drape Testing Method

PubMed Central

Mei, Zaihuan; Yang, Jingzhi; Zhou, Ting; Zhou, Hua

2015-01-01

In most cases, fabrics such as curtains, skirts, suit pants and so on are draped under their own gravity parallel to fabric plane while the gravity is perpendicular to fabric plane in traditional drape testing method. As a result, it does not conform to actual situation and the test data is not convincing enough. To overcome this problem, this paper presents a novel method which simulates the real mechanical conditions and ensures the gravity is parallel to the fabric plane. This method applied a low-cost Kinect Sensor device to capture the 3-dimensional (3D) drape profile, thus we obtained the drape degree parameters and aesthetic parameters by 3D reconstruction and image processing and analysis techniques. The experiment was conducted on our self-devised drape-testing instrument by choosing different kinds of weave structure fabrics as our testing samples and the results were compared with those of traditional method and subjective evaluation. Through regression and correlation analysis we found that this novel testing method was significantly correlated with the traditional and subjective evaluation method. We achieved a new, non-contact 3D measurement method for drape testing, namely unidirectional fabric drape testing method. This method is more suitable for evaluating drape behavior because it is more in line with actual mechanical conditions of draped fabrics and has a well consistency with the requirements of visual and aesthetic style of fabrics. PMID:26600387

Adaptive finite element method for turbulent flow near a propeller

NASA Astrophysics Data System (ADS)

Pelletier, Dominique; Ilinca, Florin; Hetu, Jean-Francois

1994-11-01

This paper presents an adaptive finite element method based on remeshing to solve incompressible turbulent free shear flow near a propeller. Solutions are obtained in primitive variables using a highly accurate finite element approximation on unstructured grids. Turbulence is modeled by a mixing length formulation. Two general purpose error estimators, which take into account swirl and the variation of the eddy viscosity, are presented and applied to the turbulent wake of a propeller. Predictions compare well with experimental measurements. The proposed adaptive scheme is robust, reliable and cost effective.

Load Testing, Finite Element Analysis, and Design of Steel Traffic-Signal Poles

DOT National Transportation Integrated Search

1993-07-01