Search NASA⌕ Search

SEARCH · Search NASA

Results for “nodal formulation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Senstitivty analysis and optimization of nodal point placement for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location - for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, J. I.↗
NASA NTRS1987

Sensitivity derivatives and optimization of nodal point locations for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location; for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, Jocelyn I.↗
NASA NTRS1987

Sensitivity analysis and optimization of nodal point placement for vibration reduction

A method is developed for sensitivity analysis and optimization of nodal point locations in connection with vibration reduction. A straightforward derivation of the expression for the derivative of nodal locations is given, and the role of the derivative in assessing design trends is demonstrated. An optimization process is developed which uses added lumped masses on the structure as design variables to move the node to a preselected location - for example, where low response amplitude is required or to a point which makes the mode shape nearly orthogonal to the force distribution, thereby minimizing the generalized force. The optimization formulation leads to values for added masses that adjust a nodal location while minimizing the total amount of added mass required to do so. As an example, the node of the second mode of a cantilever box beam is relocated to coincide with the centroid of a prescribed force distribution, thereby reducing the generalized force substantially without adding excessive mass. A comparison with an optimization formulation that directly minimizes the generalized force indicates that nodal placement gives essentially a minimum generalized force when the node is appropriately placed.

Pritchard, J. I.↗
DOE OSTI2020

Optimization-based, property-preserving finite element methods for scalar advection equations and their connection to Algebraic Flux Correction

In this paper, we continue our efforts to exploit optimization and control ideas as a common foundation for the development of property-preserving numerical methods. Here we focus on a class of scalar advection equations whose solutions have fixed mass in a given Eulerian region and constant bounds in any Lagrangian volume. Our approach separates discretization of the equations from the preservation of their solution properties by treating the latter as optimization constraints. This relieves the discretization process from having to comply with additional restrictions and makes stability and accuracy the sole considerations in its design. A property-preserving solution is then sought as a state that minimizes the distance to an optimally accurate but not property-preserving target solution computed by the scheme, subject to constraints enforcing discrete proxies of the desired properties. Furthermore, we consider two such formulations in which the optimization variables are given by the nodal solution values and suitably defined nodal fluxes, respectively. A key result of the paper reveals that a standard Algebraic Flux Correction (AFC) scheme is a modified version of the second formulation obtained by shrinking its feasible set to a hypercube. In conclusion, we present numerical studies illustrating the optimization-based formulations and comparing them with AFC

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗
NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗
DOE OSTI2024

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗
NASA NTRS1989

Specialty functions singularity mechanics problems

The focus is in the development of more accurate and efficient advanced methods for solution of singular problems encountered in mechanics. At present, finite element methods in conjunction with special functions, boolean sum and blending interpolations are being considered. In dealing with systems which contain a singularity, special finite elements are being formulated to be used in singular regions. Further, special transition elements are being formulated to couple the special element to the mesh that models the rest of the system, and to be used in conjunction with 1-D, 2-D and 3-D elements within the same mesh. Computational simulation with a least squares fit is being utilized to construct special elements, if there is an unknown singularity in the system. A novel approach is taken in formulation of the elements in that: (1) the material properties are modified to include time, temperature, coordinate and stress dependant behavior within the element; (2) material properties vary at nodal points of the elements; (3) a hidden-symbolic computation scheme is developed and utilized in formulating the elements; and (4) special functions and boolean sum are utilized in order to interpolate the field variables and their derivatives along the boundary of the elements. It may be noted that the proposed methods are also applicable to fluids and coupled problems.

Sarigul, Nesrin↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2007

Flow Separation Side Loads Excitation of Rocket Nozzle FEM

Modern rocket nozzles are designed to operate over a wide range of altitudes, and are also built with large aspect ratios to enable high efficiencies. Nozzles designed to operate over specific regions of a trajectory are being replaced in modern launch vehicles by those that are designed to operate from earth to orbit. This is happening in parallel with modern manufacturing and wall cooling techniques allowing for larger aspect ratio nozzles to be produced. Such nozzles, though operating over a large range of altitudes and ambient pressures, are typically designed for one specific altitude. Above that altitude the nozzle flow is 'underexpanded' and below that altitude, the nozzle flow is 'overexpanded'. In both conditions the nozzle produces less than the maximum possible thrust at that altitude. Usually the nozzle design altitude is well above sea level, leaving the nozzle flow in an overexpanded state for its start up as well as for its ground testing where, if it is a reusable nozzle such as the Space Shuttle Main Engine (SSME), the nozzle will operate for the majority of its life. Overexpansion in a rocket nozzle presents the critical, and sometimes design driving, problem of flow separation induced side loads. To increase their understanding of nozzle side loads, engineers at MSFC began an investigation in 2000 into the phenomenon through a task entitled "Characterization and Accurate Modeling of Rocket Engine Nozzle Side Loads", led by A. Brown. The stated objective of this study was to develop a methodology to accurately predict the character and magnitude of nozzle side loads. The study included further hot-fire testing of the MC-l engine, cold flow testing of subscale nozzles, CFD analyses of both hot-fire and cold flow nozzle testing, and finite element (fe.) analysis of the MC-1 engine and cold flow tested nozzles. A follow on task included an effort to formulate a simplified methodology for modeling a side load during a two nodal diameter fluid/structure interaction for a single moment in time.

Smalley, Kurt B.↗
DOE OSTI2022

Implementation of hybrid finite element method based transport solver in GRIFFIN

A new transport solver option based on the hybrid FEM (HFEM) was implemented in GRIFFIN, the MOOSE-based reactor analysis code, as an effort to support routine core design calculations for advanced reactor applications. The HFEM formulation with P{sub N} (spherical harmonics expansion), akin to the variational nodal method, is effective for solving a spatially homogenized problem with strong transport effect. The residual and Jacobian evaluations of the HFEM weak form were derived and successfully implemented in GRIFFIN, having the diffusion and the PN options available in the new HFEM based transport solver. The performance was tested with the simplified ABTR benchmark problems. The results indicate that the HFEM-based transport solver is a feasible option for solving problems with spatially homogenized and strong streaming by providing superior accuracy with a proper p-refinement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Evaluating Manual Sampling Locations for Regulatory and Emergency Response

Drinking water systems commonly use manual or grab sampling to monitor water quality, identify or confirm issues, and verify that corrective or emergency response actions have been effective. In this paper, the effectiveness of regulatory sampling locations for emergency response is explored. An optimization formulation based on the literature was used to identify manual sampling locations to maximize overall nodal coverage of the system. Results showed that sampling locations could be effective in confirming incidents for which they were not designed. When evaluating sampling locations optimized for emergency response against regulatory scenarios, the average performance was reduced by 3%–4%, while using optimized regulatory sampling locations for emergency response reduced performance by 7%–10%. Secondary constraints were also included in the formulation to ensure geographical and water age diversity with minimal impact on the performance. This work highlighted that regulatory sampling locations provide value in responding to an emergency for these networks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Analytical integration of the tractions induced by non-singular dislocations on an arbitrary shaped triangular quadratic element

Here, an analytical model is proposed to evaluate the nodal force induced by a segment of dislocation upon an arbitrary shaped triangular element. This calculation is required in hybrid methods that associate dislocation dynamics to boundary or finite element to solve simultaneously the evolution of large ensembles of dislocations with complex boundary conditions. Nodal forces are defined as the triple integration of the unbalanced traction field induced by a straight dislocation upon the surface of the element. Following our previous approach (Queyreau et al 2014 Modelling Simul. Mater. Sci. Eng. 22 035004) on a simpler geometry and in the case of linear isotropic elasticity, triple integrals are solved by sequences of integration by parts that exhibit recurrence relations. The traction field is defined and finite everywhere even at the core of dislocations, thanks to the use of the non-singular stress expression formulated by Cai et al (2006 J. Mech. Phys. Solids 54 561–587). The nodal force expressions can be used when considering both a single convolution or double convolution of the Green's function with the core distribution. A solution is also proposed for the case of a semi-infinite segment through the study of the asymptotic behavior of the analytical expressions. The proposed approach is exact and very computationally efficient. The choice of an arbitrary shaped triangular element and quadratic shape functions allow the consideration of complex geometries and comply with automatic meshing procedures. These analytical expressions could also be employed to estimate dislocation interactions with interfaces.

42 ENGINEERING↗
NASA NTRS1984

On the development of hierarchical solution strategies for nonlinear finite element formulations

This paper develops a hierarchical type solution scheme which can handle the field equations associated with nonlinear finite element simulations. The overall procedure possesses various levels of application namely degree of freedom, nodal, elemental, substructural as well as global. In particular iteration, updating, assembly and solution control occurs at the various hierarchical levels. Due to the manner of formulation, the degree of matrix inversion depends on the size of the various hierarchical partitioned groups. In this context, degree of freedom partitioning requires no inversion. To benchmark the overall scheme, the results of several numerical examples are presented.

Padovan, J.↗
NASA NTRS2014

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗
NASA NTRS1995

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗
NASA NTRS1995

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗
NASA NTRS2018

Combining Progressive Nodal Release with the Virtual Crack Closure Technique to Model Fatigue Delamination Growth Without Re-Meshing

The present work summarizes an approach to model mixed-mode 3D fatigue crack growth using the Virtual Crack Closure Technique (VCCT) without requiring re-meshing. It is demonstrated that the proposed approach can be used to simulate crack shapes that do not conform to the underlying mesh. The proposed approach relies solely on Paris Law characterization data to model delamination growth. Growth is determined as a post-processing step at the end of each increment, and hence no convergence issues associated with the progressive nodal release are encountered. This approach can be readily applied using standard solid element formulations and is implemented via an interface user element in Abaqus/Standard.

Carvalho, N. V. De↗