Search NASA⌕ Search

SEARCH · Search NASA

Results for “embedded 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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

New Developments in the Embedded Statistical Coupling Method: Atomistic/Continuum Crack Propagation

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain has been enhanced. The concurrent MD-FEM coupling methodology uses statistical averaging of the deformation of the atomistic MD domain to provide interface displacement boundary conditions to the surrounding continuum FEM region, which, in turn, generates interface reaction forces that are applied as piecewise constant traction boundary conditions to the MD domain. The enhancement is based on the addition of molecular dynamics-based cohesive zone model (CZM) elements near the MD-FEM interface. The CZM elements are a continuum interpretation of the traction-displacement relationships taken from MD simulations using Cohesive Zone Volume Elements (CZVE). The addition of CZM elements to the concurrent MD-FEM analysis provides a consistent set of atomistically-based cohesive properties within the finite element region near the growing crack. Another set of CZVEs are then used to extract revised CZM relationships from the enhanced embedded statistical coupling method (ESCM) simulation of an edge crack under uniaxial loading.

Saether, E.↗
NASA NTRS1999

Structural Embeddings: Mechanization with Method

The most powerful tools for analysis of formal specifications are general-purpose theorem provers and model checkers, but these tools provide scant methodological support. Conversely, those approaches that do provide a well-developed method generally have less powerful automation. It is natural, therefore, to try to combine the better-developed methods with the more powerful general-purpose tools. An obstacle is that the methods and the tools often employ very different logics. We argue that methods are separable from their logics and are largely concerned with the structure and organization of specifications. We, propose a technique called structural embedding that allows the structural elements of a method to be supported by a general-purpose tool, while substituting the logic of the tool for that of the method. We have found this technique quite effective and we provide some examples of its application. We also suggest how general-purpose systems could be restructured to support this activity better.

Munoz, Cesar↗
NASA NTRS1987

Method for hygromechanical characterization of graphite/epoxy composite

An experimental method is described for measuring hygroscopic swelling strains and mechanical strains of moisture-conditioned composite specimens. The method consists of embedding encapsulated strain gages in the midplane of the composite laminate; thus it does not interfere with normal moisture diffusion. It is particularly suited for measuring moisture swelling coefficients and for mechanical testing of moisture-conditioned specimens at high strain rates. Results obtained by the embedded gage method were shown to be more reliable and reproducible than those obtained by surface gages, dial gages, or extensometers.

Yaniv, Gershon↗
NASA NTRS2004

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗
DOE OSTI2026

The Extended Embedded Self-Shielding Method in SCALE 6.3/Polaris

The SCALE transport lattice code, Polaris, has been previously developed to generate few-group homogenized cross sections for whole-core nodal diffusion simulators in which the embedded self-shielding method (ESSM) is used for resonance self-shielding calculations to process cross sections. Although the ESSM capability has been very successful in light-water reactor analysis, it may require enhancements in computational efficiency; treatment of spatially dependent resonance self-shielding effects; and handling of interrelated resonance effects among fuel, cladding, and control rod materials. Therefore, this study focuses on improving computational efficiency by using a Dancoff-based Wigner–Seitz approximation combined with a material-based resonance categorization, through which a spatially dependent ESSM capability is developed to accurately estimate self-shielded cross sections inside the fuel. Benchmark results show that the new capability significantly enhances computational efficiency and accuracy for spatially dependent local zones within the fuel and through depletion.

ESSM↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗
NASA NTRS1994

The luminosity functions of embedded stellar clusters. 1: Method of solution and analytic results

We describe a method for computing the history of the luminosity function in a young cluster still forming within a molecular cloud complex. Our method, which utilizes detailed results from stellar evolution theory, assumes that clusters arise from the continuous collapse of dense cloud cores over a protracted period of time. It is also assumed that stars reaching the main sequence are distributed in mass according to a prescribed initial mass function (IMF). We keep track separately of the contributions to the luminosity function from the populations of protostars, pre-main-sequence stars, and main-sequence stars. We derive expressions for the fractional contribution of these populations to both the total number of stars produced and the total cluster luminosity. In our model, the number of protostars rises quickly at first, but then levels off to a nearly constant value, which it maintains until the dispersal of the cloud complex. The number fraction of protostars always decreases with time. Averaged over the life of the parent cloud, this fraction is typically a few percent. The protostar mass distribution can be expressed as an integral over the IMF.

Fletcher, Andre B.↗
NASA NTRS1983

A conservative type-dependent full potential method for the treatment of supersonic flows with embedded subsonic regions

A nonlinear method based on the full potential equation in conservation form, cast in an arbitrary coordinate system, has been developed to treat predominantly supersonic flows with embedded subsonic regions. This type of flow field occurs frequently near the fuselage-canopy junction area and wing leading edge regions for a moderately swept fighter configuration. The method uses the theory of characteristics to accurately monitor the type-dependent flow field. A conservative switching scheme is developed to transition from the supersonic marching algorithm to a subsonic relaxation procedure, and vice versa. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. Results are shown for a few configurations, including a wing-body-wake realistic fighter model having embedded subsonic regions.

Shankar, V.↗
NASA NTRS2000

Immunohistochemical myofiber typing and high-resolution myofibrillar lesion detection in LR white embedded muscle

We have developed a method of fixing, embedding, sectioning, and staining that allows high-resolution detection of myofibrillar structure and myosin immunocytochemical muscle fiber typing in serial semithin sections of LR White plastic embedded muscle at the light microscopic level. Traditional approaches, such as cryostat sections, permit fiber typing, but small myofibrillar lesions (1-3 sarcomeres) are difficult to detect because of section thickness. Semithin sections of hydrophobic resins do not stain well either histochemically or immunocytochemically. Electron microscopy can resolve lesions and discriminate fiber types based on morphology, but the sampling area is small. Our goal was to develop a rapid method for defining both fiber type and high-resolution primary myofibrillar lesion damage. Mild fixation (1-4% paraformaldehyde, 0. 05-0.1% glutaraldehyde) and embedment in a hydrophilic resin (LR White) were used. Myofibrillar structure was extremely well preserved at the light microscopic (LM) level, and lesions could be readily resolved in Toluidine blue stained 500-nm sections. Fiber type was defined by LM immunomyosin staining of serial plastic semithin sections, which demonstrated reciprocal staining patterns for "fast (Sigma M4276) and "total" (skeletal muscle) myosins (Sigma M7523). Copyright 2000 Wiley-Liss, Inc.

NASA Program Fundamental Space Biology↗
NASA NTRS1994

Characterization on defect solids

The purpose of this research has been to develop new semi-empirical techniques to describe ceramics, alloys, and metal/ceramic interfaces with applications in mind that support the materials aspects of the high speed civil transport program (HSCT). HSCT requires methods that aid in the design of alloys and ceramics for new high strength, high temperature applications. Current theoretical methods are not capable of carrying out this mission. Hence, new accurate and more efficient theoretical techniques are needed to facilitate the design of new materials and the examination of their properties. This program concentrated on modeling ceramics, but also dealt with alloy properties to a lesser degree. The primary accomplishment of this research was the development of a new equation of state (EOS) that models the energy/bond length relationship and includes terms that represent charge transfer between cations and anions. The new EOS has been used by researchers at the Naval Research Laboratory to develop a new semi-empirical method for calculating properties of ceramics and metal ceramic interfaces based on the embedded atom method. The results of the director's discretionary fund (DDF) research describing the binding energy relation was the enabling theoretical basis. Also the workers at Cleveland State and NASA Lewis have detailed an approach based on equivalent crystal theory as part of the DDF which is currently under development. In addition, initial results for using the Harris functional method have been developed. Finally, the techniques developed by us have proved useful for the study of high pressure properties of solids and have been used extensively by researchers in this field.

Schlosser, Herbert↗
NASA NTRS2004

Adjoint Algorithm for CAD-Based Shape Optimization Using a Cartesian Method

Adjoint solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape optimization. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (geometric parameters that control the shape). More recently, emerging adjoint applications focus on the analysis problem, where the adjoint solution is used to drive mesh adaptation, as well as to provide estimates of functional error bounds and corrections. The attractive feature of this approach is that the mesh-adaptation procedure targets a specific functional, thereby localizing the mesh refinement and reducing computational cost. Our focus is on the development of adjoint-based optimization techniques for a Cartesian method with embedded boundaries.12 In contrast t o implementations on structured and unstructured grids, Cartesian methods decouple the surface discretization from the volume mesh. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin et developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation. Central to automated shape optimization algorithms is the issue of geometry modeling and control. The need to optimize complex, "real-life" geometry provides a strong incentive for the use of parametric-CAD systems within the optimization procedure. In previous work, we presented an effective optimization framework that incorporates a direct-CAD interface. In this work, we enhance the capabilities of this framework with efficient gradient computations using the discrete adjoint method. We present details of the adjoint numerical implementation, which reuses the domain decomposition, multigrid, and time-marching schemes of the flow solver. Furthermore, we explain and demonstrate the use of CAD in conjunction with the Cartesian adjoint approach. The final paper will contain a number of complex geometry, industrially relevant examples with many design variables to demonstrate the effectiveness of the adjoint method on Cartesian meshes.

Nemec, Marian↗
NASA NTRS1998

A Domain-Decomposed Multi-Level Method for Adaptively Refined Cartesian Grids with Embedded Boundaries

The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.

Aftosmis, M. J.↗
DOE OSTI2025

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1973

Analysis of embedded shock waves calculated by relaxation methods.

The requirements for uniqueness of the calculated jump conditions across embedded shock waves are investigated for type-dependent difference systems used in transonic flow studies. A mathematical analysis shows that sufficient conditions are (1) the equations should be differenced in conservative form and (2) a special difference operator should be used when switching from a hyperbolic to an elliptic operator. The latter results in a consistency condition on the integral equations, rather than the differential, at these points. Calculated jump conditions for several embedded and detached shock waves are analyzed in the physical and hodograph planes. Comparisons are made with previous results, a time-dependent calculation, and data.

Murman, E. M.↗
NASA NTRS1990

Adaptive grid embedding for the two-dimensional Euler equations

A numerical algorithm is presented for solving the two-dimensional flux-split Euler equations using a multigrid method with adaptive grid embedding. The method uses an unstructured data set along with a system of pointers for communication on the irregularly shaped grid topologies. An explicit two-stage time advancement scheme is implemented. A multigrid algorithm is used to provide grid level communication and to accelerate the convergence of the solution to steady state. Results are presented for an NACA 0012 airfoil in a freestream with Mach numbers of 0.95 and 1.054. Excellent resolution of the shock structures is obtained with the adaptive grid embedding method with significantly fewer grid points than the comparable structured grid.

Warren, Gary P.↗