SEARCH · Search NASA

Results for “MATHEMATICAL SOLUTIONS”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1985

Dynamic Creep Buckling: Analysis of Shell Structures Subjected to Time-dependent Mechanical and Thermal Loading

The objective of the present research is to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. A complete, true ab-initio rate theory of kinematics and kinetics for continuum and curved thin structures, without any restriction on the magnitude of the strains or the deformations, was formulated. The time dependence and large strain behavior are incorporated through the introduction of the time rates of metric and curvature in two coordinate systems: fixed (spatial) and convected (material). The relations between the time derivative and the covariant derivative (gradient) were developed for curved space and motion, so the velocity components supply the connection between the equations of motion and the time rates of change of the metric and curvature tensors.

Simitses, G. J.↗

NASA NTRS1991

Formulation of the nonlinear analysis of shell-like structures, subjected to time-dependent mechanical and thermal loading

The object of the research reported herein was to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermomechanical loads. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) can be anticipated and must be considered in developing the mathematical model. The methodology is demonstrated through different problems of extension, shear, and of planar curved beams. Moreover, importance of the inclusion of large strain is clearly demonstrated, through the chosen applications.

Simitses, George J.↗

DOE OSTI2020

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

A finite element program for postbuckling calculations (PSTBKL)

The object of the research reported herein was to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermochemical loads. This report describes the computer program resulting from the research. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) have been anticipated and are considered in developing the mathematical model. The methodology is demonstrated through different problems of extension, shear, and of planar curved beams. Moreover, importance of the inclusion of large strains is clearly demonstrated, through the chosen applications.

Simitses, G. T.↗

NASA NTRS1985

Analysis of shell type structures subjected to time dependent mechanical and thermal loading

A general mathematical model and solution methodologies for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic or static thermomechanical loads is considered. Among the system responses, which are associated with these load conditions, are thermal buckling, creep buckling and ratchetting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model.

Simitses, G. J.↗

DOE OSTI2023

Semi-Analytical Leakage Solutions for Aquifers (SALSA) v1

SALSA is a semi-analytical modeling tool that can provide assessments of pressure perturbations and brine leakage because of fluid injection and extraction activities in sedimentary basins. Sedimentary basins typically contain vertical sequences of near-horizontal aquifers separated by less-permeable aquitards. In the context of geologic carbon sequestration (GCS), while aquitards are relied upon to contain injected buoyant free-phase CO2 and limit inter-aquifer brine flow, the large number of leaky wells (e.g., unused/unsuccessful exploration wells, or improperly plugged water or oil/gas wells) that are present in sedimentary basins creates concerns for potential leakage and groundwater contamination. The mathematical theory and solution method included in SALSA is presented in our recent manuscript under review. SALSA is the first to account for the coupled leakage through aquitards and leaky wells, with geologic pressurization sources. SALSA based on semi-analytical solutions can be computationally very efficient for problems containing many injection and leaky wells in multilayered sedimentary systems, compared to any numerical simulation software requiring local mesh refinement around each well to obtain accurate results. Therefore, stakeholders (e.g., regulators, operators) for site screening, injection and post-injection pressure behavior and assessing leakage risks can use SALSA as a fast-predictive tool in GCS applications in multilayered aquifer systems.

Cihan, Abdullah↗

DOE OSTI2024

Generating Euler Diagrams Through Combinatorial Optimization

Abstract Can a given set system be drawn as an Euler diagram? We present the first method that correctly decides this question for arbitrary set systems if the Euler diagram is required to represent each set with a single connected region. If the answer is yes, our method constructs an Euler diagram. If the answer is no, our method yields an Euler diagram for a simplified version of the set system, where a minimum number of set elements have been removed. Further, we integrate known wellformedness criteria for Euler diagrams as additional optimization objectives into our method. Our focus lies on the computation of a planar graph that is embedded in the plane to serve as the dual graph of the Euler diagram. Since even a basic version of this problem is known to be NP‐hard, we choose an approach based on integer linear programming (ILP), which allows us to compute optimal solutions with existing mathematical solvers. For this, we draw upon previous research on computing planar supports of hypergraphs and adapt existing ILP building blocks for contiguity‐constrained spatial unit allocation and the maximum planar subgraph problem. To generate Euler diagrams for large set systems, for which the proposed simplification through element removal becomes indispensable, we also present an efficient heuristic. We report on experiments with data from MovieDB and Twitter. Over all examples, including 850 non‐trivial instances, our exact optimization method failed only for one set system to find a solution without removing a set element. However, with the removal of only a few set elements, the Euler diagrams can be substantially improved with respect to our wellformedness criteria.

Computer Science↗

NASA NTRS1969

Diffraction of shock wave by moving thin wing.

Flow field pressure distribution due to plane shock wave impinging by thin wing moving in opposite direction, discussing mathematical formulation, analytic solution and applications

Gunzburger, M.↗

NASA NTRS1973

Steady, subsonic, lifting surface theory for wings with swept, partial span, trailing edge control surfaces

A method for computing the lifting pressure distribution on a wing with partial span, swept control surfaces is presented. This method is valid within the framework of linearized, steady, potential flow theory and consists of using conventional lifting surface theory in conjuction with a flap pressure mode. The cause of a numerical instability that can occur during the quadrature of the flap pressure mode is discussed, and an efficient technique to eliminate the instability is derived. This technique is valid for both the flap pressure mode and regular pressure modes and could be used to improve existing lifting surface methods. Examples of the use of the flap pressure mode and comparisons among this method, other theoretical methods, and experiments are given. Discrepancies with experiment are indicated and candidate causes are presented. It is concluded that the method can lead to an efficient and accurate solution of the mathematical problem when a partial span, trailing edge flap is involved.

Medan, R. T.↗

NASA NTRS1976

A simple solution of sound transmission through an elastic wall to a rectangular enclosure, including wall damping and air viscosity effects

A simple solution to the problem of the acoustical coupling between a rectangular structure, its air content, and an external noise source is presented. This solution is a mathematical expression for the normalized acoustic pressure inside the structure. Numerical results for the sound-pressure response for a specified set of parameters are also presented.

Nahavandi, A. N.↗

NASA NTRS1976

Acoustoelasticity

Internal sound fields are considered. Specifically, the interaction between the (acoustic) sound pressure field and the (elastic) flexible wall of an enclosure is discussed. Such problems frequently arise when the vibrating walls of a transportation vehicle induce a significant internal sound field. Cabin noise in various flight vehicles and the internal sound field in an automobile are representative examples. A mathematical model, simplified solutions, and numerical results and comparisons with representative experimental data are briefly considered. An overall conclusion is that reasonable grounds for optimism exist with respect to available theoretical models and their predictive capability.

Dowell, E. H.↗

NASA NTRS1980

Acoustoelasticity - Review and extension

Research on the interaction between the sound pressure field and the flexible wall of an enclosure is reviewed. The topics discussed include: a mathematical model, simplified solutions for sinusoidal excitation, damping, random excitation, numerical results, including comparisons with representative experimental data, and design analysis methodology. Results of various authors using finite element techniques and modal analysis suggest that effective analytical models are available.

Dowell, E. H.↗

NASA NTRS1984

Thermal shock in a circumferentially cracked hollow cylinder with cladding

An theoretical analysis is presented which demonstrates the effect of cladding on the thermal resistance of a circumferentially cracked hollow cylinder. The cladding is assumed to be bonded to the inner wall of the hollow cylinder. The axisymmetric circumferential crack may be either embedded in the cylinder wall or may be an edge crack which passes through the clad and opens into the inner wall of the hollow cylinder. The problem is formulated mathematically and a solution is found which is in the form of a single integral equation. The integral equation is solved numerically and yields estimates of transient temperature distributions, thermal stresses in the uncracked cylinder, and stress intensity factors as a function of time for various cladding thickness to cylinder wall thickness ratios. It is shown that yielding of the clad under certain conditions can result in a reduction in the magnitude of the stress intensity factor for the crack tip in the elastic base material.

Nied, H. F.↗

NASA NTRS1987

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗

NASA NTRS1988

Non-isothermal elastoviscoplastic analysis of planar curved beams

The development of a general mathematical model and solution methodologies, to examine the behavior of thin structural elements such as beams, rings, and arches, subjected to large nonisothermal elastoviscoplastic deformations is presented. Thus, geometric as well as material type nonlinearities of higher order are present in the analysis. For this purpose a complete true abinito rate theory of kinematics and kinetics for thin bodies, without any restriction on the magnitude of the transformation is presented. A previously formulated elasto-thermo-viscoplastic material constitutive law is employed in the analysis. The methodology is demonstrated through three different straight and curved beams problems.

Simitses, G. J.↗

NASA NTRS1988

Non-isothermal buckling behavior of viscoplastic shell structures

Described are the mathematical model and solution methodologies for analyzing the structural response of thin, metallic elasto-viscoplastic shell structures under large thermomechanical loads and their non-isothermal buckling behavior. Among the system responses associated with these loads and conditions are snap-through, buckling, thermal buckling, and creep buckling. This geometric and material nonlinearities (of high order) can be anticipated and are considered in the model and the numerical treatment.

Riff, Richard↗

NASA NTRS1989

Computational Methods for Structural Mechanics and Dynamics, part 1

The structural analysis methods research has several goals. One goal is to develop analysis methods that are general. This goal of generality leads naturally to finite-element methods, but the research will also include other structural analysis methods. Another goal is that the methods be amenable to error analysis; that is, given a physical problem and a mathematical model of that problem, an analyst would like to know the probable error in predicting a given response quantity. The ultimate objective is to specify the error tolerances and to use automated logic to adjust the mathematical model or solution strategy to obtain that accuracy. A third goal is to develop structural analysis methods that can exploit parallel processing computers. The structural analysis methods research will focus initially on three types of problems: local/global nonlinear stress analysis, nonlinear transient dynamics, and tire modeling.

Stroud, W. Jefferson↗