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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Laminar mixing in a small floating zone

The relationship between the flow and solute fields during steady mass transfer of a dilute component is analyzed for multi-cellular rotating flows in the floating zone process of semiconductor growth. When the recirculating flows are weak in relation to the rate of crystal growth, a closed-form solution clearly shows the link between the convection pattern in the melt and the solute distribution across the surface of the growing solid. In the limit of strong convection, finite element calculations demonstrate the tendency of the composition to become uniform over the majority of the melt. The solute segregation in the product crystal is greatest when the recirculating motion is comparable to the rate of crystal growth, and points to the danger in attempting to grow compositionally uniform materials from a nearly convectionless melt.

Harriott, George M.↗
NASA NTRS1989

Sublaminate analysis of interlaminar fracture in composites. I - Analytical model

A simple analysis method based on a transverse shear-deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed-form expressions for the interlaminar shear stresses ahead of the crack, the total strain-energy release rate, and the strain-energy release-rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite-element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.

Armanios, Erian A.↗
NASA NTRS1987

Mixing, transport and combustion in sprays

The state-of-the-art in analytical methods for sprays and drop/turbulence interactions in combustion processes is presented with reference to dilute sprays, and dilute dispersed flows, that contain well-defined dispersed-phase elements in volume fractions lower than 1 percent. Attention is given to the near-injector dense spray region. The methods tested are those of locally homogeneous flow, deterministic separated flow, and stochastic separated flow; the first of these was ineffective due to its ignoring finite interphase transport rates, and the second failed by ignoring turbulent dispersion. Stochastic separated flow methodology yielded encouraging results through its treatment of both finite interphase transport rates and dispersed-phase/turbulence interactions.

Faeth, G. M.↗
DOE OSTI2022

Stable and unstable supersonic stagnation of an axisymmetric rotating magnetized plasma

The Naval Research Laboratory ‘Mag Noh problem’, described in this paper, is a self-similar magnetized implosion flow, which contains a fast magnetohydrodynamic (MHD) outward propagating shock of constant velocity. We generalize the classic Noh ( OSTI Tech. Rep. 577058, 1983) problem to include azimuthal and axial magnetic fields as well as rotation. Our family of ideal MHD solutions is five parametric, each solution having its own self-similarity index, gas gamma, magnetization, the ratio of axial to the azimuthal field and rotation. While the classic Noh problem must have a supersonic implosion velocity to create a shock, our solutions have an interesting three-parametric special case with zero initial velocity in which magnetic tension, instead of implosion flow, creates the shock at $t=0+$ . Our self-similar solutions are indeed realized when we solve the initial value MHD problem with the finite volume MHD code Athena. We numerically investigated the stability of these solutions and found both stable and unstable regions in parameter space. Stable solutions can be used to test the accuracy of numerical codes. Unstable solutions have also been widely used to test how codes reproduce linear growth, transition to turbulence and the practically important effects of mixing. Now we offer a family of unstable solutions featuring all three elements relevant to magnetically driven implosions: convergent flow, magnetic field and a shock wave.

Mechanics↗
DOE OSTI2026

JIMWLK on a quantum computer

We propose a method for solving the Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner (JIMWLK) evolution equation on quantum computers. Our approach exploits the reformulation of the JIMWLK equation as a Lindblad master equation governing the rapidity evolution of the hadronic density matrix, as established in prior work. To render the problem tractable for quantum simulation, we introduce several approximations: the two-dimensional transverse plane is reduced to a one-dimensional radial lattice by assuming azimuthal symmetry of the jump operators; the gauge group is restricted to SU(2); and the infinite Wilson lines of the JIMWLK equation are replaced by finite Wilson links along the light-cone direction. The resulting bosonic Hilbert space is truncated using the electric field basis familiar from Hamiltonian lattice gauge theory, with states restricted to angular momenta 𝑗 ≤ 𝑗 max . We derive the matrix elements of the JIMWLK Lindblad jump operators in this basis. As a benchmark, we demonstrate rapid convergence of the fundamental dipole expectation value with 𝑗 max for both pure and mixed Gaussian initial density matrices. For the simplest truncation, 𝑗 max =1/2, we implement the Lindblad evolution using a quantum simulation algorithm verified with the Qiskit statevector simulator by decomposing the non-unitary evolution operator into a linear combination of unitaries. This work establishes a concrete pathway toward quantum simulation of high-energy quantum chromodynamics evolution equations, with direct relevance to the physics program of the Electron-Ion Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1986

Simple formulas for strain-energy release rates with higher order and singular finite elements

A general finite element procedure for obtaining strain-energy release rates for crack growth in isotropic materials is presented. The procedure is applicable to two-dimensional finite element analyses and uses the virtual crack-closure method. The procedure was applied to non-singular 4-noded (linear), 8-noded (parabolic), and 12-noded (cubic) elements and to quarter-point and cubic singularity elements. Simple formulas for strain-energy release rates were obtained with this procedure for both non-singular and singularity elements. The formulas were evaluated by applying them to two mode I and two mixed mode problems. Comparisons with results from the literature for these problems showed that the formulas give accurate strain-energy release rates.

Raju, I. S.↗
NASA NTRS1988

Convection and mass-transport in laser-induced chemical vapor deposition

Gas flow and energy and species transport in laser-induced chemical vapor deposition (LICVD) of amorphous silicon films by silane pyrolysis are analyzed by finite element analysis of a two-dimensional model for the process. Spatial nonuniformity of the deposited film is shown to result from diffusion controlled transport of products between the beam and substrate. Deposition profiles are affected by buoyancy-driven convection only at increased gas pressures. Horizontal orientation of the reactor with respect to gravity is optimal because the stagnation-like flow, that results adjacent to the substrate, enhances mixing, and smoothes the film profile.

Patnaik, S.↗
NASA NTRS1966

Wind-tunnel boundary interference for V/STOL TESTING

The wake skew angle used in applying the theory of NASA TR R-124 to data correction should be such that the angular deflection of the wake vorticity from the horizontal is one-half that calculated from momentum theory at the lifting element. This usage is in contrast to that of the original paper which used the angle of the mass flow. Because of large-scale recirculation effects, there is a finite lower limit to the test speed at which reliable and correctable data can be obtained in closed wind tunnels. Although a zero-correction wind tunnel for V/STOL testing has not yet been achieved, it is shown that the use of suitably mixed wind-tunnel boundaries can alleviate boundary effects on V/STOL data.

V/STOL AIRCRAFT↗
NASA NTRS1985

A first-order Green's function approach to supersonic oscillatory flow: A mixed analytic and numeric treatment

A frequency domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. In this range the effects of the nonlinear terms in the unsteady supersonic compressible velocity potential equation are negligible and therefore these terms will be omitted. The Green's function method is employed in order to convert the potential flow differential equation into an integral one. This integral equation is then discretized, through standard finite element technique, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration (e.g., finite-thickness wing, wing-body-tail) is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. The long range goal is to develop a comprehensive theory for unsteady supersonic potential aerodynamic which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗
NASA NTRS1990

A Poisson equation formulation for pressure calculations in penalty finite element models for viscous incompressible flows

The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.

Sohn, J. L.↗
NASA NTRS1998

Fracture Mechanics Analysis of Stitched Stiffener-Skin Debonding

An analysis based on plate finite elements and the virtual crack closure technique has been implemented to study the effect of stitching on mode I and mode II strain energy release rates for debond configurations. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered, however, the two compliances and failure loads were assumed to be independent. Both a double cantilever beam (mode I) and a mixed mode skin-stiffener debond configuration were studied. In the double cantilever beam configurations, G(sub I) began to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds. In the mixed-mode skin-stiffener configurations, G(sub I) showed a similar behavior as in the double cantilever beam configurations, however, G(sub u), continued to increase with increasing debond length.

Glaessgen, E. H.↗
NASA NTRS2014

A Numerical and Experimental Study of Damage Growth in a Composite Laminate

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

McElroy, Mark↗
NASA NTRS2017

Computationally Efficient Beam Elements for Accurate Stresses in Sandwich Laminates and Laminated Composites with Delaminations

Laminated composites are prone to delamination failure due to the lack of reinforcement through the thickness. Therefore, during the design process the initiation and propagation of delaminations should be accounted for as early as possible. This paper presents computationally efficient nine degree-of-freedom (dof) and eight-dof shear locking-free beam elements using the mixed form of the refined zigzag theory (RZT(sup(m))). The corresponding nine-dof and eight-dof elements use the anisoparametric and constrained anisoparametric interpolation schemes, respectively, to eliminate shear locking in slender beams. The advantage of the present element over previous RZT beam elements is that no post-processing is required to accurately model the transverse shear stress while maintaining the computational efficiency of a low-order beam element. Comparisons with high-fidelity finite element models and three-dimensional elasticity solutions show that the elements can robustly and accurately predict the displacement field, axial stress and transverse shear stress through the thickness of a sandwich beam or a composite laminate with an embedded delamination. In fact, the accuracy and computational efficiency of predicting stresses in laminates with embedded delaminations make the present elements attractive choices for RZT-based delamination initiation and propagation methodologies available in the literature.

R. M. J. Groh↗
NASA NTRS2014

TPSAS-NF1676L-19616-DND

The present study has three goals: (1) perform an experiment where a simple laminate damage process can be characterized in high detail; (2) evaluate the performance of existing commercially available laminate damage simulation tools by modeling the experiment; (3) observe and understand the underlying physics of damage in a composite honeycomb sandwich structure subjected to low-velocity impact. A quasi-static indentation experiment has been devised to provide detailed information about a simple mixed-mode damage growth process. The test specimens consist of an aluminum honeycomb core with a cross-ply laminate facesheet supported on a stiff uniform surface. When the sample is subjected to an indentation load, the honeycomb core provides support to the facesheet resulting in a gradual and stable damage growth process in the skin. This enables real time observation as a matrix crack forms, propagates through a ply, and then causes a delamination. Finite element analyses were conducted in ABAQUS/Explicit(TradeMark) 6.13 that used continuum and cohesive modeling techniques to simulate facesheet damage and a geometric and material nonlinear model to simulate core crushing. The high fidelity of the experimental data allows a detailed investigation and discussion of the accuracy of each numerical modeling approach.

Mark McElroy↗
NASA NTRS1980

ACCESS 3. Approximation concepts code for efficient structural synthesis: User's guide

A user's guide is presented for ACCESS-3, a research oriented program which combines dual methods and a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and dual algorithms of mathematical programming are applied in the design optimization procedure. This program retains all of the ACCESS-2 capabilities and the data preparation formats are fully compatible. Four distinct optimizer options were added: interior point penalty function method (NEWSUMT); second order primal projection method (PRIMAL2); second order Newton-type dual method (DUAL2); and first order gradient projection-type dual method (DUAL1). A pure discrete and mixed continuous-discrete design variable capability, and zero order approximation of the stress constraints are also included.

Fleury, C.↗
NASA NTRS1981

International Conference on Numerical Methods in Fluid Dynamics, 7th, Stanford University, Stanford and Moffett Field, CA, June 23-27, 1980, Proceedings

Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.

Reynolds, W. C.↗
NASA NTRS1988

Effect of vertical magnetic field on convection and segregation in vertical Bridgman crystal growth

A previous finite-element analysis of vertical Bridgman growth for dilute and nondilute alloys is extended to include the effect of a vertically-aligned magnetic field in the limit of zero magnetic Reynolds number. Calculations are presented for growth of a dilute gallium-germanium alloy in a vertically stabilized Bridgman-Stockbarger system and in a furnace with a uniform temperature gradient imposed along the ampoule. Steady cellular convection driven by radial temperature gradients causes good axial and radial mixing in both systems without a magnetic field. A weak magnetic field decreases the intensity of convection and the effectiveness of solute mixing. The radial nonuniformity is greatest for an intermediate field strength. Stronger fields suppress flow recirculation completely, and lead to uniform solute segregation across the crystal and to diffusion-controlled axial segregation.

Kim, Do Hyun↗
DOE OSTI2024

Simultaneous shape and topology optimization of inflatable soft robots

Simultaneous shape and topology optimization is used to design pressure-activated inflatable soft robots. The pressure loaded boundary is meshed conformingly and shape optimized, while the morphology of the robot is topology optimized. The design objective is to exert maximum force on an object, i.e. to produce soft “grippers”. The robot’s motion is modeled using nearly incompressible finite deformation hyperelasticity. To ensure stability of the robot, the buckling load factors obtained via linearized buckling analyses are constrained. The finite element method is used to evaluate the optimization cost and constraint functions and the adjoint method is employed to compute their sensitivities. The numerical examples produce pressure-driven soft robots with varying complexity. We also compare our simultaneous optimization results to those obtained via sequential topology and then shape optimization.

42 ENGINEERING↗