SEARCH · Engineering Papers

Results for “finite element problems”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗

DOE OSTI2019

A Classical Field Theory Formulation for the Numerical Solution of Time Harmonic Electromagnetic Fields

Finite element representations of Maxwell's equations pose unusual challenges inherent to the variational representation of the “curl-curl” equation for the fields. We present a variational formulation which solves for the four-potential instead, based on classical field theory. Borrowing from quantum electrodynamics, we modify the Lagrangian by adding an implicit gauge-fixing term. This reformulation explicitly accounts for Gauss’ law through the coupling between $\phi$ and $\rho$ , and enables the use of nodal basis functions instead of edge elements for time-harmonic problems. Here, we demonstrate how this formulation, adhering to the deeper underlying symmetries of the four-dimensional covariant field description, provides a highly general, robust numerical framework.

43 PARTICLE ACCELERATORS↗

NASA NTRS1980

Computational strategy for the solution of large strain nonlinear problems using the Wilkins explicit finite-difference approach

The STEALTH code system, which solves large strain, nonlinear continuum mechanics problems, was rigorously structured in both overall design and programming standards. The design is based on the theoretical elements of analysis while the programming standards attempt to establish a parallelism between physical theory, programming structure, and documentation. These features have made it easy to maintain, modify, and transport the codes. It has also guaranteed users a high level of quality control and quality assurance.

Hofmann, R.↗

NASA NTRS1996

Overset Grid Methods for Multidisciplinary Applications in Rotorcraft Problems

A methodology for the coupling of an advanced computational fluid dynamics method based on an overset grid flow-solver and an advanced computational structural dynamics method based on a finite element analysis is presented. Various procedures for the fluid-structure interactions modeling along with their limitations are also discussed. The flight test data for the four-bladed UH-60A Blackhawk helicopter rotor is chosen for the validation of the results. Convergence and accuracy are tested by numerical experiments with a single-bladed rotor. A comparison of airload predictions with flight test data as well as with a rigid blade case is presented. Grid and interpolation related issues for this aeroelastic application are described.

Ahmad, J. U.↗

DOE OSTI2018

PolyMat: an efficient Matlab code for multi-material topology optimization

We introduce a Matlab implementation of topology optimization for compliance minimization on unstructured polygonal finite element meshes that efficiently accommodates many materials and many volume constraints. Leveraging the modular structure of the educational code, PolyTop, we extend it to the multi-material version, PolyMat, with only a few modifications. First, a design variable for each candidate material is defined in each finite element. Next, we couple a Discrete Material Optimization interpolation with the existing penalization and introduce a new parameter such that we can employ continuation and smoothly transition from a convex problem without any penalization to a non-convex problem in which material mixing and intermediate densities are penalized. Mixing that remains due to the density filter operation is eliminated via continuation on the filter radius. To accommodate flexibility in the volume constraint definition, the constraint function is modified to compute multiple volume constraints and the design variable update is modified in accordance with the Zhang-Paulino-Ramos Jr. (ZPR) update scheme, which updates the design variables associated with each constraint independently. The formulation allows for volume constraints controlling any subset of the design variables, i.e., they can be defined globally or locally for any subset of the candidate materials. Borrowing ideas for mesh generation on complex domains from PolyMesher, we conclude which design variables are associated with each local constraint of arbitrary geometry. A number of examples are presented to demonstrate the many material capability, the flexibility of the volume constraint definition, the ease with which we can accommodate passive regions, and how we may use local constraints to break symmetries or achieve graded geometries.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Finite element solution of MHD power-law fluid with slip velocity effect and non-uniform heat source/sink

The laminar steady boundary layer flow of an incompressible MHD Power-law fluid past a continuously moving surface is investigated numerically. The study involves the influence of surface slip and non-uniform heat source/sink on the flow and heat transfer. The governing boundary layer flow equations are transformed into non-dimensional, non-linear coupled ordinary differential equations with the help of suitable similarity transformations. The Galerkin finite element method is implemented to crack the resulting system. The impact of different involved physical parameters is exhibited on the dimensionless velocity profile, temperature distributions and rate of heat transfer in graphical and tabular forms for pseudoplastic and dilatant fluids. The local Nusselt number is found to be the decreasing function of slip parameter, temperature and space-dependent heat sink parameter whereas it increases with increasing values of temperature and space dependent heat source parameter. The problem has important application in attaining the sustainable heat transfer rate for the cooling of fluids, especially in heat exchangers used frequently in chemical industry, in order to increase the trustworthiness of a system, as it removes high heat loads from these systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Resolving Rapid Variation in Energy for Particle Transport

Resolving the rapid variation in energy in neutron and thermal radiation transport is needed for the predictive simulation capability in high-energy density physics applications. Energy variation is difficult to resolve due to rapid variations in cross sections and opacities caused by quantized energy levels in the nuclei and electron clouds. In recent work, we have developed a new technique to simultaneously capture slow and rapid variations in the opacities and the solution using homogenization theory, which is similar to multiband (MB) and to the finite-element with discontiguous support (FEDS) method, but does not require closure information. We demonstrated the accuracy and efficiency of the method for a variety of problems. We are researching how to extend the method to problems with multiple materials and the same material but with different temperatures and densities. In this highlight, we briefly describe homogenization theory and some results.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS1973

Effect of out-of-planeness of membrane quadrilateral finite elements.

Investigation of the effect of using plane quadrilateral membrane elements for modeling nonplanar structures. The effect is assessed by analyzing a simplified finite element model with the aid of the structure network analysis program. The results obtained indicate that the use of planar quadrilateral membrane elements for modeling bending problems can lead to large errors if the four points that define the quadrilateral are not in the same plane.

Haftka, R. T.↗

NASA NTRS1974

Time-fixed rendezvous by impulse factoring with an intermediate timing constraint

A method is presented for factoring a two-impulse orbital transfer into a three- or four-impulse transfer which solves the rendezvous problem and satisfies an intermediate timing constraint. Both the time of rendezvous and the intermediate time of a alinement are formulated as any element of a finite sequence of times. These times are integer multiples of a constant plus an additive constant. The rendezvous condition is an equality constraint, whereas the intermediate alinement is an inequality constraint. The two timing constraints are satisfied by factoring the impulses into collinear parts that vectorially sum to the original impulse and by varying the resultant period differences and the number of revolutions in each orbit. Five different types of solutions arise by considering factoring either or both of the two impulses into two or three parts with a limit for four total impulses. The impulse-factoring technique may be applied to any two-impulse transfer which has distinct orbital periods.

Green, R. N.↗

NASA NTRS1981

Static and dynamic characteristics of large deployable space reflectors

A linear numerical model of the structural characteristics of deployable reflectors was developed. Due to cyclic symmetry of the reflector structure about its axis, only one of many segments is modeled using finite elements. The succeeding segments satisfy continuity of displacement and slope at the interface between consecutive segments. This process leads to (N/2+1) static or dynamic problems of smaller order and bandwidth where N is the number of segments. The solution of each reduced problem leads to motions having a distinct circumferential wave number. The dynamic coupling to the feed support structure is studied adopting modal synthesis.

El-Raheb, M.↗

NASA NTRS1984

Application of Iosipescu specimen geometry to determination of shear strength in unidirectional composites

The Iosipescu (1967) test method, in conjunction with a finite element analysis of anisotropic, unidirectionally reinforced composite specimens, can be used to derive composite shear strength values despite the small size of the samples used. When such a specimen is loaded parallel to the fiber direction, fracture occurs by shear within a narrow region near the center of the specimen; this is consistent with the maximum shear area determined by finite element analysis. The resulting fracture surface, as observed, is characteristic of shear failure, occurring either within the matrix or at the fiber-matrix interface. Problems of compressive and mixed mode failure which accrue to short beam shear specimens are by the present method precluded.

Hurwitz, F. I.↗

NASA NTRS1985

Numerical Techniques in Acoustics

This is the compilation of abstracts of the Numerical Techniques in Acoustics Forum held at the ASME's Winter Annual Meeting. This forum was for informal presentation and information exchange of ongoing acoustic work in finite elements, finite difference, boundary elements and other numerical approaches. As part of this forum, it was intended to allow the participants time to raise questions on unresolved problems and to generate discussions on possible approaches and methods of solution.

Baumeister, K. J.↗

NASA NTRS1987

Evaluation of on-line pulse control for vibration suppression in flexible spacecraft

On-line pulse control for vibration suppression in a flexible spacecraft was evaluated. A continuous beam vs. a truss was modeled. A linear finite element model was used to determine the truss characteristics. Control issues outlined are ED pulse actuator development, pseudo pulse algorithm development, and large nonlinear simulation problems.

Bekey, G. A.↗

NASA NTRS1989

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

NASA NTRS1989

Finite element solvers for incompressible fluid flows and heat transfer

Two different finite-element solvers for incompressible viscous flow, i.e., the mixed interpolation method and the SIMPLE-type iterative method, are compared and tested with some benchmark problems. The advantages of the SIMPLE-type iterative method are the decoupling of the governing equations and the use of equal-order interpolation functions for both velocity and pressure. Even though there is a significant difference between the two methods in terms of the pressure field, similar solutions are obtained for the velocity field.

Sohn, Jeong L.↗

NASA NTRS1990

Static and dynamic aeroelastic characterization of an aerodynamically heated generic hypersonic aircraft configuration

This work-in-progress presentation describes an ongoing research activity at the NASA Langley Research Center to develop analytical methods for the prediction of aerothermoelastic stability of hypersonic aircraft including active control systems. The objectives of this research include application of aerothermal loads to the structural finite element model, determination of the thermal effects on flutter, and assessment of active controls technology applied to overcome any potential adverse aeroelastic stability or response problems due to aerodynamic heating- namely flutter suppression and ride quality improvement. For this study, a generic hypersonic aircraft configuration was selected which incorporates wing flaps, ailerons and all-moveable fins to be used for active control purposes. The active control systems would use onboard sensors in a feedback loop through the aircraft flight control computers to move the surfaces for improved structural dynamic response as the aircraft encounters atmospheric turbulence.

Heeg, Jennifer↗

NASA NTRS1993

Probabilistic nonlinear finite element analysis of composite structures

A probabilistic finite element analysis procedure for laminated composite shells is developed. A total Lagrangian finite element formulation, employing a degenerated three-dimensional laminated composite shell element with the full Green-Lagrange strains and first-order shear deformable kinematics, is used. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed, and results are presented in the form of mean and variance of the structural response. Reliability calculations are made by using the first-order reliability method combined with sensitivity derivatives from the finite element analysis. Both ply-level and micromechanics-level random variables are incorporated, the latter by means of the Aboudi micromechanics model. Two sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. In general, the procedure is quite effective in determining the response statistics and reliability for linear and geometric nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

NASA NTRS1993

Transverse cracking and stiffness reduction in composite laminates

A study of transverse cracking mechanism in composite laminates is presented using a singular hybrid finite element model. The model provides the global structural response as well as the precise local crack-tip stress fields. An elasticity basis for the problem is established by employing Lekhnitskii's complex variable potentials and method of eigenfunction expansion. Stress singularities associated with the transverse crack are obtained by decomposing the deformation into the symmetric and antisymmetric modes and proper boundary conditions. A singular hybrid element is thereby formulated based on the variational principle of a modified hybrid functional to incorporate local crack singularities. Axial stiffness reduction due to transverse cracking is studied. The results are shown to be in very good agreement with the existing experimental data. Comparison with simple shear lag analysis is also given. The effects of stress intensity factors and strain energy density on the increase of crack density are analyzed. The results reveal that the parameters approach definite limits when crack densities are saturated, an evidence of the existence of characteristic damage state.

Yuan, F. G.↗