SEARCH · Search NASA

Results for “model simplification”

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 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Variable thickness shear layer aerodynamics revisited

The constant boundary-layer thickness (BLT) figuring in the Ventres (1975) Kernel function can be replaced by a slowly varying BLT, in the case of shear layers of slowly varying thickness. A simplification of the extension by Chi (1976) of Ventres' solution is put forward. The Kernel function in this instance relates pressure on the lifting surface to downwash over the surface. The results should also apply, formally, to three-dimensional compressible unsteady flows, but the accuracy in assuming slowly varying shear BLT remains to be determined. All variants of the shear layer model fail when the shear layer thickness varies rapidly.

Dowell, E. H.↗

NASA NTRS1978

Fracture and crack growth in orthotropic laminates

An approximate solution is developed for the determination of the interlaminar normal and shear stresses in the vicinity of a crack in a three dimensional composite containing unidirectional linearly elastic fibers in an infinite linearly elastic matrix. In order to reduce the complexity of the formulation, certain assumptions are made as to the physically significant stresses to be retained. These simplifications reduce the partial differential equations of elasticity to differential-difference equations which are tractable using Fourier transform techniques. The potential for damaged or debonded zones to be generated by an embedded crack is discussed, and stress concentration factors for fibers near the crack are given. Detailed comparisons are made between the present solution, the analogous two dimensional problem, and corresponding shear-lag models.

Goree, J. G.↗

NASA NTRS1956

Comparison of Theoretical Stresses and Deflections of Multicell Wings with Experimental Results Obtained from Plastic Models

The experimental deflections and stresses of six plastic multicell-wing models of unswept, delta, and swept plan form are presented and compared with previously published theoretical results obtained by the electrical analog method. The comparisons indicate that the theory is reliable except for the evaluation of stresses in the vicinity of the leading edge of delta wings and the leading and trailing edges of swept wings. The stresses in these regions are questionable, apparently because of simplifications employed in idealizing the actual structure for theoretical purposes and because of local effects of concentrated loads.

LOADS AND STRESSES, STRUCTURAL - TORSION↗

NASA NTRS1977

Pure field theories and MACSYMA algorithms

A pure field theory attempts to describe physical phenomena through singularity-free solutions of field equations resulting from an action principle. The physics goes into forming the action principle and interpreting specific results. Algorithms for the intervening mathematical steps are sketched. Vacuum general relativity is a pure field theory, serving as model and providing checks for generalizations. The fields of general relativity are the 10 components of a symmetric Riemannian metric tensor; those of the Einstein-Straus generalization are the 16 components of a nonsymmetric. Algebraic properties are exploited in top level MACSYMA commands toward performing some of the algorithms of that generalization. The light cone for the theory as left by Einstein and Straus is found and simplifications of that theory are discussed.

Ament, W. S.↗

DOE OSTI2020

Parallel algorithms for finding connected components using linear algebra

Finding connected components is one of the most widely used operations on a graph. Optimal serial algorithms for the problem have been known for half a century, and many competing parallel algorithms have been proposed over the last several decades under various different models of parallel computation. This paper presents a class of parallel connected-component algorithms designed using linear-algebraic primitives. These algorithms are based on a PRAM algorithm by Shiloach and Vishkin and can be designed using standard GraphBLAS operations. Here, we demonstrate two algorithms of this class, one named LACC for Linear Algebraic Connected Components, and the other named FastSV which can be regarded as LACC’s simplification. With the support of the highly-scalable Combinatorial BLAS library, LACC and FastSV outperform the previous state-of-the-art algorithm by a factor of up to 12x for small to medium scale graphs. For large graphs with more than 50B edges, LACC and FastSV scale to 4K nodes (262K cores) of a Cray XC40 supercomputer and outperform previous algorithms by a significant margin. This remarkable performance is accomplished by (1) exploiting sparsity that was not present in the original PRAM algorithm formulation, (2) using high-performance primitives of Combinatorial BLAS, and (3) identifying hot spots and optimizing them away by exploiting algorithmic insights.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Limited modified gravity

In this work, we systematically assess several limiting cases of modified gravity, where particular theoretical or observational conditions hold. This framework includes the well known scalar-tensor gravity and No Slip Gravity and No Run Gravity, and we extend it to three new limits: only Run, Only Light, and Only Growth Gravities. These limits give simplifications that allow deeper understanding of modified gravity, including demonstration that gravitational effects on light and matter can have opposite signs in their deviation from general relativity. We also show observational predictions for the different cosmic structure growth rates fσ8 and the ratio of gravitational wave standard siren luminosity to photon standard candle luminosity distance relations, defining a new statistic DG that emphasizes their complementarity and ability to distinguish models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Post-Marcus Theory and Simulation of Interfacial Charge Transfer Dynamics in Organic Semiconducting Materials (Final Report)

The research program developed, validated and applies predictive computational protocols for calculating charge transfer (CT) rates in complex molecular systems, including molecular dyads and triads in liquid solution and solid-state organic semiconducting (OSC) materials. We have established a transformative computational scheme that goes beyond widely used simplifications, to achieve realistic descriptions of CT processes. The approach properly addresses the contribution of molecular environment at ambient conditions to CT processes. Our approach achieves unique insight on investigated CT processes in relevant experimental efforts. The collaborative team included three principal investigators (PIs), with complimentary expertise in classical molecular dynamics simulations and data science (Cheung), state-of-the-art electronic structure calculations (Dunietz), and cutting-edge theory and simulation techniques for modeling energy, charge and coherence transfer dynamics in molecular systems (Geva).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

Quantitative computer representation of propellant processing

With the technology currently available for the manufacture of propellants, it is possible to control the variance of the total specific impulse obtained from the rocket boosters to within approximately five percent. Though at first inspection this may appear to be a reasonable amount of control, when it is considered that any uncertainty in the total kinetic energy delivered to the spacecraft translates into a design with less total usable payload, even this degree of uncertainty becomes unacceptable. There is strong motivation to control the variance in the specific impulse of the shuttle's solid boosters. Any small gains in the predictability and reliability of the booster would lead to a very substantial payoff in earth-to-orbit payload. The purpose of this study is to examine one aspect of the manufacture of solid propellants, namely, the mixing process. The traditional approach of computational fluid mechanics is notoriously complex and time consuming. Certain simplifications are made, yet certain fundamental aspects of the mixing process are investigated as a whole. It is possible to consider a mixing process in a mathematical sense as an operator, F, which maps a domain back upon itself. An operator which demonstrates good mixing should be able to spread any subset of the domain completely and evenly throughout the whole domain by successive applications of the mixing operator, F. Two and three dimensional models are developed and graphical visualization two and three dimensional mixing processes are presented.

Hicks, M. D.↗

NASA NTRS1996

Free-Surface and Contact Line Motion of Liquid in Microgravity

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management in a low-gravity environment. Specific issues include: (1) the mechanic's of liquid masses in large-amplitude motions, (2) the influence of bounding surfaces on the motion, and (3) the ability of such surfaces to control liquid motion by wetting forces, especially when they are augmented by various surface treatments. Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is divided into an 'outer' or inviscid potential-flow region and one or more inner, or viscous dominated, regions. Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior.

Schwartz, Leonard W.↗

NASA NTRS2003

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗

DOE OSTI2021

Dynamical vertex approximation for many-electron systems with spontaneously broken SU(2) symmetry

We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.

2-dimensional systems↗

NASA NTRS1996

Self-Consistent Simulation of the Brownian Stage of Dust Growth

It is a widely accepted view that in proto-planetary accretion disks the collision and following sticking of dust particles embedded in the gas eventually leads to the formation of planetesimals (coagulation). For the smallest dust grains, Brownian motion is assumed to be the dominant source of their relative velocities leading to collisions between these dust grains. As the dust grains grow they eventually couple to the turbulent motion of the gas which then drives the coagulation much more efficiently. Many numerical coagulation simulations have been carried out to calculate the fractal dimension of the aggregates, which determines the duration of the ineffective Brownian stage of growth. Predominantly on-lattice and off-lattice methods were used. However, both methods require simplification of the astrophysical conditions. The aggregates found by those methods had a fractal dimension of approximately 2 which is equivalent to a constant, mass-independent friction time. If this value were valid for the conditions in an accretion disk, this would mean that the coagulation process would finally 'freeze out' and the growth of a planetesimal would be impossible within the lifetime of an accretion disk. In order to investigate whether this fractal dimension is model independent, we simulate self-consistently the Brownian stage of the coagulation by an N-particle code. This method has the advantage that no further assumptions about homogeneity of the dust have to be made. In our model, the dust grains are considered as aggregates built up of spheres. The equation of motion of the dust grains is based on the probability density for the diffusive transport within the gas atmosphere. Because of the very low number density of the dust grains, only 2-body-collisions have to be considered. As the Brownian stage of growth is very inefficient, the system is to be simulated over long periods of time. In order to find close particle pairs of the system which are most likely to undergo a collision, we use a particle-in-cell (PIC) method for the early stages of the simulation where the system is still very homogeneous and a tree method later when the particles are more clustered.

Kempf, S.↗

DOE OSTI2023

Stellarator Simplification using Permanent Magnets (PM4Stell)

This Report describes the design and prototyping of an array of rare-Earth permanent magnets to form a stellarator. This effort was motivated by the hypothesis that the usage of permanent magnets, rather than electromagnetic coils with complex geometry, could reduce the cost of stellarator construction and thereby make increase the feasibility of the stellarator as a technology for a fusion-based power plant. In this project, we have developed novel methods for specifying the positions, shapes, and polarizations of the magnets in the array, and have developed designs for mounting structures and tooling for assembly. We have also performed detailed finite-element modeling to qualify the accuracy of the magnetic field produced by the magnet array as designed, and to confirm that the structure can withstand the forces between the magnets. We have also developed techniques for measuring the magnetic field produced by the array once constructed, as well as for correcting errors in the field arising from misalignments and offsets within the tolerances for mounting and fabrication. Finally, we have constructed a tabletop prototype of a section of the array to qualify the concept for assembling and mounting magnets within the array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Stellarator Simplification using Permanent Magnets (PM4Stell)

47 OTHER INSTRUMENTATION↗

NASA NTRS2015

Space Radiation Transport Code Development: 3DHZETRN

The space radiation transport code, HZETRN, has been used extensively for research, vehicle design optimization, risk analysis, and related applications. One of the simplifying features of the HZETRN transport formalism is the straight-ahead approximation, wherein all particles are assumed to travel along a common axis. This reduces the governing equation to one spatial dimension allowing enormous simplification and highly efficient computational procedures to be implemented. Despite the physical simplifications, the HZETRN code is widely used for space applications and has been found to agree well with fully 3D Monte Carlo simulations in many circumstances. Recent work has focused on the development of 3D transport corrections for neutrons and light ions (Z < 2) for which the straight-ahead approximation is known to be less accurate. Within the development of 3D corrections, well-defined convergence criteria have been considered, allowing approximation errors at each stage in model development to be quantified. The present level of development assumes the neutron cross sections have an isotropic component treated within N explicit angular directions and a forward component represented by the straight-ahead approximation. The N = 1 solution refers to the straight-ahead treatment, while N = 2 represents the bi-directional model in current use for engineering design. The figure below shows neutrons, protons, and alphas for various values of N at locations in an aluminum sphere exposed to a solar particle event (SPE) spectrum. The neutron fluence converges quickly in simple geometry with N > 14 directions. The improved code, 3DHZETRN, transports neutrons, light ions, and heavy ions under space-like boundary conditions through general geometry while maintaining a high degree of computational efficiency. A brief overview of the 3D transport formalism for neutrons and light ions is given, and extensive benchmarking results with the Monte Carlo codes Geant4, FLUKA, and PHITS are provided for a variety of boundary conditions and geometries. Improvements provided by the 3D corrections are made clear in the comparisons. Developments needed to connect 3DHZETRN to vehicle design and optimization studies will be discussed. Future theoretical development will relax the forward plus isotropic interaction assumption to more general angular dependence.

Wilson, John W.↗

NASA NTRS2004

Simulation of Flow Through Breach in Leading Edge at Mach 24

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of holes through the leading edge into a vented cavity. The simulations were generated relatively quickly and early in the investigation by making simplifications to the leading edge cavity geometry. These simplifications in the breach simulations enabled: 1) A very quick grid generation procedure; 2) High fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a 2 inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the 6 inch diameter breach the boundary layer is fully ingested.

Gnoffo, Peter A.↗

DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗

NASA NTRS1990

Analysis of aircraft tires via semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.

Noor, Ahmed K.↗