SEARCH · Search NASA

Results for “fully linear model”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS2007

An Investigation of High-Cycle Fatigue Models for Metallic Structures Exhibiting Snap-Through Response

A study is undertaken to develop a methodology for determining the suitability of various high-cycle fatigue models for metallic structures subjected to combined thermal-acoustic loadings. Two features of this problem differentiate it from the fatigue of structures subject to acoustic loading alone. Potentially large mean stresses associated with the thermally pre- and post-buckled states require models capable of handling those conditions. Snap-through motion between multiple post-buckled equilibrium positions introduces very high alternating stress. The thermal-acoustic time history response of a clamped aluminum beam structure with geometric and material nonlinearities is determined via numerical simulation. A cumulative damage model is employed using a rainflow cycle counting scheme and fatigue estimates are made for 2024-T3 aluminum using various non-zero mean fatigue models, including Walker, Morrow, Morrow with true fracture strength, and MMPDS. A baseline zero-mean model is additionally considered. It is shown that for this material, the Walker model produces the most conservative fatigue estimates when the stress response has a tensile mean introduced by geometric nonlinearity, but remains in the linear elastic range. However, when the loading level is sufficiently high to produce plasticity, the response becomes more fully reversed and the baseline, Morrow, and Morrow with true fracture strength models produce the most conservative fatigue estimates.

Przekop, Adam↗

DOE OSTI2022

A hybrid meshfree discretization to improve the numerical performance of peridynamic models

Efficient and accurate calculation of spatial integrals is of major interest in the numerical implementation of peridynamics (PD). The standard way to perform this calculation is a particle-based approach that discretizes the strong form of the PD governing equation. This approach has rapidly been adopted by the PD community since it offers some advantages. Additionally, it is computationally cheaper than other available schemes, can conveniently handle material separation, and effectively deals with nonlinear PD models. Nevertheless, PD models are still computationally very expensive compared with those based on the classical continuum mechanics theory, particularly for large-scale problems in three dimensions. This results from the nonlocal nature of the PD theory which leads to interactions of each node of a discretized body with multiple surrounding nodes. Here, we propose a new approach to significantly boost the numerical efficiency of PD models. We propose a discretization scheme that employs a simple collocation procedure and is truly meshfree; i.e., it does not depend on any background integration cells. In contrast to the standard scheme, the proposed scheme requires a much smaller set of neighboring nodes (keeping the same physical length scale) to achieve a specific accuracy and is thus computationally more efficient. Our new scheme is applicable to the case of linear PD models and within neighborhoods where the solution can be approximated by smooth basis functions. Therefore, to fully exploit the advantages of both the standard and the proposed schemes, a hybrid discretization is presented that combines both approaches within an adaptive framework. The high performance of the developed framework is illustrated by several numerical examples, including brittle fracture and corrosion problems in two and three dimensions.

42 ENGINEERING↗

DOE OSTI2016

Thermal Safety Analyses for the Production of Plutonium-238 at the High Flux Isotope Reactor

There has been a considerable effort over the previous few years to demonstrate and optimize the production of plutonium-238 ( 238 Pu) at the High Flux Isotope Reactor (HFIR). This effort has involved resources from multiple divisions and facilities at the Oak Ridge National Laboratory (ORNL) to demonstrate the fabrication, irradiation, and chemical processing of targets containing neptunium-237 ( 237 Np) dioxide (NpO 2 )/aluminum (Al) cermet pellets. A critical preliminary step to irradiation at the HFIR is to demonstrate the safety of the target under irradiation via documented experiment safety analyses. The steady-state thermal safety analyses of the target are simulated in a finite element model with the COMSOL Multiphysics code that determines, among other crucial parameters, the limiting maximum temperature in the target. Safety analysis efforts for this model discussed in the present report include: (1) initial modeling of single and reduced-length pellet capsules in order to generate an experimental knowledge base that incorporate initial non-linear contact heat transfer and fission gas equations, (2) modeling efforts for prototypical designs of partially loaded and fully loaded targets using limited available knowledge of fabrication and irradiation characteristics, and (3) the most recent and comprehensive modeling effort of a fully coupled thermo-mechanical approach over the entire fully loaded target domain incorporating burn-up dependent irradiation behavior and measured target and pellet properties, hereafter referred to as the production model. These models are used to conservatively determine several important steady-state parameters including target stresses and temperatures, the limiting condition of which is the maximum temperature with respect to the melting point. The single pellet model results provide a basis for the safety of the irradiations, followed by parametric analyses in the initial prototypical designs that were necessary due to the limiting fabrication and irradiation data available. The calculated parameters in the final production target model are the most accurate and comprehensive, while still conservative. Over 210 permutations in irradiation time and position were evaluated, and are supported by the most recent inputs and highest fidelity methodology. The results of these analyses show that the models presented in this report provide a robust and reliable basis for previous, current and future experiment safety analyses. In addition, they reveal an evolving knowledge of the steady-state behavior of the NpO 2 /Al pellets under irradiation for a variety of target encapsulations and potential conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1987

A mapped finite difference study of noise propagation in nonuniform ducts with mean flow

The primary objective of this work is to study noise propagation in acoustically lined variable area ducts with mean fluid flow. The method of study is numerical in nature and involves a body-fitted grid mapping procedure in conjunction with a factored-implicit finite difference technique. The mean fluid flow model used is two-dimensional, inviscid, irrotational, incompressible, and nonheat conducting. Fully-coupled solutions of the linearized gasdynamic equations are obtained for both positive and negative Mach numbers as well as for hard and soft wall conditions. The factored-implicit finite difference technique used did give rise to short wavelength perturbations, but these were dampened by the introduction of higher order artificial dissipation terms into the scheme. Results compared favorably with available numerical and experimental data.

Raad, Peter E.↗

NASA NTRS1994

Nonlinear evolution of the coronal magnetic field under reconnective relaxation

Recently, Vekstein et al. (Vekstein, Priest, & Steele 1993) have developed a model for coronal heating in which the corona responds to photospheric footpoint motions by small-scale reconnection events that bring about a relaxed state while conserving magnetic helicity but not field-line connectivity. Vekstein et al. consider a partially open field configuration in which magnetic helicity is ejected to infinity on open field lines but retained in the closed-field region. Under this scheme, they describe the evolution of an initially potential field, in response to helicity injection, in the linear regime. The present work uses numerical calculations to extend the model of Vekstein et al. into the fully nonlinear regime. The results show a rise and bulging of the field lines of the closed-field region with increasing magnetic helicity, to a point where further solutions are impossible. We interpret these solution-sequence endpoints as indicating a possible loss of equilibrium, in the sense that a relaxed equilibrium state may no longer be available to the corona when sufficient helicity has been injected. The rise and bulging behavior is reminiscent of what is observed in a helmet streamer just before the start of a coronal mass ejection (CME), and so our model suggests that a catastrophic loss of magnetic equilibrium might be the initiation mechanism for CMEs. We also find that some choices of boundary conditions can result in qualitative changes in the magnetic topology, with the appearance of magnetic islands. Whether or not this behavior occurs depends on the relative strengths of the fields in the closed- and open-field regions; in particular, island formation is most likely when the open field (which is potential) is strong and thus acts to confine the force-free closed field. Finally, we show that the energy released through reconnective relaxation can be a substantial fraction of the magnetic energy injected into the corona through footpoint motions and may be sufficient for heating the corona above active regions.

Wolfson, R.↗

NASA NTRS1981

Linear polarization of light by two wheat canopies measured at many view angles

The linear polarization and reflection of visible light by wheat as a function of sun-view directions, crop development stage, and wavelength were examined. Two-hundred spectra were taken continuously in wave-lengths from 0.45 to 0.72 Micron in 33 view directions using an Exotech model 20C spectroradiometer six meters above two wheat canopies in the boot and fully headed maturity stages. The analysis results show that the amount of linearly polarized light from the wheat canopies is greatest in the blue spectral region and decreases gradually with increasing wavelength. The results also show that the linearly polarized light from the canopies is generally greatest in the azimuth direction of the Sun and tends toward zero as the view direction tends toward the direction of the hot spot or anti-solar point. It is demonstrated that the single, angle of incidence of sunlight on the leaf, explains almost all of the variation of the amount of polarized light with Sun-view direction.

Bauer, M. E.↗

NASA NTRS1990

Direct use of linear time-domain aerodynamics in aeroservoelastic analysis: Aerodynamic model

The work presented here is the first part of a continuing effort to expand existing capabilities in aeroelasticity by developing the methodology which is necessary to utilize unsteady time-domain aerodynamics directly in aeroservoelastic design and analysis. The ultimate objective is to define a fully integrated state-space model of an aeroelastic vehicle's aerodynamics, structure and controls which may be used to efficiently determine the vehicle's aeroservoelastic stability. Here, the current status of developing a state-space model for linear or near-linear time-domain indicial aerodynamic forces is presented.

Woods, J. A.↗

NASA NTRS1995

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗

DOE OSTI2015

Lattice model of linear telechelic polymer melts. I. Inclusion of chain semiflexibility in the lattice cluster theory

The lattice cluster theory (LCT) for the thermodynamics of polymer systems has recently been reformulated to treat strongly interacting self-assembling polymers composed of fully flexible linear telechelic chains [J. Dudowicz and K. F. Freed, J. Chem. Phys. 136, 064902 (2012)]. Here, we further extend the LCT for linear telechelic polymer melts to include a description of chain semiflexibility, which is treated by introducing a bending energy penalty whenever a pair of consecutive bonds from a single chain lies along orthogonal directions. An analytical expression for the Helmholtz free energy is derived for the model of semiflexible linear telechelic polymer melts. The extension provides a theoretical tool for investigating the influence of chain stiffness on the thermodynamics of self-assembling telechelic polymers, and for further exploring the influence of self-assembly on glass formation in such systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Deep Learning for Multigroup Cross-Section Representation in Two-Step Core Calculations

Here we investigate using deep learning, a type of machine-learning algorithm employing multiple layers of artificial neurons, for the mathematical representation of multigroup cross sections for use in the Griffin reactor multiphysics code for two-step deterministic neutronics calculations. A three-dimensional fuel element typical of a high-temperature gas reactor as well as a two-dimensional sodium-cooled fast reactor lattice are modeled using the Serpent Monte Carlo code, and multigroup macroscopic cross sections are generated for various state parameters to produce a training data set and a separate validation data set. A fully connected, feedforward neural network is trained using the open-source PyTorch machine-learning framework, and its accuracy is compared against the standard piecewise linear interpolation model. Additionally, we provide in this work a generic technique for propagating the cross-section model errors up to the k eff using sensitivity coefficients with the first-order uncertainty propagation rule. Quantifying the eigenvalue error due to the cross-section regression errors is especially practical for appropriately selecting the mathematical representation of the cross sections. We demonstrate that the artificial neural network model produces lower errors and therefore enables better accuracy relative to the piecewise linear model when the cross sections exhibit nonlinear dependencies; especially when a coarse grid is employed, where the errors can be halved by the artificial neural network. However, for linearly dependent multigroup cross sections as found for the sodium-cooled fast reactor case, a simpler linear regression outperforms deeper networks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1992

A fully nonlinear, mixed spectral and finite difference model for thermally driven, rotating flows

Finite difference in time and the meridional plane, in conjunction with a spectral technique in the azimuthal direction, are used to approximate the Navier-Stokes equations in a model that can simulate a variety of thermally driven rotating flows in cylindrical and spherical geometries. Axisymmetric flow, linearized waves relative to a fixed or changing axisymmetric flow, nonlinear waves without wave-wave interaction, and fully nonlinear 3D flow, can in this way be calculated. A reexamination is conducted of the steady baroclinic wave case previously treated by Williams (1971) and Quon (1976).

Miller, Timothy L.↗

DOE OSTI2020

Linear point and sound horizon as purely geometric standard rulers

The baryon acoustic oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the purely geometric BAO methods. To date two different purely geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable (e.g., the BAO). We define the linear point and sound horizon and we fully characterize and compare the two rulers’ cosmological-parameter dependence. We find that they are both geometrical (i.e., independent of the primordial cosmological parameters) within the required accuracy, and that they have the same parameter dependence for a wide range of parameter values. We estimate the rulers’ best-fit values and errors, given the cosmological constraints obtained by the Planck Satellite team from their measurements of the cosmic microwave background temperature and polarization anisotropies. We do this for three different cosmological models encompassed by the purely geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by 0.5σ when we do not fix the spatial geometry to be flat in Λ CDM. Finally, this points toward a sensitivity of the rulers to different cosmological models when they are estimated from the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Valorization of Waste Lipids through Hydrothermal Catalytic Conversion to Liquid Hydrocarbon Fuels with in Situ Hydrogen Production

We demonstrate hydrothermal (300 degrees C, 10 MPa) catalytic conversion of real waste lipids (e.g., waste vegetable oil, sewer trap grease) to liquid hydrocarbon fuels without net need for external chemical inputs (e.g., H2 gas, methanol). A supported bimetallic catalyst (Pt-Re/C; 5 wt % of each metal) previously shown to catalyze both aqueous phase reforming of glycerol (a triacylglyceride lipid hydrolysis coproduct) to H2 gas and conversion of oleic and stearic acid, model unsaturated and saturated fatty acids, to linear alkanes was applied to process real waste lipid feedstocks in water. For reactions conducted with an initially inert headspace gas (N2), waste vegetable oil (WVO) was fully converted into linear hydrocarbons (C15-C17) and other hydrolyzed byproducts within 4.5 h, and H2 gas production was observed. Addition of H2 to the initial reactor headspace accelerated conversion, but net H2 production was still observed, in agreement with results obtained for aqueous mixtures containing model fatty acids and glycerol. Conversion to liquid hydrocarbons with net H2 production was also observed for a range of other waste lipid feedstocks (animal fat residuals, sewer trap grease, dry distiller's grain oil, coffee oil residual). These findings demonstrate potential for valorization of waste lipids through conversion to hydrocarbons that are more compatible with current petroleum-based liquid fuels than the biodiesel and biogas products of conventional waste lipid processing technologies.

09 BIOMASS FUELS↗

DOE OSTI2017

Extended MHD modeling of nonlinear instabilities in fusion and space plasmas

A number of different sub-projects where pursued within this DOE early career project. The primary focus was on using fully nonlinear, curvilinear, extended MHD simulations of instabilities with applications to fusion and space plasmas. In particular, we performed comprehensive studies of the dynamics of the double tearing mode in different regimes and confi gurations, using Cartesian and cyclindrical geometry and investigating both linear and non-linear dynamics. In addition to traditional extended MHD involving Hall term and electron pressure gradient, we also employed a new multi-fluid moment model, which shows great promise to incorporate kinetic effects, in particular off-diagonal elements of the pressure tensor, in a fluid model, which is naturally computationally much cheaper than fully kinetic particle or Vlasov simulations. We used our Vlasov code for detailed studies of how weak collisions effect plasma echos. In addition, we have played an important supporting role working with the PPPL theory group around Will Fox and Amitava Bhattacharjee on providing simulation support for HED plasma experiments performed at high-powered laser facilities like OMEGA-EP in Rochester, NY. This project has support a great number of computational advances in our fluid and kinetic plasma models, and has been crucial to winning multiple INCITE computer time awards that supported our computational modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Subexponentially Growing Hilbert Space and Nonconcentrating Distributions in a Constrained Spin Model

Motivated by recent experiments with two-component Bose–Einstein condensates, we study fully-connected spin models subject to an additional constraint. The constraint is responsible for the Hilbert space dimension to scale only linearly with the system size. We discuss the unconventional statistical physical and thermodynamic properties of such a system, in particular the absence of concentration of the underlying probability distributions. As a consequence, expectation values are less suitable to characterize such systems, and full distribution functions are required instead. Sharp signatures of phase transitions do not occur in such a setting, but transitions from singly peaked to doubly peaked distribution functions of an “order parameter” may be present.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Test-Based Modeling of Photovoltaic Inverter Impact on Distribution Systems

Distributed PV has important effects on overcurrent, overvoltage and islanding protection. These effects are more difficult to analyze because the faulted behavior of PV inverters is less well understood than for rotating machines. Full transient models of the inverter can be accurate, but that approach requires inverter design details and simulation software that doesn’t easily interoperate with utility protection system software. Further, anomalous model behaviors often emerge at different grid conditions, and these are difficult to fix in a proprietary model. Testing provides another approach to characterize inverters, validate models or answer questions of application. This earlier work found major impacts from the inverter design, firmware, operating point and balance-of-system components. Our objective is to characterize PV inverters from black-box testing, leading to standard test reports that could provide parameter sets for widely-implemented models in software tools. The tool vendors can then develop standard inverter models, as they did for rotating machines and exciters. Besides fully detailed transient models, inverters can be represented with steady-state current sources, dynamically controlled current sources or non-linear frequency-dependent black boxes. Type-test data could support the last three kinds of model. Short circuit and sag tests would encompass the basic fault types, reporting the maximum transient and rms currents over all power levels and points-on-wave within four time ranges : 1) 0 to 1 cycles, for close and latch consideration 2) 1 to 4 cycles, for equipment withstand 3) 4 to 30 cycles, for interrupt ratings 4) 30 cycles and beyond, for backup protection analysis Open circuit tests would vary output power and power factor, and report the maximum rms voltage from 1 to 10 cycles after the event. Without models, the type-test data could also provide upper limits on the expected overvoltages and overcurrents.

Wieserman, Laura M.↗

DOE OSTI2023

FEOTS v0.0.0: a new offline code for the fast equilibration of tracers in the ocean

Abstract. In this paper we introduce a new software framework for the offline calculation of tracer transport in the ocean. The Fast Equilibration of Ocean Tracers Software (FEOTS) is an end-to-end set of tools to efficiently calculate tracer distributions on a global or regional sub-domain using transport operators diagnosed from a comprehensive ocean model. To the best of our knowledge, this is the first application of a transport matrix model to an eddying ocean state. While a Newton–Krylov-based equilibration capability is still under development and not presented here, we demonstrate in this paper the transient modeling capabilities of FEOTS in an application focused on the Argentine Basin, where intense eddy activity and the Zapiola Anticyclone lead to strong mixing of water masses. The demonstration illustrates progress in developing offline passive tracer simulation capabilities, while highlighting the challenges of the impulse response functions approach in capturing tracer transports by a non-linear advection scheme. Our future work will focus on improving the computational efficiency of the code to reduce time-to-solution, using different basis functions to better represent non-linear advection operators, applying FEOTS to a parent model with unstructured grids (Ocean Model for Prediction Across Scales, MPAS-Ocean), and fully implementing a Newton–Krylov steady-state solver.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

Solar Spicules: A Review of Recent Models and Targets for Future Observations

Since their discovery over 100 years ago, there have been many suggestions for the origin and development of solar spicules. Because the velocities of spicules are comparable to the sound and Alfven speeds of the low chromosphere, linear theory cannot fully describe them. Consequently, detailed tests of theoretical ideas had to await the development of computing power that only became available during the 1970s. This work reviews theories for spicules and spicule-like features over approximately the past 25 years, with an emphasis on the models based on nonlinear numerical simulations. These models have given us physical insight into wave propagation in the solar atmosphere, and have helped elucidate how such waves, and associated shock waves, may be capable of creating motions and structures on magnetic flux tubes in the lower solar atmosphere. So far, however, it has been difficult to reproduce the most- commonly -quoted parameters for spicules with these models, using what appears to be the most suitable input parameters. A key impediment to developing satisfactory models has been the lack of reliable observational information, which is a consequence of the small angular size and transient lifetime of spicules. I close with a list of key observational questions to be addressed with space-based satellites, such as the currently operating TRACE satellite, and especially the upcoming Solar-B mission. Answers to these questions will help determine which, if any, of the current models correctly explains spicules.

Sterling, Alphonse C.↗