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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

An analytical and experimental investigation of edge delamination in laminates subjected to tension, bending, and torsion

An integrated two-dimensional finite element was developed to calculate interlaminar stresses and strain energy release rates for the study of delamination in composite laminates subjected to uniaxial tension, bending, and torsion loads. Addressed are the formulation, implementation, and verification of the model. Parametric studies were conducted on the effect of Poisson's ratio mismatch between plies and the stacking sequence on interlaminar stress, and on the effect of delamination opening height and delamination length, due to bending, on strain energy release rate for various laminates. A comparison of strain energy release rates in all-graphite and graphite/glass hybrid laminates is included. The preliminary results of laminates subjected to torsion are also included. Fatigue tension tests were conducted on Mode 1 and mixed mode edge-delamination coupons to establish the relationship between fatigue load vs. onset of delamination cycle. The effect on the fatigue delamination onset of different frequencies (1 and 5 Hz) was investigated for glass, graphite,and their hybrid laminates. Although a 20 percent increase in the static onset-of-delamination strength and a 10 percent increase in ultimate strength resulted from hybridizing the all-graphite laminate with a 90 deg glass ply, the fatigue onset is lower in the hybrid laminate than in the all-graphite laminate.

Chan, Wen S.↗
NASA NTRS1986

Elimination sequence optimization for SPAR

SPAR is a large-scale computer program for finite element structural analysis. The program allows user specification of the order in which the joints of a structure are to be eliminated since this order can have significant influence over solution performance, in terms of both storage requirements and computer time. An efficient elimination sequence can improve performance by over 50% for some problems. Obtaining such sequences, however, requires the expertise of an experienced user and can take hours of tedious effort to affect. Thus, an automatic elimination sequence optimizer would enhance productivity by reducing the analysts' problem definition time and by lowering computer costs. Two possible methods for automating the elimination sequence specifications were examined. Several algorithms based on the graph theory representations of sparse matrices were studied with mixed results. Significant improvement in the program performance was achieved, but sequencing by an experienced user still yields substantially better results. The initial results provide encouraging evidence that the potential benefits of such an automatic sequencer would be well worth the effort.

Hogan, Harry A.↗
NASA NTRS1995

Modeling complex chemical effects in turbulent nonpremixed combustion

Virtually all of the energy derived from the consumption of combustibles occurs in systems which utilize turbulent fluid motion. Since combustion is largely related to the mixing of fluids and mixing processes are orders of magnitude more rapid when enhanced by turbulent motion, efficiency criteria dictate that chemically powered devices necessarily involve fluid turbulence. Where combustion occurs concurrently with mixing at an interface between two reactive fluid bodies, this mode of combustion is called nonpremixed combustion. This is distinct from premixed combustion where flame-fronts propagate into a homogeneous mixture of reactants. These two modes are limiting cases in the range of temporal lag between mixing of reactants and the onset of reaction. Nonpremixed combustion occurs where this lag tends to zero, while premixed combustion occurs where this lag tends to infinity. Many combustion processes are hybrids of these two extremes with finite non-zero lag times. Turbulent nonpremixed combustion is important from a practical standpoint because it occurs in gas fired boilers, furnaces, waste incinerators, diesel engines, gas turbine combustors, and afterburners etc. To a large extent, past development of these practical systems involved an empirical methodology. Presently, efficiency standards and emission regulations are being further tightened (Correa 1993), and empiricism has had to give way to more fundamental research in order to understand and effectively model practical combustion processes (Pope 1991). A key element in effective modeling of turbulent combustion is making use of a sufficiently detailed chemical kinetic mechanism. The prediction of pollutant emission such as oxides of nitrogen (NO(x)) and sulphur (SO(x)) unburned hydrocarbons, and particulates demands the use of detailed chemical mechanisms. It is essential that practical models for turbulent nonpremixed combustion are capable of handling large numbers of 'stiff' chemical species equations.

Smith, Nigel S. A.↗
DOE OSTI2021

Evolution of the electronic structure in Ta 2 NiSe 5 across the structural transition revealed by resonant inelastic x-ray scattering

We utilized high-energy-resolution resonant inelastic X-ray scattering (RIXS) at both the Ta and Ni L 3 -edges to map out element-specific particle-hole excitations in Ta 2 NiSe 5 across the phase transition. Our results reveal a momentum dependent gap-like feature in the low energy spectrum, which agrees well with the band gap in element-specific joint density of states calculations based on ab initio estimates of the electronic structure in both the low temperature monoclinic and the high temperature orthorhombic structure. Below T c , the RIXS energy-momentum map shows a minimal gap at the Brillouin zone center (~0.16 eV), conforming that Ta 2 NiSe 5 possesses a direct band gap in its low temperature ground state. However, inside the gap, no signature of anticipated collective modes with an energy scale comparable to the gap size can be identified. Upon increasing the temperature to above T c , whereas the gap at the zone center closes, the RIXS map at finite momenta still possesses the gross features of the low temperature map, suggesting a substantial mixing between the Ta and Ni orbits in the conduction and valence bands, which does not change substantially across the phase transition. Our experimental observations and comparison to the theoretical calculations lend further support that the phase transition and the corresponding gap opening in Ta 2 NiSe 5 is largely structural by nature with possible minor contribution from the putative exciton condensate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Studies of implicit and explicit solution techniques in transient thermal analysis of structures

Studies aimed at an increase in the efficiency of calculating transient temperature fields in complex aerospace vehicle structures are reported. The advantages and disadvantages of explicit and implicit algorithms are discussed and a promising set of implicit algorithms with variable time steps, known as GEARIB, is described. Test problems, used for evaluating and comparing various algorithms, are discussed and finite element models of the configurations are described. These problems include a coarse model of the Space Shuttle wing, an insulated frame tst article, a metallic panel for a thermal protection system, and detailed models of sections of the Space Shuttle wing. Results generally indicate a preference for implicit over explicit algorithms for transient structural heat transfer problems when the governing equations are stiff (typical of many practical problems such as insulated metal structures). The effects on algorithm performance of different models of an insulated cylinder are demonstrated. The stiffness of the problem is highly sensitive to modeling details and careful modeling can reduce the stiffness of the equations to the extent that explicit methods may become the best choice. Preliminary applications of a mixed implicit-explicit algorithm and operator splitting techniques for speeding up the solution of the algebraic equations are also described.

Adelman, H. M.↗
NASA NTRS2007

Analysis of Composite Panel-Stiffener Debonding Using a Shell/3D Modeling Technique

Interlaminar fracture mechanics has proven useful for characterizing the onset of delaminations in composites and has been used successfully primarily to investigate onset in fracture toughness specimens and laboratory size coupon type specimens. Future acceptance of the methodology by industry and certification authorities, however, requires the successful demonstration of the methodology on the structural level. For this purpose, a panel was selected that is reinforced with stiffeners. Shear loading causes the panel to buckle, and the resulting out-of-plane deformations initiate skin/stiffener separation at the location of an embedded defect. A small section of the stiffener foot, web and noodle as well as the panel skin in the vicinity of the delamination front were modeled with a local 3D solid model. Across the width of the stiffener foot, the mixedmode strain energy release rates were calculated using the virtual crack closure technique. A failure index was calculated by correlating the results with a mixed-mode failure criterion of the graphite/epoxy material. Computed failure indices were compared to corresponding results where the entire web was modeled with shell elements and only a small section of the stiffener foot and panel were modeled locally with solid elements. Including the stiffener web in the local 3D solid model increased the computed failure index. Further including the noodle and transition radius in the local 3D solid model changed the local distribution across the width. The magnitude of the failure index decreased with increasing transition radius and noodle area. For the transition radii modeled, the material properties used for the noodle area had a negligible effect on the results. The results of this study are intended to be used as a guide for conducting finite element and fracture mechanics analyses of delamination and debonding in complex structures such as integrally stiffened panels.

Krueger, Ronald↗
NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 2

Two families of parametrized mixed variational principles for linear electromagnetodynamics are constructed. The first family is applicable when the current density distribution is known a priori. Its six independent fields are magnetic intensity and flux density, magnetic potential, electric intensity and flux density and electric potential. Through appropriate specialization of parameters the first principle reduces to more conventional principles proposed in the literature. The second family is appropriate when the current density distribution and a conjugate Lagrange multiplier field are adjoined, giving a total of eight independently varied fields. In this case it is shown that a conventional variational principle exists only in the time-independent (static) case. Several static functionals with reduced number of varied fields are presented. The application of one of these principles to construct finite elements with current prediction capabilities is illustrated with a numerical example.

Felippa, Carlos A.↗
NASA NTRSDate not supplied

Edge Based Viscous Method for Node-Centered Formulations

This paper presents a novel, efficient, conservative, edge-based method for evaluation of mean flow viscous fluxes and turbulence-model diffusion terms of the Reynolds-averaged Navier-Stokes equations on tetrahedral grids. The new method is implemented in a practical, node-centered, finite-volume computational fluid dynamics solver. The baseline finite-volume scheme that is equivalent to a second-order accurate finite-element Galerkin approximation of viscous stresses is reformulated. The order of operations to compute the cell-based Green-Gauss gradients is changed to combine the operations by edge, which leads to an equivalent formulation on tetrahedral grids, improves efficiency, and preserves the compact discretization stencil based on the nearest neighbors. The computational results presented in this paper verify the implementation of this edge-based method by comparing its accuracy and iterative convergence with those of the well verified and validated baseline formulation. Efficiency gains for residual and Jacobian evaluations result in significant reduction of time to solution. This novel edge-based formulation on tetrahedra can be seamlessly combined with the baseline formulation on cells of other types for computing solutions on mixed-element grids.

Edge Based↗
NASA NTRS1997

Flowfield-Dependent Mixed Explicit-Implicit (FDMEL) Algorithm for Computational Fluid Dynamics

Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Garcia, S. M.↗
DOE OSTI2021

Thermomechanical Microstructural Predictions of Fracture Nucleation of Zircaloy-4 Alloys With δ and ε Hydride Distributions

Abstract A crystalline dislocation-density formulation that was incorporated with a nonlinear finite-element (FE) method was utilized to understand and to predict the thermomechanical behavior of an hexagonal closest packed (h.c.p.) zircaloy system with hydrides with either face-centered cubic (f.c.c.) or body-centered cubic (b.c.c.) hydrides. This formulation was then used with a recently developed fracture methodology that is adapted for finite inelastic strains and multiphase crystalline systems to understand how different microstructurally based fracture modes nucleate and propagate. The interrelated microstructural characteristics of the different crystalline hydride and matrix phases with the necessary orientation relationships (ORs) have been represented, such that a detailed physical understanding of fracture nucleation and propagation can be predicted for the simultaneous thermomechanical failure modes of hydride populations and the matrix. The effects of volume fraction, morphology, crystalline structure, and orientation and distribution of the hydrides on simultaneous and multiple fracture modes were investigated for radial, circumferential, and mixed distributions. Another key aspect was accounting for temperatures changes due to the effects of thermal conduction and dissipated plastic work and their collective effects on fracture. For hydrided aggregates subjected to high temperatures, thermal softening resulted in higher ductility due to increased dislocation-density activity, which led to higher shear strain accumulation and inhibited crack nucleation and growth. The predictions provide validated insights into why circumferential hydrides are more fracture-resistant than radial hydrides for different volume fractions and thermomechanical loading conditions.

Engineering↗
NASA NTRS2008

Towards Verification of Unstructured-Grid Solvers

New methodology for verification of finite-volume computational methods using unstructured grids is presented. The discretization order properties are studied in computational windows, easily constructed within a collection of grids or a single grid. Tests are performed within each window and address a combination of problem-, solution-, and discretization/grid-related features affecting discretization error convergence. The windows can be adjusted to isolate particular elements of the computational scheme, such as the interior discretization, the boundary discretization, or singularities. Studies can use traditional grid-refinement computations within a fixed window or downscaling, a recently-introduced technique in which computations are made within windows contracting toward a focal point of interest. Grids within the windows are constrained to be consistently refined, allowing a meaningful assessment of asymptotic error convergence on unstructured grids. Demonstrations of the method are shown, including a comparative accuracy assessment of commonly-used schemes on general mixed grids and the identification of local accuracy deterioration at boundary intersections. Recommendations to enable attainment of design-order discretization errors for large-scale computational simulations are given.

Thomas, James L.↗
DOE OSTI2025

Design of a Structure for Assembly and Cooling the Magnet of the Next-Generation 45 GHz ECR Ion Source MARS-D

Here, the current Electron Cyclotron Resonance Ion Sources (ECRISs), constructed with Nb-Ti wires and the conventional racetrack-and-solenoid structure, have achieved operating frequencies up to 28 GHz and utilized about 90% of the critical current of the Nb-Ti wire. A Mixed Axial and Radial field System Demonstrator (MARS-D) is being developed at Lawrence Berkeley National Laboratory (LBNL). This system, which consists of an innovative hexagonal Closed-Loop Coil (CLC) and a set of solenoids, can generate higher magnetic fields (up to 150% ) while requiring only about 50% of the superconducting wire, enabling Nb-Ti wires to be used in the next-generation 45 GHz ECRIS. However, the assembly and cooling of such an efficient and compact magnet are particularly challenging due to the small radial gap between the CLC and solenoids, as well as the tight operating temperature margin. To address these challenges, a structure was developed that combines a three-section radially split solenoid mandrel with a series of shrink-fit reinforcement rings and cooling channels. This paper presents the detailed structure, manufacturing method, assembly procedure, impregnation method, mechanical Finite Element Analysis (FEA) comparison, and thermal FEA comparison.

ECR ion sources↗
NASA NTRS2001

Fatigue Life Methodology for Bonded Composite Skin/Stringer Configurations

A methodology is presented for determining the fatigue life of composite structures based on fatigue characterization data and geometric nonlinear finite element (FE) analyses. To demonstrate the approach, predicted results were compared to fatigue tests performed on specimens which represented a tapered composite flange bonded onto a composite skin. In a first step, tension tests were performed to evaluate the debonding mechanisms between the flange and the skin. In a second step, a 2D FE model was developed to analyze the tests. To predict matrix cracking onset, the relationship between the tension load and the maximum principal stresses transverse to the fiber direction was determined through FE analysis. Transverse tension fatigue life data were used to -enerate an onset fatigue life P-N curve for matrix cracking. The resulting prediction was in good agreement with data from the fatigue tests. In a third step, a fracture mechanics approach based on FE analysis was used to determine the relationship between the tension load and the critical energy release rate. Mixed mode energy release rate fatigue life data were used to create a fatigue life onset G-N curve for delamination. The resulting prediction was in good agreement with data from the fatigue tests. Further, the prediction curve for cumulative life to failure was generated from the previous onset fatigue life curves. The results showed that the methodology offers a significant potential to Predict cumulative fatigue life of composite structures.

Krueger, Ronald↗
DOE OSTI2021

Variable Laser Ignition Pathways in Al/Pt Reactive Multilayers across 10 Decades of Pulse Duration

Pulsed laser irradiation is used to investigate the local initiation of rapid, self-propagating formation reactions in Al/Pt multilayers. The single pulse direct laser ignition of these 1.6 μm thick freestanding foils was characterized over 10 decades of pulse duration (10 ms to 150 fs). Finite element, reactive heat transport modeling of the near-threshold conditions has identified three distinct ignition pathways. For milli- to microsecond pulses, ignition occurs following sufficient absorption of laser energy to enable diffusion of Al and Pt between layers such that the heat released from the corresponding exothermic reaction overcomes conductive losses outside the laser-irradiated zone. When pulse duration is decreased into the nanosecond regime, heat is concentrated near the surface such that the Al locally melts, and a portion of the top-most bilayers react initially. The favorable kinetics and additional heat enable ignition. Further reducing pulse duration to hundreds of femtoseconds leads to a third ignition pathway. In this work, while much of the energy from these pulses is lost to ablation, the remaining heat beneath the crater can be sufficiently concentrated to drive a transverse self-propagating reaction, wherein the heat released from mixing at each interface occurs under kinetic conditions capable of igniting the subsequent layer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Ventless pressure control of two-phase propellant tanks in microgravity

This work studies pressurization and pressure control of a large liquid hydrogen storage tank. A finite element model is developed that couples a lumped thermodynamic formulation for the vapor region with a complete solution of the Navier-Stokes and energy equations for the flow and temperature fields in the liquid. Numerical results show that buoyancy effects are strong, even in microgravity, and can reposition a vapor bubble that is initially at the center of the tank to a region near the tank wall in a relatively short time. Long-term tank pressurization with the vapor bubble at the tank wall shows that after an initial transient lasting about a week, the final rate of pressure increase agrees with a purely thermodynamic analysis of the entire tank. However, the final pressure levels are quite different from thermodynamic predictions. Numerical results also show that there is significant thermal stratification in the liquid due to the effects of natural convection. A subcooled jet is used to provide simultaneous cooling and mixing in order to bring the tank pressure back down to its initial value. Three different jet speeds are examined. Although the lowest jet speed is ineffective at controlling the pressure because of insufficient penetration into the liquid region, the highest jet speed is shown to be quite effective at disrupting thermal stratification and reducing the tank pressure in reasonable time.

Weightlessness↗
NASA NTRS2002

Influence of 2D Finite Element Modeling Assumptions on Debonding Prediction for Composite Skin-stiffener Specimens Subjected to Tension and Bending

The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane strain elements as well as three different generalized plane strain type approaches were performed. The computed deflections, skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with lamination length. For more accurate predictions, however, a three-dimensional analysis is required.

Krueger, Ronald↗
NASA NTRS2004

Comparison of 2D Finite Element Modeling Assumptions with Results From 3D Analysis for Composite Skin-Stiffener Debonding

The influence of two-dimensional finite element modeling assumptions on the debonding prediction for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane-strain elements as well as three different generalized plane strain type approaches were performed. The computed skin and flange strains, transverse tensile stresses and energy release rates were compared to results obtained from three-dimensional simulations. The study showed that for strains and energy release rate computations the generalized plane strain assumptions yielded results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress assumption gave the best agreement. Based on this study it is recommended that results from plane stress and plane strain models be used as upper and lower bounds. The results from generalized plane strain models fall between the results obtained from plane stress and plane strain models. Two-dimensional models may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and mixed mode ratios with delamination length. For more accurate predictions, however, a three-dimensional analysis is required.

Krueger, Ronald↗
DOE OSTI2024

Imaginary phonon modes and phonon-mediated superconductivity in Y 2 C 3

For Y 2 C 3 with a superconducting critical temperature ($T_c$) ~18 K, zone-center imaginary optical phonon modes have been found for the high-symmetry I-43d structure due to C dimer wobbling motion and electronic instability from a flat band near Fermi energy. After lattice distortion to the more stable lowest-symmetry P1 structure, these stabilized low-energy phonon modes with a mixed C and Y character carry a strong electronphonon coupling to give rise to the observed sizable $T_c$. Our work shows that compounds with the calculated dynamical instability should not be simply excluded in the high-throughput search for new phonon-mediated superconductors. Moreover, we have studied the phase stability of the I-43d structure by calculating the enthalpy of different structural motifs of binary compounds containing group-IV elements at the 2:3 composition and also exploring the energy landscapes via ab initio molecular dynamics near and out of the I-43d structure. Here, our results show that the I-43d-type structures with C dimers are preferred in the low to medium pressure range. Because of the wobbling motion of the C dimers, there are many local energy minimums with degenerated energies. Thus, the ensemble average of many I-43d-distorted structures with C dimer wobbling motion at finite temperature still gives an overall I-43d structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗