Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mixed 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 181 records · Page 10

Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities such as line generation for general unstructured grids and optimization of the iterative solver. The efficiency of the new developments is assessed in a massively parallel environment. The assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Three benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere cylinder configuration, a transonic flow around the ONERA M6 wing, and a supersonic flow through a square duct are considered.

Mohagna J Pandya↗

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SCARE: A post-processor program to MSC/NASTRAN for the reliability analysis of structural ceramic components

A computer program was developed for calculating the statistical fast fracture reliability and failure probability of ceramic components. The program includes the two-parameter Weibull material fracture strength distribution model, using the principle of independent action for polyaxial stress states and Batdorf's shear-sensitive as well as shear-insensitive crack theories, all for volume distributed flaws in macroscopically isotropic solids. Both penny-shaped cracks and Griffith cracks are included in the Batdorf shear-sensitive crack response calculations, using Griffith's maximum tensile stress or critical coplanar strain energy release rate criteria to predict mixed mode fracture. Weibull material parameters can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and fracture data. The reliability prediction analysis uses MSC/NASTRAN stress, temperature and volume output, obtained from the use of three-dimensional, quadratic, isoparametric, or axisymmetric finite elements. The statistical fast fracture theories employed, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

SCARE - A postprocessor program to MSC/NASTRAN for reliability analysis of structural ceramic components

A computer program was developed for calculating the statistical fast fracture reliability and failure probability of ceramic components. The program includes the two-parameter Weibull material fracture strength distribution model, using the principle of independent action for polyaxial stress states and Batdorf's shear-sensitive as well as shear-insensitive crack theories, all for volume distributed flaws in macroscopically isotropic solids. Both penny-shaped cracks and Griffith cracks are included in the Batdorf shear-sensitive crack response calculations, using Griffith's maximum tensile stress or critical coplanar strain energy release rate criteria to predict mixed mode fracture. Weibull material parameters can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and fracture data. The reliability prediction analysis uses MSC/NASTRAN stress, temperature and volume output, obtained from the use of three-dimensional, quadratic, isoparametric, or axisymmetric finite elements. The statistical fast fracture theories employed, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

An analysis of delamination in angle-ply fiber reinforced composites

An analytical investigation of the mechanics and failure modes of delamination initiated from a service-induced surface flaw in angle-ply fiber reinforced composites is presented. The analysis employs a hybrid-stress finite element method including a crack-tip singular element with its field variables expressed by Muskhelishvili's complex stress functions. Solutions are obtained for delaminated composites with various laminate parameters. The results elucidate unique and important characteristics of delamination crack-tip response and interlaminar stress transfer mechanisms. Of particular interest are the mixed-mode stress intensity factors associated with the delamination crack. The influence of ply orientation on KI and KII and their effects on subsequent crack extension are discussed.

Wang, S. S.↗

An analysis of delamination in angle-ply fiber-reinforced composites

A study of the mechanics and failure modes of delamination initiated from a surface flaw in angle-ply fiber-reinforced composites is presented. The analysis employs a hybrid-stress finite-element method including a crack-tip singular element with its field variables expressed by Muskhelishvili's complex stress functions. Solutions are obtained for the delaminated composites with various laminate parameters. The results elucidate unique and important characteristics of delamination crack-tip response and interlaminar stress transfer mechanisms. Of particular interest are the mixed-mode stress-intensity factors associated with the delamination crack. The influence of ply orientation on KI and KII and their effects on subsequent crack extension are discussed.

Wang, S. S.↗

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Accuracy, Scalability, and Efficiency of Mixed-Element USM3D for Benchmark Three-Dimensional Flows

The unstructured, mixed-element, cell-centered, finite-volume flow solver USM3D is enhanced with new capabilities including parallelization, line generation for general unstructured grids, improved discretization scheme, and optimized iterative solver. The paper reports on the new developments to the flow solver and assesses the accuracy, scalability, and efficiency. The USM3D assessments are conducted using a baseline method and the recent hierarchical adaptive nonlinear iteration method framework. Two benchmark turbulent flows, namely, a subsonic separated flow around a three-dimensional hemisphere-cylinder configuration and a transonic flow around the ONERA M6 wing are considered.

Pandya, Mohagna J.↗

Scale–Resolving Simulations of a Fundamental Trailing–Edge Cooling Slot Using a Discontinuous–Galerkin Spectral–Element Method

The accurate prediction of turbulent mixing in high-pressure turbines that incorporate various airfoil surface-cooling strategies is becoming increasing critical to the design of modern gas turbine engines where the quest for improved efficiency is driving compressor overall pressure ratios and turbine inlet temperatures to much higher levels than ever before. In the present paper, a recently developed computational capability for accurate and efficient scale–resolving simulations of turbomachinery is extended to study the turbulent mixing mechanism of a simplified abstraction of an airfoil trailing-edge cooling slot - a plane wall jet with finite lip thickness discharging into an ambient flow. The computational capability is based on an entropy– stable, discontinuous–Galerkin approach that extends to arbitrarily high orders of spatial and temporal accuracy. The numerical results show that the present simulations capture the trends observed in the experiments. Discrepancies between the simulations and experiments are believed to be due to differences in the inflow profiles and tunnel side–wall effects. The thick lip configuration leads to a thicker wake and higher unsteadiness in the wall jet compared to the thin lip. A detailed comparison of the turbulent flowfields is presented to highlight differences arising due to lip thickness variations.

Scale Resolving↗

Mixed-mode fracture in unidirectional graphite epoxy composite laminates with central notch

Mixed-mode matrix fracture in central notched off-axis unidirectional composite laminates was investigated. A limited number of unidirectional tensile type specimens with a central, horizontal, notch were tested. Crack initiation and propagation were examined under various local stress fields that were controlled by fiber orientations. The tested specimens were simulated using a two dimensional finite element method with constant strain loading. The strain energy release rates along the crack were evaluated via crack closure technique. The variation of critical strain energy rates with off-axis angle was studied. The results from single (one-sided) and double (two-sided) crack simulations were presented and compared.

Binienda, Wieslaw K.↗

Re-analysis procedure based on the mixed formulation

An analysis procedure is presented for large-scale structural systems, with large numbers of degrees of freedom and design variables. The procedure uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. Other elements of the procedure include: (1) lumping the design variables into a single tracing parameter; (2) operator splitting or additive decomposition of different arrays in the finite element equations into the corresponding arrays of the original structure plus correction terms; and (3) application of a reduction method through the use of the finite element method and the classical Bubnov-Galerkin technique. The re-analysis procedure is applied to the linear static and free vibration problems of plate and shell structures.

Noor, Ahmed K.↗

Mixed-Mode Bending Method for Delamination Testing

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.

Reeder, James R.↗

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes

Cell-centered and node-centered approaches have been compared for unstructured finite-volume discretization of inviscid fluxes. The grids range from regular grids to irregular grids, including mixed-element grids and grids with random perturbations of nodes. Accuracy, complexity, and convergence rates of defect-correction iterations are studied for eight nominally second-order accurate schemes: two node-centered schemes with weighted and unweighted least-squares (LSQ) methods for gradient reconstruction and six cell-centered schemes two node-averaging with and without clipping and four schemes that employ different stencils for LSQ gradient reconstruction. The cell-centered nearest-neighbor (CC-NN) scheme has the lowest complexity; a version of the scheme that involves smart augmentation of the LSQ stencil (CC-SA) has only marginal complexity increase. All other schemes have larger complexity; complexity of node-centered (NC) schemes are somewhat lower than complexity of cell-centered node-averaging (CC-NA) and full-augmentation (CC-FA) schemes. On highly anisotropic grids typical of those encountered in grid adaptation, discretization errors of five of the six cell-centered schemes converge with second order on all tested grids; the CC-NA scheme with clipping degrades solution accuracy to first order. The NC schemes converge with second order on regular and/or triangular grids and with first order on perturbed quadrilaterals and mixed-element grids. All schemes may produce large relative errors in gradient reconstruction on grids with perturbed nodes. Defect-correction iterations for schemes employing weighted least-square gradient reconstruction diverge on perturbed stretched grids. Overall, the CC-NN and CC-SA schemes offer the best options of the lowest complexity and secondorder discretization errors. On anisotropic grids over a curved body typical of turbulent flow simulations, the discretization errors converge with second order and are small for the CC-NN, CC-SA, and CC-FA schemes on all grids and for NC schemes on triangular grids; the discretization errors of the CC-NA scheme without clipping do not converge on irregular grids. Accurate gradient reconstruction can be achieved by introducing a local approximate mapping; without approximate mapping, only the NC scheme with weighted LSQ method provides accurate gradients. Defect correction iterations for the CC-NA scheme without clipping diverge; for the NC scheme with weighted LSQ method, the iterations either diverge or converge very slowly. The best option in curved geometries is the CC-SA scheme that offers low complexity, second-order discretization errors, and fast convergence.

Diskin, Boris↗

A Benchmark Example for Delamination Growth Predictions Based on the Single Leg Bending Specimen Under Fatigue Loading

Analysis benchmarking is used to evaluate new algorithms for automated VCCT-based delamination growth analysis. First, existing benchmark cas s based on the Single Leg Bending (SLB) specimen for crack propagation prediction under quasi-static loading are summarized. Second, the development of new SLB-based benchmark cases to assess the static and fatigue growth prediction capabilities under mixed-mode I/II conditions is discussed in detail. Additionally, a scheme is proposed to interpolate between known fatigue delamination growth rates to obtain values for mixed-mode ratios for which data has not been defined in the input. Further, a comparison is presented, in which the benchmark cases are used to assess new analysis tools in ABAQUS/Standard FD03. These recently implemented tools yield results that are in good agreement with the benchmark examples. The ability to assess the implementation of new methods in one finite element code illustrates the value of establishing benchmark solutions.

Krueger, Ronald↗

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

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.

Raju, I. S.↗

A decentralized linear quadratic control design method for flexible structures

A decentralized suboptimal linear quadratic control design procedure which combines substructural synthesis, model reduction, decentralized control design, subcontroller synthesis, and controller reduction is proposed for the design of reduced-order controllers for flexible structures. The procedure starts with a definition of the continuum structure to be controlled. An evaluation model of finite dimension is obtained by the finite element method. Then, the finite element model is decomposed into several substructures by using a natural decomposition called substructuring decomposition. Each substructure, at this point, still has too large a dimension and must be reduced to a size that is Riccati-solvable. Model reduction of each substructure can be performed by using any existing model reduction method, e.g., modal truncation, balanced reduction, Krylov model reduction, or mixed-mode method. Then, based on the reduced substructure model, a subcontroller is designed by an LQ optimal control method for each substructure independently. After all subcontrollers are designed, a controller synthesis method called substructural controller synthesis is employed to synthesize all subcontrollers into a global controller. The assembling scheme used is the same as that employed for the structure matrices. Finally, a controller reduction scheme, called the equivalent impulse response energy controller (EIREC) reduction algorithm, is used to reduce the global controller to a reasonable size for implementation. The EIREC reduced controller preserves the impulse response energy of the full-order controller and has the property of matching low-frequency moments and low-frequency power moments. An advantage of the substructural controller synthesis method is that it relieves the computational burden associated with dimensionality. Besides that, the SCS design scheme is also a highly adaptable controller synthesis method for structures with varying configuration, or varying mass and stiffness properties.

Su, Tzu-Jeng↗