SEARCH · Search NASA

Results for “Initial conditions problem”

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 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Micromechanical predictions of crack initiation, propagation and crack growth resistance in boron/aluminum composites

An elastoplastic, axisymmetric finite element model has been used to predict the initiation and propagation of a crack in a composite model consisting of a single broken boron fiber embedded in an annular sheath of aluminum matrix. The accuracy of the axisymmetric finite element model for crack problems has been established by solving the classical problem of a penny-shaped crack in a thick cylindrical rod under axial tension. Also, the stress intensity factors predicted by the present numerical model are compared with continuum results. A constant displacement boundary condition applied during an increment of crack growth permits a substantial amount of stable crack growth in the matrix material. The concept of Crack Growth Resistance Curves (KR-curves) has been used to determine the point of crack instability

Mahishi, J. M.↗

NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗

NASA NTRS1993

Transient liquid-crystal technique used to produce high-resolution convective heat-transfer-coefficient maps

In this transient technique the preheated isothermal model wall simulates the classic one-dimensional, semi-infinite wall heat transfer conduction problem. By knowing the temperature of the air flowing through the model, the initial temperature of the model wall, and the surface cooling rate measured at any location with time (using the fast-response liquid-crystal patterns recorded on video tape), the heat transfer coefficient can be calculated for the color isothermal pattern produced. Although the test was run transiently, the heat transfer coefficients are for the steady-state case. The upstream thermal boundary condition was considered to be isothermal. This transient liquid-crystal heat-transfer technique was used in a transient air tunnel in which a square-inlet, 3-to-1 exit transition duct was placed. The duct was preheated prior to allowing room temperature air to be suddenly drawn through it. The resulting isothermal contours on the duct surfaces were revealed using a surface coating of thermochromic liquid crystals that display distinctive colors at particular temperatures. A video record was made of the temperature and time data for all points on the duct surfaces during each test. The duct surfaces were uniformly heated using two heating systems: the first was an automatic temperature-controlled heater blanket completely surrounding the test duct like an oven, and the second was an internal hot-air loop through the inside of the test duct. The hot-air loop path was confined inside the test duct by insulated heat dams located at the inlet and exit ends of the test duct. A recirculating fan moved hot air into the duct inlet, through the duct, out of the duct exit, through the oven, and back to the duct inlet. The temperature nonuniformity of the test duct model wall was held very small. Test results are reported for two inlet Reynolds numbers of 200,000 and 1,150,000 (based on the square-inlet hydraulic diameter) and two free-stream turbulence intensities of about 1 percent, which is typical of wind tunnels, and up to 20 percent (using a grid), which is typical of real engine conditions.

Hippensteele, Steven A.↗

NASA NTRSDate not supplied

The Containment Assurance Risk Framework of the Mars Sample Return Program

The Mars Sample Return campaign aims at bringing rock and atmospheric samples from Mars to Earth through a series of robotic missions. These missions would collect the samples being cached and deposited on Martian soil by the Perseverance rover, place them in a container, and launch them into Martian orbit for subsequent capture by an orbiter that would bring them back. Given there exists a non-zero probability that the samples contain biological material, precautions are being taken to design systems that would break the chain of contact between Mars and Earth. These include techniques such as sterilization of Martian particles, redundant containment vessels, and a robust reentry capsule capable of accurate landings without a parachute. Requirements exist that the probability of containment not assured of Martian-contaminated material into Earth’s biosphere be less than one in a million. To demonstrate compliance with this strict requirement, a statistical framework was developed to assess the likelihood of containment loss during each sample return phase and make a statement about the total combined mission probability of containment not assured. The work presented here describes this framework, which considers failure modes or fault conditions that can initiate failure sequences ultimately leading to containment not assured. Reliability estimates are generated from databases, design heritage, component specifications, or expert opinion in the form of probability density functions or point estimates and provided as inputs to the mathematical models that simulate the different failure sequences. The probabilistic outputs are then combined following the logic of several fault trees to compute the ultimate probability of containment not assured. Given the multidisciplinary nature of the problem and the different types of mathematical models used, the statistical tools needed for analysis are required to be computationally efficient. While standard Monte Carlo approaches are used for fast models, a multi-fidelity approach to rare event probabilities is proposed for expensive models. In this paradigm, inexpensive low-fidelity models are developed for computational acceleration purposes while the expensive high-fidelity model is kept in the loop to retain accuracy in the results. This work presents an example of end-to-end application of this framework highlighting the computational benefits of a multi-fidelity approach.

Giuseppe Cataldo↗

NASA NTRS1993

Dynamics of the atmosphere of Mars: Review of some outstanding problems

Theoretical and modeling studies based on fragmentary data are responsible for our current level of understanding of the dynamics of Mars' atmosphere at the planetary scale. Largely because the boundary conditions and governing parameters are such that the large-scale circulation regime is similar to that or Earth, it has been possible to derive a remarkable amount of information from a relatively small number of observations. Nevertheless, the level of detail available from theoretical and modeling studies far exceeds the validation capability of existing data. Global observations from Mars Observer instruments are expected to go far to remedy this mismatch. This review will present and examine several outstanding dynamical questions, which Mars Observer measurements should help to answer. These include the following: What are the temporal and three-dimensional spatial distributions of thermal tides? What is the distribution of momentum flux divergence due to internal gravity waves? Are normal modes an important component of the general circulation? If so, how do they vary with season? Do normal modes contribute to the initiation of planetary-scale dust storms? What are the details of the global circulation response to large and rapid injections of dust? Current understanding of these problems, the potential role of Mars Observer observations for addressing them, and the implications of answers for our understanding of the present climate and past climate regimes will be discussed.

Leovy, Conway B.↗

NASA NTRS2003

Macrovoid Defect Growth during Evaporative Casting of Polymeric Membranes

Macrovoid (MV) formation is a significant problem in evaporatively cast polymeric membranes. MVs are large, elongated or teardrop-shaped pores (~10-50 micron) that can impair membrane structural integrity. Although MVs have been extensively studied, there is no general agreement on the mechanisms governing MV growth. Recently, our research group has formulated the solutocapillary convection (SC) hypothesis, which contends that MV growth involves three principal forces: a Marangoni force generated by surface tension gradients within the MV interface, a viscous drag force, and a gravitationally induced body force. Two sets of complementary experiments were conducted to test the SC hypothesis. Ground-based videomicroscopy flow-visualization (VMFV) was utilized to measure the flow velocities at the MV-casting solution interface and deep within the casting solution. The measurements were performed with casting solutions containing 10 wt% cellulose acetate (CA), 30 wt% H2O, 60 wt% acetone, and 200- ppm TiO2 particles for flow visualization, and the surface tension was controlled by surfactant addition. Qualitatively, the experiments indicated that MV growth occurs in three distinct phases: (1) a very rapid initial growth period, (2) a much slower growth phase, and (3) absorption of selected MVs into the expanding demixed region. The presence of tracer particles inside the MVs suggests the presence of a convective flow, which transfers the particles from the bulk solution to the MV interior. Although the VMFV experiments did not establish any surfactant effect on the interfacial velocities, a statistically significant effect on the MV number density was observed. In the second set of experiments, membranes were cast aboard a KC-135 aircraft under 0-g and 2-g conditions. Despite careful attention to the design and fabrication of the membrane casting apparatus (MCA), several problems were encountered, the most significant of which was the contamination of the casting solution by the activated carbon particles used for solvent absorption.

Greenberg, A. R.↗

DOE OSTI2026

Automatic Calibration and Health Monitoring of Infrastructure Sensors

Smart transportation infrastructure relies on networks of heterogeneous sensors - cameras, radars, and lidars - continuously monitoring traffic conditions. However, executing the initial spatial calibration of multiple sensors and the subsequent health monitoring presents significant operational challenges. Environmental factors, mechanical vibrations, and gradual drift cause spatial misalignment, degrading fusion performance and tracking accuracy. Traditional calibration approaches require manual intervention with specialized targets or survey equipment, resulting in service interruptions and high maintenance costs. This work presents an automated framework for initial calibration and continuous health monitoring without human intervention or service disruption. Our approach addresses two critical problems: (1) detecting when sensors become miscalibrated during operation, and (2) automatically re-establishing spatial alignment using only operational traffic data. The health monitoring component analyzes measurement innovations - differences between sensor observations and predicted object states - to detect systematic biases indicative of calibration drift. By computing bias magnitude, directional consistency, and rejection rates, the system identifies miscalibrations as small as 0.5 meters. Unlike traditional methods requiring known calibration targets, our diagnostic operates continuously on live traffic observations, enabling early detection before fusion quality degrades. The automatic recalibration algorithm leverages overlapping sensor fields-of-view and temporal correlation of vehicle observations. Using graph-based optimization, the system automatically discovers which sensor pairs observe common regions, estimates pairwise spatial transformations using RANSAC-based robust estimation, and jointly optimizes all sensor poses through bundle adjustment. The framework handles practical deployment challenges, including different sensor sampling rates (1-10 Hz), varying installation positions, unknown orientations, and limited overlap regions (>10%). When approximate sensor positions are available from installation surveys (+/-1m accuracy), the algorithm additionally estimates sensor orientations, refining both position and rotation to sub-meter and sub-degree accuracy. We validate the framework on multi-hour traffic datasets from six heterogeneous sensors with sampling rates ranging from 1 Hz to 10 Hz. Results demonstrate successful calibration even with sparse overlap (<20%) and automatic detection of miscalibrations exceeding 0.8 meters. This work enables a "deploy-and-forget" sensor infrastructure that maintains calibration autonomously, reducing maintenance costs while improving tracking accuracy. The techniques generalize beyond transportation to any multi-sensor monitoring application requiring robust spatial alignment, including smart cities, industrial monitoring, and surveillance systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1973

Nonlinear longitudinal combustion instability in rocket motors.

A new analytical technique for the solution of nonlinear longitudinal combustion instability problems in rocket combustors is developed. Using relatively little computation time, this technique is capable of predicting the transient and limit cycle behavior of the combustion instability oscillations as well as the disturbance amplitude required to trigger an instability in a linearly stable motor. The limit cycle waveforms are found to exhibit shock wave characteristics for most unstable engine operating conditions. It is shown that the characteristics of the resulting instability are independent of the nature of the initial disturbance and they depend solely upon the engine operating conditions and the characteristics of the unsteady combustion process.

Lores, M. E.↗

NASA NTRS1975

The mixture problem in computer mapping of terrain: Improved techniques for establishing spectral signature, atmospheric path radiance, and transmittance

The results of LANDSAT and Skylab research programs on the effects of the atmosphere on computer mapping of terrain include: (1) the concept of a ground truth map needs to be drastically revised; (2) the concept of training areas and test areas is not as simple as generally thought because of the problem of pixels that represent a mixture of terrain classes; (3) this mixture problem needs to be more widely recognized and dealt with by techniques of calculating spectral signatures of mixed classes, or by other methods; (4) atmospheric effects should be considered in computer mapping of terrain and in monitoring changes; and (5) terrain features may be used as calibration panels on the ground, from which atmospheric conditions can be determined and monitored. Results are presented of a test area in mountainous terrain of south-central Colorado for which an initial classification was made using simulated mixture-class spectral signatures and actual LANDSAT-1-MSS data.

Smedes, H. W.↗

NASA NTRS1977

Wide range operation of advanced low NOx combustors for supersonic high-altitude aircraft gas turbines

An initial rig program tested the Jet Induced Circulation (JIC) and Vortex Air Blast (VAB) systems in small can combustor configurations for NOx emissions at a simulated high altitude, supersonic cruise condition. The VAB combustor demonstrated the capability of meeting the NOx goal of 1.0 g NO2/kg fuel at the cruise condition. In addition, the program served to demonstrate the limited low-emissions range available from the lean, premixed combustor. A follow-on effort was concerned with the problem of operating these lean, premixed combustors with acceptable emissions at simulated engine idle conditions. Various techniques have been demonstrated that allow satisfactory operation on both the JIC and VAB combustors at idle with CO emissions below 20 g/kg fuel. The VAB combustor was limited by flashback/autoignition phenomena at the cruise conditions to a pressure of 8 atmospheres. The JIC combustor was operated up to the full design cruise pressure of 14 atmospheres without encountering an autoignition limitation although the NOx levels, in the 2-3 g NO2/kg fuel range, exceeded the program goal.

Roberts, P. B.↗

NASA NTRS1980

The primer vector in linear, relative-motion equations

Primer vector theory is used in analyzing a set of linear, relative-motion equations - the Clohessy-Wiltshire equations - to determine the criteria and necessary conditions for an optimal, N-impulse trajectory. Since the state vector for these equations is defined in terms of a linear system of ordinary differential equations, all fundamental relations defining the solution of the state and costate equations, and the necessary conditions for optimality, can be expressed in terms of elementary functions. The analysis develops the analytical criteria for improving a solution by (1) moving any dependent or independent variable in the initial and/or final orbit, and (2) adding intermediate impulses. If these criteria are violated, the theory establishes a sufficient number of analytical equations. The subsequent satisfaction of these equations will result in the optimal position vectors and times of an N-impulse trajectory. The solution is examined for the specific boundary conditions of (1) fixed-end conditions, two-impulse, and time-open transfer; (2) an orbit-to-orbit transfer; and (3) a generalized rendezvous problem. A sequence of rendezvous problems is solved to illustrate the analysis and the computational procedure.

Source record↗

NASA NTRS1983

Experiments using atmospheric forcing from a FGGE analysis to drive an upper ocean model

Results are presented of a series of numerical experiments in which an upper ocean model is driven by surface heat fluxes and stress fields derived from the FGGE SOP-1 GLAS analysis/forecast system (Halem, et al., 1982). The model results show that most changes in the mixed layer height and horizontal velocity occurs in the first days. On the other hand, changes in the temperature field take a longer time to develop. In the most realistic case (real initial conditions, instantaneous forcing fields from the atmospheric analysis), the resulting changes in temperature were larger than observed and the correlation between observed and predicted changes was poor. The deficiency in the forecast of SST changes may be due to several factors: lack of sufficient ocean resolution, improper initialization, lack of feedback between the ocean and the atmosphere and the absence of transports by the strong boundary currents, and perhaps unrealistic surface fluxes of heat and momentum. Unless these problems are alleviated it is not reasonable to perform coupled atmospheric ocean forecasts.

Camerlengo, A.↗

NASA NTRS2003

A Physical Validation Program for the GPM Mission

The GPM mission is currently planned for start in the late 2007 - early 2008 time frame. Its main scientific goal is to help answer pressing scientific problems arising within the context of global and regional water cycling. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like core satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe -- continuously. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. A major requirement before the retrieved rainfall information generated by the GPM mission can be used effectively by prognostic models to improve weather forecasts, hydrometeorological forecasts, and climate model reanalysis simulations is a capability to quantify the error characteristics of the retrievals. A solution for this problem has been upheld in past precipitation missions because of the lack of suitable error modeling systems incorporated into the validation programs and data distribution systems. An overview of how NASA intends to overcome this problem for the GPM mission using a physically-based error modeling approach within a multi-faceted validation program is described. The solution is to first identify specific user requirements and then determine the most stringent of these requirements that embodies all essential error characterization information needed by the entire user community. In the context of NASA s scientific agenda for the GPM mission, the most stringent user requirement is found within the data assimilation community. The fundamental theory of data assimilation vis-a-vis ingesting satellite precipitation information into the pre-forecast initializations is based on quantifying the conditional bias and precision errors of individual rain retrievals, and the space-time structure of the precision error (i.e., the spatial-temporal error covariance). By generating the hardware and software capability to produce this information in a near real-time fashion, and to couple the derived quantitative error properties to the actual retrieved rainrates, all key validation users can be satisfied. The talk will describe the essential components of the hardware and software systems needed to generate such near real-time error properties, as well as the various paradigm shifts needed within the validation community to produce a validation program relevant to the precipitation user community.

Smith, Eric A.↗

DOE OSTI2024

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

An experiment for electron-hadron scattering at the LHC

Abstract Novel considerations are presented on the physics, apparatus and accelerator designs for a future, luminous, energy frontier electron-hadron ( eh ) scattering experiment at the LHC in the thirties for which key physics topics and their relation to the hadron-hadron HL-LHC physics programme are discussed. Demands are derived set by these physics topics on the design of the LHeC detector, a corresponding update of which is described. Optimisations on the accelerator design, especially the interaction region (IR), are presented. Initial accelerator considerations indicate that a common IR is possible to be built which alternately could serve eh and hh collisions while other experiments would stay on hh in either condition. A forward-backward symmetrised option of the LHeC detector is sketched which would permit extending the LHeC physics programme to also include aspects of hadron-hadron physics. The vision of a joint eh and hh physics experiment is shown to open new prospects for solving fundamental problems of high energy heavy-ion physics including the partonic structure of nuclei and the emergence of hydrodynamics in quantum field theory while the genuine TeV scale DIS physics is of unprecedented rank.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Optimal polynomial based quantum eigenstate filtering with application to solving quantum linear systems

We present a quantum eigenstate filtering algorithm based on quantum signal processing (QSP) and minimax polynomials. The algorithm allows us to efficiently prepare a target eigenstate of a given Hamiltonian, if we have access to an initial state with non-trivial overlap with the target eigenstate and have a reasonable lower bound for the spectral gap. We apply this algorithm to the quantum linear system problem (QLSP), and present two algorithms based on quantum adiabatic computing (AQC) and quantum Zeno effect respectively. Both algorithms prepare the final solution as a pure state, and achieves the near optimal $\mathcal{\tilde{O}}$(dκlog(1/ϵ)) query complexity for a d -sparse matrix, where κ is the condition number, and ϵ is the desired precision. Neither algorithm uses phase estimation or amplitude amplification.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1978

UH cosmic rays - Possible origin in massive stars

The origin of the Z greater than 28, ultraheavy (UH) cosmic rays in supernova explosions of massive stars, at least about 10 solar masses, is considered. For Z greater than 70, the UH data are dominated by an r-process source distribution, but for the elements just beyond iron, Z from 29 to 36, the data cannot be explained by any single process of nucleosynthesis. This problem is solved naturally in a massive-star model by secondary neutron-capture reactions occurring during core helium burning (a limited s-process) and during explosive carbon burning. Interstellar-propagation calculations have been performed with these episodes of synthesis as source distributions, and the results offer an explanation for the current UH cosmic-ray data. Further, the heavy-element synthesis during explosive carbon burning is re-examined by using more realistic initial conditions given by the post-helium-burning configuration of the star. These results are compared with earlier work and the UH cosmic-ray data. Some effects of preferential acceleration, based upon ionization potential, are considered, and experimental tests for this model are discussed.

Wefel, J. P.↗

NASA NTRS1990

Improved boundary layer heat transfer calculations near a stagnation point

A thermal design of a solar receiver has been developed for the solutions of problems involving phase-change thermal energy storage and natural convection loss. Two dimensional axisymmetrical solidification and melting of materials contained between two concentric cylinders of finite length has been studied for thermal energy storage analysis. For calculation of free convection loss inside receiver cavity, two dimensional axisymmetrical, laminar, transient free convection including radiation effects has been studied using integral/finite difference method. Finite difference equations are derived for the above analysis subject to constant or variable material properties, initial conditions, and boundary conditions. The validity of the analyses has been substantiated by comparing results of the present general method with available analytic solutions or numerical results reported in the literature. Both explicit and implicit schemes are tested in phase change analysis with different number of nodes ranging from 4 to 18. The above numerical methods have been applied to the existing solar receiver analyzing computer code as additional subroutines. The results were computed for one of the proposed Brayton cycle solar receiver models running under the actual environmental conditions. Effect of thermal energy storage on the thermal behavior of the receiver has been estimated. Due to the thermal energy storage, about 65% reduction on working gas outlet temperature fluctuation has been obtained; however, maximum temperature of thermal energy storage containment has been increased about 18%. Also, effect of natural convection inside a receiver cavity on the receiver heat transfer has been analyzed. The finding indicated that thermal stratification occurs during the sun time resulting in higher receiver temperatures at the outlet section of the gas tube, and lower temperatures at the inlet section of the gas tube when compared with the results with no natural convection. Due to heat supply from the air during the shade time, minimum temperature has been increased, while maximum temperature has been reduced due to convection loss to air. Consequently, cyclic temperature fluctuation has been reduced 29% for working gas and 16% for thermal energy storage containment. On the other hand, despite the presence of the natural convection the time-averaged temperatures for receiver components were found to be similar for two cases with/without natural convection (maximum difference was 1.8%).

Ahn, Kyung Hwan↗