SEARCH · Search NASA

Results for “virtual element methods”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Strain energy release rate analysis of the end-notched flexure specimen using the finite-element method

Two-dimensional finite-element analysis of the end-notched flexure specimen was performed using 8-node isoparametric, parabolic elements to evaluate compliance and mode II strain energy release rates, G sub II. The G sub II values were computed using two different techniques: the virtual crack-closure technique (VCCT) and the rate of change of compliance with crack length (compliance derivative method). The analysis was performed for various crack-length-to-semi-span (a/L) ratios ranging from 0.2 to 0.9. Three material systems representing a wide range of material properties were analyzed. The compliance and strain energy release rates of the specimen calculated with the present finite-element analysis agree very well with beam theory equations including transverse shear. The G sub II values calculated using the compliance derivative method compared extremely well with those calculated using the VCCT. The G sub II values obtained by the compliance derivative method using the top or bottom beam deflections agreed closely with each other. The strain energy release rates from a plane-stress analysis were higher than the plane-strain values by only a small percentage, indicating that either assumption may be used in the analysis. The G sub II values for one material system calculated from the finte-element analysis agreed with one solution in the literature and disagreed with the other solution in the literature.

Salpekar, S. A.↗

NASA NTRSDate not supplied

A Framework for Optimization-Based ISRU Tool Design Using Discrete Element Modeling

Novel robotic excavation technologies are needed to perform in-situ resource utilization (ISRU) tasks at levels required to sustain a long-term presence on the lunar surface. Developing and testing multiple iterations of functional hardware is time and cost prohibitive, thus slowing down the pace of progress and delaying humanity’s settlement of the Moon. High-fidelity, physics-based simulation can reduce the time and effort required to develop and deploy robotic systems [1]. We have adopted this approach to create high-fidelity models of robotic test hardware to enable rapid virtual design and optimization of excavation technologies [2]. Such models can leverage modern computational tools like Discrete Element Method (DEM) simulations that can be coupled with automated design approaches like topology optimization to reduce the amount of prototyping and physical testing needed to realize useful tools.

ISRU↗

NASA NTRS1996

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

NASA NTRS1987

Physically based sensitivity derivatives for finite-element programs

Attention is given to a virtual-work approach to the definition of sensitivity in terms of physical quantities; sensitivity fields are treated as a solution to loadings that consist of body forces, initial strains, and initial stresses. The method can accordingly be implemented in a finite element package, with only minimal knowledge of its programming detail and without access to its source code. The application of the approach to one- and two-dimensional problems is demonstrated. The sensitivity of the structural response for truss structures was computed with respect to the cross-sectional areas of the members.

Barthelemy, Bruno↗

NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗

NASA NTRS2002

Formation metrology and control for large separated optics space telescopes

In this paper we present formation flying performance analysis initial results for a representative large space telescope composed of separated optical elements [Mett 02]. A virtual-structure construct (an equivalent rigid body) is created by unique metrology and control that combines both centralized and decentralized methods. The formation may be in orbit at GEO for super-resolution Earth observation, as in the case of Figure 1, or it may be in an Earth-trailing orbit for astrophysics, Figure 2. Extended applications are envisioned for exo-solar planet interferometric imaging by a formation of very large separated optics telescopes, Figure 3. Space telescopes, with such large apertures and f/10 to f/100 optics, are not feasible if connected by massive metering structures. Instead, the new virtual-structure paradigm of information and control connectivity between the formation elements provides the necessary spatial rigidity and alignment precision for the telescope.

formation flying virtual-structure telescopes↗

NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Finite element crash modeling of composite structures is limited by the inability of current commercial crash codes to accurately model delamination growth. Efforts are made to implement and assess delamination modeling techniques using a current finite element crash code, MSC/DYTRAN. Three methods are evaluated, including a straightforward method based on monitoring forces in elements or constraints representing an interface; a cohesive fracture model proposed in the literature; and the virtual crack closure technique commonly used in fracture mechanics. Results are compared with dynamic double cantilever beam test data from the literature. Examples show that it is possible to accurately model delamination propagation in this case. However, the computational demands required for accurate solution are great and reliable property data may not be available to support general crash modeling efforts. Additional examples are modeled including an impact-loaded beam, damage initiation in laminated crushing specimens, and a scaled aircraft subfloor structures in which composite sandwich structures are used as energy-absorbing elements. These examples illustrate some of the difficulties in modeling delamination as part of a finite element crash analysis.

Fleming, David C.↗

NASA NTRS1999

Delamination Modeling of Composites for Improved Crash Analysis

Fleming, David C.↗

NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗

DOE OSTI2021

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method (DEM) has been intensively used to study the constitutive behaviour of granular materials. However, to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear. This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials. Three groups of models with spherical and clumped particles are investigated from four perspectives: (i) deviatoric stress and volumetric behaviour; (ii) critical state behaviour; (iii) stress-dilatancy relationship; and (iv) the evolution of principal stress ratio against axial strain. The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses, volumetric behaviour, and critical state behaviour observed in experiments for granular materials. On the other hand, some qualitative deviations between experiments and the investigated DEM simulations are also observed, in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain, which are proven to be critical for constitutive modelling. The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective. It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

58 GEOSCIENCES↗

NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗

DOE OSTI2019

Method and apparatus for fast context cloning in a data processing system

A data processing system includes a memory system, a first processing element, a first address translator that maps virtual addresses to system addresses, a second address translator that maps system address to physical addresses, and a task management unit. A first program task uses a first virtual memory space that is mapped to a first system address range using a first table. The context of the first program task includes an address of the first table and is cloned by creating a second table indicative of a mapping from a second virtual address space to a second range of system addresses, where the second range is mapped to the same physical addresses as the first range until a write occurs, at which time memory is allocated and the mapping of the second range is updated. The cloned context includes an address of the second table.

Beard, Jonathan Curtis↗

NASA NTRS1990

Evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

Tan, P. W.↗

DOE OSTI2016

Collective communications apparatus and method for parallel systems

A collective communication apparatus and method for parallel computing systems. For example, one embodiment of an apparatus comprises a plurality of processor elements (PEs); collective interconnect logic to dynamically form a virtual collective interconnect (VCI) between the PEs at runtime without global communication among all of the PEs, the VCI defining a logical topology between the PEs in which each PE is directly communicatively coupled to a only a subset of the remaining PEs; and execution logic to execute collective operations across the PEs, wherein one or more of the PEs receive first results from a first portion of the subset of the remaining PEs, perform a portion of the collective operations, and provide second results to a second portion of the subset of the remaining PEs.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2013

The Virtual Observatory Service TheoSSA: Establishing a Database of Synthetic Stellar Flux Standards I. NLTE Spectral Analysis of the DA-Type White Dwarf G191-B2B ,,,****

Hydrogen-rich, DA-type white dwarfs are particularly suited as primary standard stars for flux calibration. State-of-the-art NLTE models consider opacities of species up to trans-iron elements and provide reliable synthetic stellar-atmosphere spectra to compare with observations. Aims. We will establish a database of theoretical spectra of stellar flux standards that are easily accessible via a web interface. Methods. In the framework of the Virtual Observatory, the German Astrophysical Virtual Observatory developed the registered service TheoSSA. It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code. In case of the DA white dwarf G191−B2B, we demonstrate that the model reproduces not only its overall continuum shape but also the numerous metal lines exhibited in its ultraviolet spectrum. Results. TheoSSA is in operation and contains presently a variety of SEDs for DA-type white dwarfs. It will be extended in the near future and can host SEDs of all primary and secondary flux standards. The spectral analysis of G191−B2B has shown that our hydrostatic models reproduce the observations best at Teff =60 000 +/- 2000K and log g=7.60 +/- 0.05.We newly identified Fe vi, Ni vi, and Zn iv lines. For the first time, we determined the photospheric zinc abundance with a logarithmic mass fraction of −4.89 (7.5 × solar). The abundances of He (upper limit), C, N, O, Al, Si, O, P, S, Fe, Ni, Ge, and Sn were precisely determined. Upper abundance limits of about 10% solar were derived for Ti, Cr, Mn, and Co. Conclusions. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of all astronomical data and cross-calibration between different instruments can be based on the same models and SEDs calculated with different model-atmosphere codes and are easy to compare.

Observatory↗

DOE OSTI2022

What, why and when to go virtual: An international analysis of early adopters of virtual building energy codes inspections

To meet greenhouse gas reduction targets, several countries are pursuing more ambitious policies in their buildings and construction sectors, such as introducing zero net energy/carbon building codes. Countries often report not having enough qualified staff for performing building energy code inspections and many are exploring faster, easier, and more reliable methods to check the compliance of buildings with their codes. Building inspections are a critical element for ensuring code compliance and they have traditionally been performed in person. However, in-person inspections can be labor and travel intensive, costly, and prone to human error. In this paper, the authors explore how virtual inspections, particularly in light of the recent COVID-19 pandemic, have impacted processes for building code compliance checks in jurisdictions and communities around the world. Here, the authors collected data on four key parameters (time and financial savings, scope of inspections, changing practices and technological innovation, and benefits to consumers) from six jurisdictions and communities in five countries (Australia, Canada, Singapore, United Arab Emirates, and the United States) to analyze the impacts of virtual inspections on code compliance checks. The analysis found the greatest value from virtual inspections in geographically dispersed regions and for cities experiencing rapid building construction. The study also explored emerging technologies that are being piloted for virtual inspections. Although many of these technologies hold promise, more resources and capacity are needed to make them viable for use in building energy code inspections.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2018

Simulation of Nonequilibrium Dynamics on a Quantum Computer

We present a hybrid quantum-classical algorithm for the time evolution of out-of-equilibrium thermal states. Here, the method depends on classically computing a sparse approximation to the density matrix and, then, time-evolving each matrix element via the quantum computer. For this exploratory study, we investigate a time-dependent Ising model with five spins on the Rigetti Forest quantum virtual machine and a one spin system on the Rigetti 8Q-Agave quantum processor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1988

A re-evaluation of finite-element models and stress-intensity factors for surface cracks emanating from stress concentrations

A re-evaluation of the 3-D finite-element models and methods used to analyze surface crack at stress concentrations is presented. Previous finite-element models used by Raju and Newman for surface and corner cracks at holes were shown to have ill-shaped elements at the intersection of the hole and crack boundaries. These ill-shaped elements tended to make the model too stiff and, hence, gave lower stress-intensity factors near the hole-crack intersection than models without these elements. Improved models, without these ill-shaped elements, were developed for a surface crack at a circular hole and at a semi-circular edge notch. Stress-intensity factors were calculated by both the nodal-force and virtual-crack-closure methods. Both methods and different models gave essentially the same results. Comparisons made between the previously developed stress-intensity factor equations and the results from the improved models agreed well except for configurations with large notch-radii-to-plate-thickness ratios. Stress-intensity factors for a semi-elliptical surface crack located at the center of a semi-circular edge notch in a plate subjected to remote tensile loadings were calculated using the improved models. The ratio of crack depth to crack length ranged form 0.4 to 2; the ratio of crack depth to plate thickness ranged from 0.2 to 0.8; and the ratio of notch radius to the plate thickness ranged from 1 to 3. The models had about 15,000 degrees-of-freedom. Stress-intensity factors were calculated by using the nodal-force method.

Tan, P. W.↗