Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

Nonlinear optimization for a tumor invasion PDE model

In this work, we introduce a methodology to approximate unknown parameters that appear on a non-linear reaction–diffusion model of tumor invasion. These equations consider that tumor-induced alteration of micro-environmental pH furnishes a mechanism for cancer invasion. A coupled system reaction–diffusion explaining this model is given by three partial differential equations for the non-dimensional spatial distribution and temporal evolution of the density of normal tissue, the neoplastic tissue growth and the excess concentration of H{sup +} ions. The tumor model parameters have a corresponding biological meaning: the reabsorption rate, the destructive influence of H{sup +} ions in the healthy tissue, the growth rate of tumor tissue and the diffusion coefficient. We propose to solve the direct problem using the Finite Element Method (FEM) and minimize an appropriate functional including both the real data (obtained via in-vitro experiments and fluorescence ratio imaging microscopy) and the numerical solution. The gradient of the functional is computed by the adjoint method.

60 APPLIED LIFE SCIENCES↗

Design of voice coil motor dynamic focusing unit for a laser scanner

Laser scanning systems have been used for material processing tasks such as welding, cutting, marking, and drilling. However, applications have been limited by the small range of motion and slow speed of the focusing unit, which carries the focusing optics. To overcome these limitations, a dynamic focusing system with a long travel range and high speed is needed. In this study, a dynamic focusing unit for a laser scanning system with a voice coil motor (VCM) mechanism is proposed to enable fast speed and a wide focusing range. The VCM has finer precision and higher speed than conventional step motors and a longer travel range than earlier lead zirconium titanate actuators. The system has a hollow configuration to provide a laser beam path. This also makes it compact and transmission-free and gives it low inertia. The VCM's magnetics are modeled using a permeance model. Its design parameters are determined by optimization using the Broyden–Fletcher–Goldfarb–Shanno method and a sequential quadratic programming algorithm. After the VCM is designed, the dynamic focusing unit is fabricated and assembled. The permeance model is verified by a magnetic finite element method simulation tool, Maxwell 2D and 3D, and by measurement data from a gauss meter. The performance is verified experimentally. The results show a resolution of 0.2 μm and travel range of 16 mm. These are better than those of conventional focusing systems; therefore, this focusing unit can be applied to laser scanning systems for good machining capability.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High order curvilinear finite elements for elastic–plastic Lagrangian dynamics

This paper presents a high-order finite element method for calculating elastic–plastic flow on moving curvilinear meshes and is an extension of our general high-order curvilinear finite element approach for solving the Euler equations of gas dynamics in a Lagrangian frame [1,2]. In order to handle transition to plastic flow, we formulate the stress–strain relation in rate (or incremental) form and augment our semi-discrete equations for Lagrangian hydrodynamics with an additional evolution equation for the deviatoric stress which is valid for arbitrary order spatial discretizations of the kinematic and thermodynamic variables. The semi-discrete equation for the deviatoric stress rate is developed for 2D planar, 2D axisymmetric and full 3D geometries. For each case, the strain rate is approximated via a collocation method at zone quadrature points while the deviatoric stress is approximated using an L{sub 2} projection onto the thermodynamic basis. We apply high order, energy conserving, explicit time stepping methods to the semi-discrete equations to develop the fully discrete method. We conclude with numerical results from an extensive series of verification tests that demonstrate several practical advantages of using high-order finite elements for elastic–plastic flow.

97 MATHEMATICS AND COMPUTING↗

Gear Tooth Root Stresses of a Very Heavily Loaded Gear Pair-Case Study: Orbiter Body Flap Actuator Pinion and Ring Gear

The space shuttle orbiter's body flap actuator gearing was assessed as a case study of the stresses for very heavily loaded external-internal gear pairs (meshing pinion and ring gear). For many applications, using the high point of single tooth contact (HPSTC) to locate the position of the tooth force is adequate for assessing the maximum tooth root stress condition. But for aerospace gearing such an approach may be inadequate for assessing the stress condition while also simultaneously minimizing mass. In this work specialized contact analyses and finite element methods were used to study gear tooth stresses of body flap actuator gears. The analytical solutions considered the elastic deformations as an inherent part of the solutions. The ratio for the maximum tooth stresses using the HPSTC approach solutions relative to the contact analysis and finite element solutions were 1.40 for the ring gear and 1.28 for the pinion gear.

finite element method↗

SU-E-J-161: Inverse Problems for Optical Parameters in Laser Induced Thermal Therapy

Purpose: Magnetic resonance-guided laser-induced thermal therapy (MRgLITT) is investigated as a neurosurgical intervention for oncological applications throughout the body by active post market studies. Real-time MR temperature imaging is used to monitor ablative thermal delivery in the clinic. Additionally, brain MRgLITT could improve through effective planning for laser fiber's placement. Mathematical bioheat models have been extensively investigated but require reliable patient specific physical parameter data, e.g. optical parameters. This abstract applies an inverse problem algorithm to characterize optical parameter data obtained from previous MRgLITT interventions. Methods: The implemented inverse problem has three primary components: a parameter-space search algorithm, a physics model, and training data. First, the parameter-space search algorithm uses a gradient-based quasi-Newton method to optimize the effective optical attenuation coefficient, μ-eff. A parameter reduction reduces the amount of optical parameter-space the algorithm must search. Second, the physics model is a simplified bioheat model for homogeneous tissue where closed-form Green's functions represent the exact solution. Third, the training data was temperature imaging data from 23 MRgLITT oncological brain ablations (980 nm wavelength) from seven different patients. Results: To three significant figures, the descriptive statistics for μ-eff were 1470 m{sup −1} mean, 1360 m{sup −1} median, 369 m{sup −1} standard deviation, 933 m{sup −1} minimum and 2260 m{sup −1} maximum. The standard deviation normalized by the mean was 25.0%. The inverse problem took <30 minutes to optimize all 23 datasets. Conclusion: As expected, the inferred average is biased by underlying physics model. However, the standard deviation normalized by the mean is smaller than literature values and indicates an increased precision in the characterization of the optical parameters needed to plan MRgLITT procedures. This investigation demonstrates the potential for the optimization and validation of more sophisticated bioheat models that incorporate the uncertainty of the data into the predictions, e.g. stochastic finite element methods.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

On controlling the electronic states of shallow donors using a finite-size metal gate

The effect of an external electric field on the states of a shallow donor near a semiconductor surface is numerically simulated. A disk-shaped metal gate is considered as an electric-field source. The wavefunctions and energies of bound states are determined by the finite-element method. The critical characteristics of electron relocation between the donor and gate are determined for various gate diameters and boundary conditions, taking into account dielectric mismatch. The empirical dependences of these characteristics on the geometrical parameters and semiconductor properties are obtained. A simple trial function is proposed, which can be used to calculate the critical parameters using the Ritz variational method.

36 MATERIALS SCIENCE↗

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

TU-H-BRA-04: A Novel Superconducting Magnet Design for Optimized Patient Access and Minimal SSD for Use in a Linac-MR Hybrid

Purpose: A prototype rotating hybrid MR imaging system and linac has been developed to allow for simultaneous imaging and radiation delivery parallel to B{sub 0}. However, the design of a compact magnet capable of rotation in a small vault with sufficient patient access and a typical clinical source-to-surface distance (SSD) is challenging. This work presents a novel superconducting magnet design that allows for a reduced SSD and ample patient access by moving the superconducting coils to the side of the yoke. The yoke and pole-plate structures are shaped to direct the magnetic flux appropriately. Methods: The surface of the pole plate for the magnet assembly is optimized. The magnetic field calculations required in this work are performed with the 3D finite element method software package Opera-3D. Each tentative design strategy is virtually modeled in this software package and externally controlled by MATLAB, with its key geometries defined as variables. The particle swarm optimization algorithm is used to optimize the variables subject to the minimization of a cost function. At each iteration, Opera-3D will solve the magnetic field solution over a field-of-view suitable for MR imaging and the degree of field uniformity will be assessed to calculate the value of the cost function associated with that iteration. Results: An optimized magnet assembly that generates a homogenous 0.2T magnetic field over an ellipsoid with large axis of 30 cm and small axes of 20 cm is obtained. Conclusion: The distinct features of this model are the minimal distance between the yoke’s top and the isocentre and the improved patient access. On the other hand, having homogeneity over an ellipsoid give us a larger field-of-view, essential for geometric accuracy of the MRI system. The increase of B{sub 0} from 0.2T in the present model to 0.5T is the subject of future work. Funding Sources: Alberta Innovates - Health Solutions (AIHS)| Disclosure and Conflict of Interest: B. Gino Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta biplanar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

Salamon, N. J.↗

Immersed Boundary Method Halo Exchange in a Hemodynamics Application

In recent years, highly parallelized simulations of blood flow resolving individual blood cells have been demonstrated. Simulating such dense suspensions of deformable particles in flow often involves a partitioned fluid-structure interaction (FSI) algorithm, with separate solvers for Eulerian fluid and Lagrangian cell grids, plus a solver - e.g., immersed boundary method - for their interaction. Managing data motion in parallel FSI implementations is increasingly important, particularly for inhomogeneous systems like vascular geometries. In this study, we evaluate the influence of Eulerian and Lagrangian halo exchanges on efficiency and scalability of a partitioned FSI algorithm for blood flow. We describe an MPI+OpenMP implementation of the immersed boundary method coupled with lattice Boltzmann and finite element methods. We consider how communication and recomputation costs influence the optimization of halo exchanges with respect to three factors: immersed boundary interaction distance, cell suspension density, and relative fluid/cell solver costs.

Gounley, John↗

Thermally induced cracking on the massive concrete structure of the NSLS II synchrotron and its engineering remediation

Synchrotron accelerator facilities such as the NSLS II require extreme stability, both transient (short-term) and quasi-static (long-term) to achieve the desired resolution performance. Consequently, even μm-level movements, particularly differential movements between locations in the concrete structure supporting the accelerator electron beam lattice (storage ring) or high sensitivity experiments (experimental floor) will lead to serious degradation of its performance. Differential settlement in the overall structure or structural movement exceeding anticipated levels will inevitably degrade the performance and will require intervention. Presented in this paper are the design philosophy of the NSLS II ring structure favoring a monolithic ring, the observed cracking behavior of the young NSLS II concrete following casting and in combination with extreme ambient temperature fluctuation, the results of a non-linear, high-fidelity numerical analysis used to emulate the observed cracking and establish the driving mechanism, the numerical analysis-based identification of the crack-arresting solution and finally the implementation of the remediation solution and the long-term performance of the adopted engineering solution. The multi-stage process revealed that computational methods such as non-linear finite element methods have the potential of providing engineering guidance even when complex structures and in combination with non-linear materials, such as steel reinforcement and concrete are involved.

42 ENGINEERING↗

A data-driven approach to modeling cancer cell mechanics during microcirculatory transport

In order to understand the effect of cellular level features on the transport of circulating cancer cells in the microcirculation, there has been an increasing reliance on high-resolution in silico models. Accurate simulation of cancer cells flowing with blood cells requires resolving cellular-scale interactions in 3D, which is a significant computational undertaking warranting a cancer cell model that is both computationally efficient yet sufficiently complex to capture relevant behavior. Given that the characteristics of metastatic spread are known to depend on cancer type, it is crucial to account for mechanistic behavior representative of a specific cancer’s cells. To address this gap, in the present work we develop and validate a means by which an efficient and popular membrane model-based approach can be used to simulate deformable cancer cells and reproduce experimental data from specific cell lines. Here, cells are modeled using the immersed boundary method (IBM) within a lattice Boltzmann method (LBM) fluid solver, and the finite element method (FEM) is used to model cell membrane resistance to deformation. Through detailed comparisons with experiments, we (i) validate this model to represent cancer cells undergoing large deformation, (ii) outline a systematic approach to parameterize different cell lines to optimally fit experimental data over a range of deformations, and (iii) provide new insight into nucleated vs. non-nucleated cell models and their ability to match experiments. While many works have used the membrane-model based method employed here to model generic cancer cells, no quantitative comparisons with experiments exist in the literature for specific cell lines undergoing large deformation. Here, we describe a phenomenological, data-driven approach that can not only yield good agreement for large deformations, but explicitly detail how it can be used to represent different cancer cell lines. This model is readily incorporated into cell-resolved hemodynamic transport simulations, and thus offers significant potential to complement experiments towards providing new insights into various aspects of cancer progression.

59 BASIC BIOLOGICAL SCIENCES↗

Numerical experiments with the extended virtual element method for the Laplace problem with strong discontinuities

The virtual element method (VEM) is a stabilized Galerkin formulation on arbitrary polytopal meshes. In the VEM, the basis functions are implicit (virtual) — they are not known nor do they need to be computed within the problem domain. Suitable projection operators are used to decompose the bilinear form at the element-level into two parts: a consistent term that reproduces a given polynomial space and a correction term that ensures stability. In this study, we consider a low-order extended virtual element method (X-VEM) that is in the spirit of the extended finite element method for crack problems. Herein, we focus on the two-dimensional Laplace crack problem. In the X-VEM, we enrich the standard virtual element space with additional discontinuous functions through the framework of partition-of-unity. The nodal basis functions in the VEM are chosen as the partition-of-unity functions, and we study means to stabilize the standard and enriched sub-matrices that constitute the element stiffness matrix. Numerical experiments are performed on the problem of a cracked membrane under mode III loading to affirm the accuracy, and to establish the optimal convergence in energy of the method.

97 MATHEMATICS AND COMPUTING↗

Finite-element analysis of a weld-penetration problem

The stress concentration factor for a weld penetration defect is calculated by the finite-element method. A stress intensity factor is computed by use of the finite-element solution and the J-integral. The results are compared with experimental results.

Rogge, T. R.↗

Scattering of Pruppacher-Pitter raindrops at 30 GHz

Optimum design of modern ground-satellite communication systems requires the knowledge of rain-induced differential attenuation, differential phase shift, and cross polarization factors. Different available analytical techniques for raindrop scattering problems were assessed. These include: (1) geometrical theory of diffraction; (2) method of moment; (3) perturbation method; (4) point matching methods; (5) extended boundary condition method; and (6) global-local finite element method. The advantages and disadvantages of each are listed. The extended boundary condition method, which was determined to yield the most scattering results, is summarized. The scattered fields for Pruppacher-Pitter raindrops with sizes ranging from 0.5 mm to 3.5 mm at 20 C and at 30 GHz for several incidence angles are tabulated.

Yeh, C.↗

On the use of harmonic expansions in magnetic field problems in NASTRAN

The use of a prolate spheroidal harmonic expansion to compute far field results in magnetics problems with the DTNSRDC version of NASTRAN is demonstrated. When field problems with infinite exterior domains are solved by the finite element method, the domain must be modeled to infinity. The density of the finite element mesh usually decreases as the distance from the structure increases, thus degrading the results in the far field.

Hurwitz, M. M.↗

Analysis of structures with rotating, flexible substructures applied to rotorcraft aeroelasticity in GRASP

Application to the General Rotorcraft Aeromechanical Stability Program (GRASP) of new methodology for structural dynamic analysis, including substructuring, frames of reference, nodes, finite elements, and constraints, is discussed. The structure is decomposed into a hierarchy of substructures, and discrete relative motion between substructures is analyzed exactly. The finite element method is used to treat deformation of continua, and the library of finite elements includes a nonlinear beam element incorporating aeroelastic effects. Analytical bases for the aeroelastic beam element and the screw constraint are considered, and the important role of geometric stiffness in the formulation is shown.

Hodges, Dewey H.↗