SEARCH · Search NASA

Results for “nonlinear heat conduction”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2010

A Fully Second Order Implicit/Explicit Time Integration Technique for Hydrodynamics Plus Nonlinear Heat Conduction Problems

We present a second order non-linearly consistent implicit/explicit time integration tech¬nique for solving hydrodynamics coupled with nonlinear heat conduction problems. The idea is to hybridize an implicit and an explicit discretization in a nonlinearly consistent way in order to achieve second order time convergent calculations. In this scope, the hydrody¬namics equations are discretized explicitly making use of the capability of well-understood explicit schemes. On the other hand , the nonlinear heat conduction is solved implicitly. Such methods are often referred to as IMEX methods [2,1,3]. The Jacobian-Free Newton Krylov (JFNK) method (e.g, [10,9]) is applied to the problem in such a way as to render a nonlinearly iterated IMEX method. We solve three test problems in order to validate the nu-merical order of the scheme. For each test, we established second order time convergence. We support these numerical results with a Modi?ed Equation Analysis (MEA) [23,22]. The set of equations studied here constitute a base model for radiation hydrodynamics.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1988

Hybrid transfinite element modeling/analysis of nonlinear heat conduction problems involving phase change

The present paper describes the applicability of hybrid transfinite element modeling/analysis formulations for nonlinear heat conduction problems involving phase change. The methodology is based on application of transform approaches and classical Galerkin schemes with finite element formulations to maintain the modeling versatility and numerical features for computational analysis. In addition, in conjunction with the above, the effects due to latent heat are modeled using enthalpy formulations to enable a physically realistic approximation to be dealt computationally for materials exhibiting phase change within a narrow band of temperatures. Pertinent details of the approach and computational scheme adapted are described in technical detail. Numerical test cases of comparative nature are presented to demonstrate the applicability of the proposed formulations for numerical modeling/analysis of nonlinear heat conduction problems involving phase change.

Tamma, Kumar K.↗

NASA NTRS1983

Pseudo-updated constrained solution algorithm for nonlinear heat conduction

This paper develops efficiency and stability improvements in the incremental successive substitution (ISS) procedure commonly used to generate the solution to nonlinear heat conduction problems. This is achieved by employing the pseudo-update scheme of Broyden, Fletcher, Goldfarb and Shanno in conjunction with the constrained version of the ISS. The resulting algorithm retains the formulational simplicity associated with ISS schemes while incorporating the enhanced convergence properties of slope driven procedures as well as the stability of constrained approaches. To illustrate the enhanced operating characteristics of the new scheme, the results of several benchmark comparisons are presented.

Tovichakchaikul, S.↗

NASA NTRS1983

Efficient linear and nonlinear heat conduction with a quadrilateral element

A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2, and 1 for the normalized stabilization parameter lead to the 5-point, 9-point finite difference, and fully integrated finite element operators, respectively, for rectangular meshes and have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.

Liu, W. K.↗

NASA NTRS1984

Efficient linear and nonlinear heat conduction with a quadrilateral element

A method is presented for performing efficient and stable finite element calculations of heat conduction with quadrilaterals using one-point quadrature. The stability in space is obtained by using a stabilization matrix which is orthogonal to all linear fields and its magnitude is determined by a stabilization parameter. It is shown that the accuracy is almost independent of the value of the stabilization parameter over a wide range of values; in fact, the values 3, 2 and 1 for the normalized stabilization parameter lead to the 5-point finite difference, 9-point finite difference and fully integrated finite element operators, respectively, for rectangular meshes; numerical experiments reported here show that the three have identical rates of convergence in the L2 norm. Eigenvalues of the element matrices, which are needed for stability limits, are also given. Numerical applications are used to show that the method yields accurate solutions with large increases in efficiency, particularly in nonlinear problems.

Liu, W. K.↗

NASA NTRS1985

Solution algorithms for nonlinear transient heat conduction analysis employing element-by-element iterative strategies

The particular problems investigated in the present study arise from nonlinear transient heat conduction. One of two types of nonlinearities considered is related to a material temperature dependence which is frequently needed to accurately model behavior over the range of temperature of engineering interest. The second nonlinearity is introduced by radiation boundary conditions. The finite element equations arising from the solution of nonlinear transient heat conduction problems are formulated. The finite element matrix equations are temporally discretized, and a nonlinear iterative solution algorithm is proposed. Algorithms for solving the linear problem are discussed, taking into account the form of the matrix equations, Gaussian elimination, cost, and iterative techniques. Attention is also given to approximate factorization, implementational aspects, and numerical results.

Winget, J. M.↗

NASA NTRS1982

Self-adaptive closed constrained solution algorithms for nonlinear conduction

Self-adaptive solution algorithms are developed for nonlinear heat conduction problems encountered in analyzing materials for use in high temperature or cryogenic conditions. The nonlinear effects are noted to occur due to convection and radiation effects, as well as temperature-dependent properties of the materials. Incremental successive substitution (ISS) and Newton-Raphson (NR) procedures are treated as extrapolation schemes which have solution projections bounded by a hyperline with an externally applied thermal load vector arising from internal heat generation and boundary conditions. Closed constraints are formulated which improve the efficiency and stability of the procedures by employing closed ellipsoidal surfaces to control the size of successive iterations. Governing equations are defined for nonlinear finite element models, and comparisons are made of results using the the new method and the ISS and NR schemes for epoxy, PVC, and CuGe.

Padovan, J.↗

DOE OSTI2022

Latent Space Dynamics Identification

LaSDI is a data-driven physical simulation software that forms a latent space for a given high-fidelity model and discovers a set of ordinary differential equations for the latent space dynamics. It allows a fast and accurate solution process, which is useful for multi-query decision making applications, such as design optimization and uncertainty quantification. The performance of the LaSDI framework is demonstrated on four different problems, i.e., 1D and 2D Burgers equations, nonlinear heat conduction, and radial advection problems. Both linear and nonlinear compression techniques, such as neural network and proper orthogonal decomposition, are used to form a latent space. A concept of local dynamics identification procedure is introduced to enable a parametric model, which enhances the accuracy level over a given parameter space.

Fries, William↗

NASA NTRS1980

Finite-element analysis of nonlinear conduction problems subject to moving fields

Through the use of a space-time warp, specialized moving finite elements are developed that can be employed to generate a nonlinear heat conduction model for situations involving traveling boundary and heat generation fields superposed on an initial state. To facilitate the solution of the resulting nonlinear finite-element formulation, a multilevel heuristic iterative solution strategy is developed. In order to demonstrate the versatility and accuracy of the moving elements and their associated nonlinear solution strategy, the results of several numerical experiments are presented.

Padovan, J.↗

DOE OSTI2022

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Solution of second order quasi-linear boundary value problems by a wavelet method

A wavelet Galerkin method based on expansions of Coiflet-like scaling function bases is applied to solve second order quasi-linear boundary value problems which represent a class of typical nonlinear differential equations. Two types of typical engineering problems are selected as test examples: one is about nonlinear heat conduction and the other is on bending of elastic beams. Numerical results are obtained by the proposed wavelet method. Through comparing to relevant analytical solutions as well as solutions obtained by other methods, we find that the method shows better efficiency and accuracy than several others, and the rate of convergence can even reach orders of 5.8.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2012

COYOTE

COYOTE is designed for the analysis of steady and transient, multi-dimensional, nonlinear heat conduction problems. The spatial discretization technique is based on the Galerkin finite element method. The code contains a library of linear and quadratic elements and incudes bars and shells. COYOTE Implicit and explicit time integration methods, with adaptive time step control, are available for transient problems and fixed-point methods are used for Conjugate gradient iterative methods are used for the matrix solution. Coupled enclosure radiation with efficient view factor computation is included. Simulation capabilities include material motion, dynamic radiation enclosures, contact boundary conditions, multipoint constraints, material (element) addition, and deletion, bulk nodes, chemically reactive materials, and coupling with other mechanics codes. The code is designed for single or multi-processor computing

Gartling, DavidK.↗

DOE OSTI2012

BISON Framework for the Analysis of Performance of Metallic Fuels

This report documents the completion of the LANL milestone, M3FT-13LA0203123-Complete Final Assessment of Metal Fuel Sensitivity & Validation Data Set. We had made progress toward understanding how to develop a kernel for the zirconium diffusion in the first documentation this summer. The progress was limited but the work had identified several issues with the implementation of the diffusion kernel and the phase diagrams. In this report we have developed potential solutions to those issues. We have worked toward developing a more robust diffusion kernel including better smoothing, meshing, and artificial diffusion flux algorithms. This report summarizes the progress we have made on the development of robust diffusion kernel for BISON. In the current work we used appropriate simple constitutive models and the BISON framework being developed at Idaho National Laboratory (INL) in order to create an initial capability to simulate the behavior of metallic fuel under irradiation. Within the MOOSE/BISON framework a new Zirconium diffusion kernel was implemented that substantially improved the reliability and stability for solving Zirconium migration problems in metallic nuclear fuels. Several benchmark problems were created, including a simple 1-D rod and several 2-D segments of a fuel rod that covered several temperature regimes. The simulations were then compared to results from the code Pedernal, which is a “high-performance solver for coupled, nonlinear heat conduction and multi-component species diffusion". Comparison between the BISON based solver and Pedernal showed exceptional agreement for all benchmark cases. Very small differences appeared at longer time frames (2,000 days) for the hotter fuel simulations; however this is believed to be due, primarily, to the different time-stepping algorithms used in the two codes. Overall the domain of the benchmark covered all components of the U-Pu-Zr phase diagram (alpha, delta, beta and gamma phases) and showed excellent agreement when compared to Pedernal based results.

Nuclear Fuel Cycle & Fuel Materials(11)↗

DOE OSTI2015

Fully-Implicit Orthogonal Reconstructed Discontinuous Petrov-Galerkin Method for Multiphysics Problems

This document is focussed on the development of reconstructed discontinuous Galerkin (RDG P n P m ) method, based on orthogonal basis functions. Orthogonality of basis functions is essential for enabling robust and efficient fully-implicit Newton-Krylov based time discretization. The method is designed for generic partial differential equations, including transient hyperbolic, parabolic and/or elliptic operators, which are attributed to many multiphysics problems. Our driving application of interest is additive manufacturing, which requires simulations of fluid-solid-gas systems, with compressibility effects and phase change (melting/ solidification). Thus, we demonstrate the method capabilities for solving nonlinear heat conduction equation coupled with multi-material Navier-Stokes based compressible fluid (in low Mach number limit) systems. We focus our attention to the method accuracy (in both time and space), as well as solvability of linear algebra involved in linear steps of the Newton-based solver. We demonstrate the accuracy and efficacy of the method, emphasizing the advantages from the orthogonality of basis functions, which makes better conditioning of underlying (approximate) Jacobian matrices, and rapid convergence of Krylov methods. iii

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

Multiply scaled constrained nonlinear equation solvers

To improve the numerical stability of nonlinear equation solvers, a partitioned multiply scaled constraint scheme is developed. This scheme enables hierarchical levels of control for nonlinear equation solvers. To complement the procedure, partitioned convergence checks are established along with self-adaptive partitioning schemes. Overall, such procedures greatly enhance the numerical stability of the original solvers. To demonstrate and motivate the development of the scheme, the problem of nonlinear heat conduction is considered. In this context the main emphasis is given to successive substitution-type schemes. To verify the improved numerical characteristics associated with partitioned multiply scaled solvers, results are presented for several benchmark examples.

Padovan, Joe↗

NASA NTRS1987

Steady and transient least square solvers for thermal problems

This paper develops a hierarchical least square solution algorithm for highly nonlinear heat transfer problems. The methodology's capability is such that both steady and transient implicit formulations can be handled. This includes problems arising from highly nonlinear heat transfer systems modeled by either finite-element or finite-difference schemes. The overall procedure developed enables localized updating, iteration, and convergence checking as well as constraint application. The localized updating can be performed at a variety of hierarchical levels, i.e., degree of freedom, substructural, material-nonlinear groups, and/or boundary groups. The choice of such partitions can be made via energy partitioning or nonlinearity levels as well as by user selection. Overall, this leads to extremely robust computational characteristics. To demonstrate the methodology, problems are drawn from nonlinear heat conduction. These are used to quantify the robust capabilities of the hierarchical least square scheme.

Padovan, Joe↗

NASA NTRS1993

Radio observations of the planets - The importance of laboratory measurements

Laboratory data on the line broadening parameters of H2S gas under Uranian/Neptunian conditions, on the far wings of the H2S and NH3 line profiles, and on the dielectric properties of CH4-, NH3-, H2S-, and NH4SH-ice are needed to constrain elemental abundances and understand the dynamics and cloud physics in the atmospheres of the giant planets. Measurements of the absorption coefficient of gaseous H2SO4 at millimeter wavelengths are needed in order to obtain a better understanding of Venus' atmosphere. To determine wind velocity fields in Venus' and Mars' atmospheres, accurate measurements of the center frequencies of the CO lines are necessary. The absorption and scattering properties of lunar soils and/or terrestrial rock powders at frequencies from approximately 1 to 200 GHz, determined in laboratory experiments, would provide a valuable addition to existing data at 450 MHz, 35 GHz, and far infrared frequencies. These data would be used to analyze the microwave spectra of planetary surfaces. Such studies may be helpful in distinguishing the effects of radiative transfer from those of nonlinear heat conduction and internal heat sources.

De Pater, Imke↗