Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixed finite element”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Theoretical and numerical difficulties in 3-D vector potential methods in finite element magnetostatic computations

This paper describes the results of application of three well known 3D magnetic vector potential (MVP) based finite element formulations for computation of magnetostatic fields in electrical devices. The three methods were identically applied to three practical examples, the first of which contains only one medium (free space), while the second and third examples contained a mix of free space and iron. The first of these methods is based on the unconstrained curl-curl of the MVP, while the second and third methods are predicated upon constraining the divergence of the MVP 10 zero (Coulomb's Gauge). It was found that the two latter methods cease to give useful and meaningful results when the global solution region contains a mix of media of high and low permeabilities. Furthermore, it was found that their results do not achieve the intended zero constraint on the divergence of the MVP.

Demerdash, N. A.↗
NASA NTRS2003

Modeling the Interactions Between Multiple Crack Closure Mechanisms at Threshold

A fatigue crack closure model is developed that includes interactions between the three closure mechanisms most likely to occur at threshold; plasticity, roughness, and oxide. This model, herein referred to as the CROP model (for Closure, Roughness, Oxide, and Plasticity), also includes the effects of out-of plane cracking and multi-axial loading. These features make the CROP closure model uniquely suited for, but not limited to, threshold applications. Rough cracks are idealized here as two-dimensional sawtooths, whose geometry induces mixed-mode crack- tip stresses. Continuum mechanics and crack-tip dislocation concepts are combined to relate crack face displacements to crack-tip loads. Geometric criteria are used to determine closure loads from crack-face displacements. Finite element results, used to verify model predictions, provide critical information about the locations where crack closure occurs.

Newman, John A.↗
NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗
NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗
NASA NTRS1992

Finite element Euler calculations of unsteady transonic cascade flows

A Galerkin finite element procedure incorporating an explicit Runge-Kutta time-stepping scheme has been developed in this work to solve unsteady transonic flow in cascades. The computational domain is discretized by a globally unstructured but locally structured blade-fitted deformable mesh. The Galerkin approximation is applied to the unsteady Euler equations based on a mixed Eulerian-Lagrangian description. The semi-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. An artificial dissipation operator of the type proposed by Jameson is adapted in the current scheme to capture shocks and suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. Results for both steady and unsteady transonic flows through cascades are presented and compared to existing finite volume solutions.

Hsiao, Chingteng↗
NASA NTRS1987

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.↗
NASA NTRS1995

Nonlinear vibrations of thin-walled composite frames

A reduced basis technique and a computational procedure are presented for generating the nonlinear vibrational response, and evaluating the first-order sensitivity coefficients of thin-walled composite frames. The sensitivity coefficients are the derivatives of the nonlinear frequency with respect to the material and lamination parameters of the frame. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The flanges and webs of the frames are modeled by using geometrically nonlinear two-dimensional shell and plate finite elements. The computational procedure can be conveniently divided into three distinct steps. The first step involves the generation of various-order perturbation vectors, and their derivatives with respect to the material and lamination parameters of the frame, using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation vectors as basis vectors, computing the amplitudes of these vectors and the nonlinear frequency of vibration, via a direct variational procedure. The third step consists of using the perturbation vectors, and their derivatives, as basis vectors and computing the sensitivity coefficients of the nonlinear frequency via a second application of the direct variational procedure. Numerical results are presented for semicircular thin-walled frames with I and J section, showing the convergence of the nonlinear frequency and the sensitivity coefficients obtained by both the reduced-basis and perturbation techniques.

Noor, Ahmed K.↗
NASA NTRS2001

Consideration of Moving Tooth Load in Gear Crack Propagation Predictions

Robust gear designs consider not only crack initiation, but crack propagation trajectories for a fail-safe design. In actual gear operation, the magnitude as well as the position of the force changes as the gear rotates through the mesh. A study to determine the effect of moving gear tooth load on crack propagation predictions was performed. Two-dimensional analysis of an involute spur gear and three-dimensional analysis of a spiral-bevel pinion gear using the finite element method and boundary element method were studied and compared to experiments. A modified theory for predicting gear crack propagation paths based on the criteria of Erdogan and Sih was investigated. Crack simulation based on calculated stress intensity factors and mixed mode crack angle prediction techniques using a simple static analysis in which the tooth load was located at the highest point of single tooth contact was validated. For three-dimensional analysis, however, the analysis was valid only as long as the crack did not approach the contact region on the tooth.

Lewicki, David G.↗
NASA NTRS1984

A triangular finite element for thin plates and shells

A finite element modeling technique which utilizes a triangular element called TSHEL with 45 degrees of freedom and seven point integration was tested for analysis of thin plate and shell structures. The element formulation is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. The effectiveness of the present modeling technique which features combined use of elements with kinematic modes and those without kinematic modes to eliminate both locking and spurious kinematic modes at the global structural levels are shown.

Lee, S. W.↗
NASA NTRS1985

Modeling of elastically mounted vertical rotor

The evaluation of the dynamic behavior of a rotating system is possible by means of modal parameters (Eigenvalues and Eigenvectors). A mixed analytical and experimetal approach is used to identify the modal parameters of a specially designed test rig. The modal identification is done both for nonrotating as well as rotating systems. These modal parameters are used to validate a developed Finite Element Model.

Dana, S. S.↗
NASA NTRS1990

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.↗
NASA NTRS1988

The origin of pristine KREEP - Effects of mixing between UrKREEP and the magmas parental to the Mg-rich cumulates

The fact that pristine KREEP basalts generally have bulk-rock molar MgO/(MgO+FeO) ratios lower than average for the lunar crust despite their extraordinarily high incompatible element contents is discussed. Pristine KREEP basalts also have primitive Ni contents compared to mare basalts with comparable REE contents. It is suggested that these facts can be explained through mixing between ancient KREEP precursor materials and primitive Mg-rich melts. Finite-difference models of this mixing, followed by anorthosite assimilation and fractional crystallization provide satisfactory fits to the composition of these basalts. It is suggested that this is a further confirmation of the magmasphere hypothesis.

Warren, Paul H.↗
NASA NTRS1993

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.↗
NASA NTRS1993

Small rocket flowfield diagnostic chambers

Instrumented and optically-accessible rocket chambers are being developed to be used for diagnostics of small rocket (less than 440 N thrust level) flowfields. These chambers are being tested to gather local fluid dynamic and thermodynamic flowfield data over a range of test conditions. This flowfield database is being used to better understand mixing and heat transfer phenomena in small rockets, influence the numerical modeling of small rocket flowfields, and characterize small rocket components. The diagnostic chamber designs include: a chamber design for gathering wall temperature profiles to be used as boundary conditions in a finite element heat flux model; a chamber design for gathering inner wall temperature and static pressure profiles; and optically-accessible chamber designs, to be used with a suite of laser-based diagnostics for gathering local species concentration, temperature, density, and velocity profiles. These chambers were run with gaseous hydrogen/gaseous oxygen (GH2/GO2) propellants, while subsequent versions will be run on liquid oxygen/hydrocarbon (LOX/HC) propellants. The purpose, design, and initial test results of these small rocket flowfield diagnostic chambers are summarized.

Morren, Sybil↗
DOE OSTI2021

CFD/FEA Co-Simulation Framework for Analysis of the Thermal Barrier Coating Design and Its Impact on the HD Diesel Engine Performance

Thermal barrier coatings (TBCs) have been investigated both experimentally and through simulation for mixing controlled combustion (MCC) concepts as a method for reducing heat transfer losses and increasing cycle efficiency, but it is still a very active research area. Early studies were inconclusive, with different groups discovering obstacles to realizing the theoretical potential. Nuanced papers have shown that coating material properties, thickness, microstructure, and surface morphology/roughness all can impact the efficacy of the thermal barrier coating and must be accounted for. Adding to the complexities, a strong spatial and temporal heat flux inhomogeneity exists for mixing controlled combustion (diesel) imposed onto the surfaces from the impinging flame jets. In support of the United States Department of Energy SuperTruck II program goal to achieve 55% brake thermal efficiency on a heavy-duty diesel engines, this study sought to develop a deeper insight into the inhomogeneous heat flux from mixing controlled combustion on thermal barrier coatings and to infer concrete guidance for designing coatings. To that end, a co-simulation approach was developed that couples high-fidelity computational fluid dynamics (CFD) modeling of in-cylinder processes and combustion, and finite element analysis (FEA) modeling of the thermal barrier-coated and metal engine components to resolve spatial and temporal thermal boundary conditions. The models interface at the surface of the combustion chamber; FEA modeling predicts the spatially resolved surface temperature profile, while CFD develops insights into the effect of the thermal barrier coating on the combustion process and the boundary conditions on the gas side. The paper demonstrates the capability of the framework to estimate cycle impacts of the temperature swing at the surface, as well as identify critical locations on the piston/thermal barrier coating that exhibit the highest charge temperature and highest heat fluxes. In addition, the FEA results include predictions of thermal stresses, thus enabling insight into factors affecting coating durability. An example of the capability of the framework is provided to illustrate its use for investigating novel coatings and provide deeper insights to guide future coating design.

42 ENGINEERING↗
NASA NTRS1979

Application of the method of integral relations to unsteady fluid flow problems with shocks

A mixed method using both the method of integral relations and the finite difference technique is developed for the solution of unsteady flow problems. The integral relations method is based upon a chosen interpolating function dependent only on the time domain. The resulting local semidiscrete finite element equations obtained are assembled into a global form. The spatial derivatives at the nodes are replaced by finite difference operators and the discretized nonlinear algebraic system is solved by an iterative scheme. Solutions are obtained for the one-dimensional gasdynamics equation, the one-dimensional wave equation and Burger's model of turbulence. Agreement with other numerical and analytical solutions is excellent for the gasdynamics problem and satisfactory for the Burger equation in cases of small viscous effects.

Wadia, A. R.↗
NASA NTRS1977

Finite element analysis of aeroacoustic jet-flap flows

A computational analysis was performed on the steady, turbulent aerodynamic flowfields associated with a jet-blown flap. For regions devoid of flow separation, a parabolic approximation to the governing time-averaged Navier-Stokes equations was applied. Numerical results are presented for the symmetry plane flow of a slot-nozzle planar jet flap geometry, including prediction of flowfield evolution within the secondary mixing region immediately downstream of the trailing edge. Using a two equation turbulence kinetic energy closure model, rapid generation and decay of large spatial gradients in mean and correlated fluctuating velocity components within the immediate wake region were predicted. Modifications to the turbulent flow structure, as induced by porous surface treatment of the flap, were evaluated. The recirculating flow within a representative discrete slot in the surface was evaluated, using the two dimensional, time-averaged Navier-Stokes equations.

Baker, A. J.↗