Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kriging”

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

A Statistics-Based Material Property Analysis to Support TPS Characterization

Accurate characterization of entry capsule heat shield material properties is a critical component in modeling and simulating Thermal Protection System (TPS) response in a prescribed aerothermal environment. The thermal decomposition of the TPS material during the pyrolysis and charring processes is poorly characterized and typically results in large uncertainties in material properties as inputs for ablation models. These material property uncertainties contribute to large design margins on flight systems and cloud re- construction efforts for data collected during flight and ground testing, making revision to existing models for entry systems more challenging. The analysis presented in this work quantifies how material property uncertainties propagate through an ablation model and guides an experimental test regimen aimed at reducing these uncertainties and characterizing the dependencies between properties in the virgin and charred states for a Phenolic Impregnated Carbon Ablator (PICA) based TPS. A sensitivity analysis identifies how the high-fidelity model behaves in the expected flight environment, while a Monte Carlo based uncertainty propagation strategy is used to quantify the expected spread in the in-depth temperature response of the TPS. An examination of how perturbations to the input probability density functions affect output temperature statistics is accomplished using a Kriging response surface of the high-fidelity model. Simulations are based on capsule configuration and aerothermal environments expected during the Mars Science Laboratory (MSL) entry sequence. We identify and rank primary sources of uncertainty from material properties in a flight-relevant environment, show the dependence on spatial orientation and in-depth location on those uncertainty contributors, and quantify how sensitive the expected results are.

Copeland, Sean R.↗

Creation of a Rapid High-Fidelity Aerodynamics Module for a Multidisciplinary Design Environment

In the traditional aerospace vehicle design process, each successive design phase is accompanied by an increment in the modeling fidelity of the disciplinary analyses being performed. This trend follows a corresponding shrinking of the design space as more and more design decisions are locked in. The correlated increase in knowledge about the design and decrease in design freedom occurs partly because increases in modeling fidelity are usually accompanied by significant increases in the computational expense of performing the analyses. When running high fidelity analyses, it is not usually feasible to explore a large number of variations, and so design space exploration is reserved for conceptual design, and higher fidelity analyses are run only once a specific point design has been selected to carry forward. The designs produced by this traditional process have been recognized as being limited by the uncertainty that is present early on due to the use of lower fidelity analyses. For example, uncertainty in aerodynamics predictions produces uncertainty in trajectory optimization, which can impact overall vehicle sizing. This effect can become more significant when trajectories are being shaped by active constraints. For example, if an optimal trajectory is running up against a normal load factor constraint, inaccuracies in the aerodynamic coefficient predictions can cause a feasible trajectory to be considered infeasible, or vice versa. For this reason, a trade must always be performed between the desired fidelity and the resources available. Apart from this trade between fidelity and computational expense, it is very desirable to use higher fidelity analyses earlier in the design process. A large body of work has been performed to this end, led by efforts in the area of surrogate modeling. In surrogate modeling, an up-front investment is made by running a high fidelity code over a Design of Experiments (DOE); once completed, the DOE data is used to create a surrogate model, which captures the relationships between input variables and responses into regression equations. Depending on the dimensionality of the problem and the fidelity of the code for which a surrogate model is being created, the initial DOE can itself be computationally prohibitive to run. Cokriging, a modeling approach from the field of geostatistics, provides a desirable compromise between computational expense and fidelity. To do this, cokriging leverages a large body of data generated by a low fidelity analysis, combines it with a smaller set of data from a higher fidelity analysis, and creates a kriging surrogate model with prediction fidelity approaching that of the higher fidelity analysis. When integrated into a multidisciplinary environment, a disciplinary analysis module employing cokriging can raise the analysis fidelity without drastically impacting the expense of design iterations. This is demonstrated through the creation of an aerodynamics analysis module in NASA’s OpenMDAO framework. Aerodynamic analyses including Missile DATCOM, APAS, and USM3D are leveraged to create high fidelity aerodynamics decks for parametric vehicle geometries, which are created in NASA’s Vehicle Sketch Pad (VSP). Several trade studies are performed to examine the achieved level of model fidelity, and the overall impact to vehicle design is quantified.

Aerodynamic Properties↗

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi ]Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post ]landing segment of the Flight Software (FSW). The vehicle crew is the prime input for touchdown detection, followed by an autonomous FSW algorithm, and finally a strictly time based backup timer. RCS thrusters must be inhibited before submersion in water to protect against possible damage due to firing these jets under water. In addition, neglecting to declare touchdown will not allow the vehicle to transition to post ]landing activities such as activating the Crew Module Up ]righting System (CMUS), resulting in possible loss of communication and difficult recovery. A previous AIAA paper gAssessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module h concluded that a strictly Inertial Measurement Unit (IMU) based detection method using an acceleration spike algorithm had the highest safety margins and shortest detection times of other methods considered. That study utilized finite element simulations of vehicle splashdown, generated by LS ]DYNA, which were expanded to a larger set of results using a Kriging surface fit. The study also used the Decelerator Systems Simulation (DSS) to generate flight dynamics during vehicle descent under parachutes. Proto ]type IMU and FSW MATLAB models provided the basis for initial algorithm development and testing. This paper documents an in ]depth trade study, using the same dynamics data and MATLAB simulations as the earlier work, to further develop the acceleration detection method. By studying the combined effects of data rate, filtering on the rotational acceleration correction, data persistence limits and values of acceleration thresholds, an optimal configuration was determined. The lever arm calculation, which removes the centripetal acceleration caused by vehicle rotation, requires that the vehicle angular acceleration be derived from vehicle body rates, necessitating the addition of a 2nd order filter to smooth the data. It was determined that using 200 Hz data directly from the vehicle IMU outperforms the 40 Hz FSW data rate. Data persistence counter values and acceleration thresholds were balanced in order to meet desired safety and performance. The algorithm proved to exhibit ample safety margin against early detection while under parachutes, and adequate performance upon vehicle splashdown. Fall times from algorithm initiation were also studied, and a backup timer length was chosen to provide a large safety margin, yet still trigger detection before CMUS inflation. This timer serves as a backup to the primary acceleration detection method. Additionally, these parameters were tested for safety on actual flight test data, demonstrating expected safety margins.

Daum, Jared↗

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi-Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post-landing segment of the Flight Software (FSW). An in-depth trade study was completed to determine optimal characteristics of the touchdown detection method resulting in an algorithm monitoring filtered, lever-arm corrected, 200 Hz Inertial Measurement Unit (IMU) vehicle acceleration magnitude data against a tunable threshold using persistence counter logic. Following the design of the algorithm, high fidelity environment and vehicle simulations, coupled with the actual vehicle FSW, were used to tune the acceleration threshold and persistence counter value to result in adequate performance in detecting touchdown and sufficient safety margin against early detection while descending under parachutes. An analytical approach including Kriging and adaptive sampling allowed for a sufficient number of finite element analysis (FEA) impact simulations to be completed using minimal computation time. The combination of a persistence counter of 10 and an acceleration threshold of approximately 57.3 ft/s2 resulted in an impact performance factor of safety (FOS) of 1.0 and a safety FOS of approximately 2.6 for touchdown declaration. An RCS termination acceleration threshold of approximately 53.1 ft/s(exp)2 with a persistence counter of 10 resulted in an increased impact performance FOS of 1.2 at the expense of a lowered under-parachutes safety factor of 2.2. The resulting tuned algorithm was then tested on data from eight Capsule Parachute Assembly System (CPAS) flight tests, showing an experimental minimum safety FOS of 6.1. The formulated touchdown detection algorithm will be flown on the Orion MPCV FSW during the Exploration Flight Test 1 (EFT-1) mission in the second half of 2014.

Daum, Jared↗

Vibroacoustic Tailoring of a Rod-Stiffened Composite Fuselage Panel with Multidisciplinary Considerations

An efficient multi-objective design tailoring procedure seeking to improve the vibroacoustic performance of a fuselage panel while maintaining or reducing weight is presented. The structure considered is the pultruded rod stitched efficient unitized structure, a highly integrated composite structure concept designed for a noncylindrical, next-generation flight vehicle fuselage. Modifications to a baseline design are evaluated within a six-parameter design space including spacing, flange width, and web height for both frame and stringer substructure components. The change in sound power radiation attributed to a design change is predicted using finite-element models sized and meshed for analyses in the 500 Hz, 1 kHz, and 2 kHz octave bands. Three design studies are carried out in parallel while considering a diffuse acoustic field excitation and two types of turbulent boundary-layer excitation. Kriging surrogate models are used to reduce the computational costs of resolving the vibroacoustic and weight objective Pareto fronts. The resulting Pareto optimal designs are then evaluated under a static pressurization ultimate load to assess structural strength and stability. Results suggest that choosing alternative configurations within the considered design space can reduce weight and improve vibroacoustic performance without compromising strength and stability of the structure under the static load condition considered, but the tradeoffs are significantly influenced by the spatial characteristics of the assumed excitation field.

Allen, Albert R.↗

Optimal Reference Strain Structure for Studying Dynamic Responses of Flexible Rockets

In the proposed paper, the optimal design of reference strain structures (RSS) will be performed targeting for the accurate observation of the dynamic bending and torsion deformation of a flexible rocket. It will provide the detailed description of the finite-element (FE) model of a notional flexible rocket created in MSC.Patran. The RSS will be attached longitudinally along the side of the rocket and to track the deformation of the thin-walled structure under external loads. An integrated surrogate-based multi-objective optimization approach will be developed to find the optimal design of the RSS using the FE model. The Kriging method will be used to construct the surrogate model. For the data sampling and the performance evaluation, static/transient analyses will be performed with MSC.Natran/Patran. The multi-objective optimization will be solved with NSGA-II to minimize the difference between the strains of the launch vehicle and RSS. Finally, the performance of the optimal RSS will be evaluated by checking its strain-tracking capability in different numerical simulations of the flexible rocket.

finite-element (FE)↗

Design of Distributed Mail-Slot Propulsion System on a Hybrid Wingbody Aircraft

This paper deals with designing a thrust distribution strategy when a Turboelectric Distributed Propulsion (TeDP) system of 16 embedded propulsors is installed on an aerodynamically optimized hybrid wing-body configuration. This HWB previously designed to satisfy conditions of trim, longitudinally static stability and specific cargo space is employed as the baseline configuration for the current study of seeking an optimal propulsion/power system. According to the nature of the entrance flow condition for each distributed propulsion passage in hybrid wing-body aircraft, the ingested boundary layer thickness differs and results in different propulsive reaction. An optimal distribution of thrust and power output is determined by how the system utilizes the propulsive characteristics of each passage. The design space and the number of design variables are selected and described accordingly. An actuator disk model is employed to model thrust generation and shaft power from the propulsor. To carry out the optimization of the propulsion/power system on a computationally expensive CFD model, a Kriging method in conjunction with a Genetic Algorithm (GA) is applied. Throughout the design process, the propulsion performances of the sampled propulsion/power system are analyzed and compared to those of a clean flow engine. The performance metrics includes mass flow rate, fan pressure ratio besides the thrust and shaft power. Minimization of total shaft power from the distributed engine is performed at multiple thrust levels. The benefit of boundary layer ingestion propulsion system is quantified via comparison of thrust equivalent, shaft power and mass flow equivalent clean flow engines with CFD based system design.

Hybrid Wingbody↗

New Approach to Characterize CO2 and CH4 Emissions over Sacramento, California Using an Airborne Aircraft Measurement

Emission fluxes of carbon dioxide (CO2) and methane (CH4) over Sacramento, California, USA, are estimated using an aircraft equipped with a cavity ring-down greenhouse gas sensor through the Alpha Jet Atmospheric eXperiment (AJAX). To better quantify the emission fluxes over the entire city and two small point sources within the area, we designed flights in a cylindrical pattern and computed the emission fluxes from 9 flights over all seasons using the Kriging method based on Gausss divergence theorem. The concentrations at the downwind side of Sacramento show somewhat persistent patterns among the 9 flights, but the magnitudes and locations of the fluxes vary depending on the individual weather conditions and seasonality on a given flight day. The net emission fluxes tend to be smaller in summer (July) than winter (November). Both trace gas mixing ratios as well as the wind speed and direction show high spatial variability both horizontally and vertically. The local maximum emission fluxes over a landfill and a rice field measured on July 29, 2015, were 2-3 times lower than the maximum measured over the entire city. These low emission fluxes suggest that not only the local source but also the transport of greenhouse gases due to diurnal and seasonal variation in winds have a strong effect on the CO2 and CH4 flux estimation of the both local emissions and city-size emissions. Our results highlight that the aircraft-based approach is effective and useful for capturing city-size emission fluxes and estimating the greenhouse gas emissions and its sources.

airborne↗

Quantification of CO2 and CH4 emissions over Sacramento, California, based on divergence theorem using aircraft measurements

Emission estimates of carbon dioxide (CO2) and methane (CH4) and the meteorological factors affecting them are investigated over Sacramento, California, using an aircraft equipped with a cavity ring-down greenhouse gas sensor as part of the Alpha Jet Atmospheric eXperiment (AJAX) project. To better constrain the emission fluxes, we designed flights in a cylindrical pattern and computed the emission fluxes from two flights using a kriging method and Gauss’s divergence theorem. Differences in wind treatment and assumptions about background concentrations affect the emission estimates by a factor of 1.5 to 7. The uncertainty is also impacted by meteorological conditions and distance from the emission sources. The vertical layer averaging affects the flux estimate, but the choice of raw wind or mass-balanced wind is more important than the thickness of the vertical averaging for massbalanced wind for both urban and local scales. The importance of vertical mass transfer for flux estimates is examined, and the difference in the total emission estimate with and without vertical mass transfer is found to be small, especially at the local scale. The total flux estimates accounting for the entire circumference are larger than those based solely on measurements made in the downwind region. This indicates that a closed-shape flight profile can better contain total emissions relative to a one-sided curtain flight because most cities have more than one point source and wind direction can change with time and altitude. To reduce the uncertainty of the emission estimate, it is important that the sampling strategy account not only for known source locations but also possible unidentified sources around the city. Our results highlight that aircraft-based measurements using a closed-shape flight pattern are an efficient and useful strategy for identifying emission sources and estimating local- and city-scale greenhouse gas emission fluxes.

Ju-Mee Ryoo↗

Multi-Fidelity Heating Prediction of Adaptable, Deployable Entry Placement Technology Vehicles

The objective of this work was to investigate a multi-fidelity modeling approach to accurately and efficiently predict the laminar and turbulent convective heating on adaptable, deployable entry placement technology vehicles in Mars entry. A previously developed co-Kriging based multi-fidelity modeling approach was used to model the laminar and turbulent convective heat fluxes at several surface locations along the vehicle, including the rib sections. The laminar convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 3% when compared to high-fidelity CFD simulations. The turbulent convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 8% when compared to high-fidelity CFD simulations. Compared to a single fidelity model, the multi-fidelity model required approximately one-third the number of high-fidelity model evaluations to obtain the same accuracy level. The computational cost of evaluating the multi-fidelity model was approximately five orders of magnitude less than one high-fidelity model simulation.

Mario Santos↗

The INSTEP Monitoring Network: Merging High-and-Low Cost Measurements to Characterize California Wildfires

Despite challenges with data quality and scope, low-cost sensor networks have skyrocketed in popularity over the last 15 years, making air quality data available on refined spatial scales. More recently, studies have leveraged both high and low-quality instruments to create stronger “hybrid” models, with most studies focusing on particulate matter. Low-cost measurements typically represent ground-level emissions only, providing context for human health issues from climate change-driven events such as wildfires. Since low-cost sensors’ capabilities are localized, daily events and microclimates tend to dominate the data rather than larger regional or atmospheric trends. Likewise, their low cost explains their high uncertainty. In contrast, some regulatory-grade instruments produce column measurements as well, providing reliable information on a broader scope. To bridge this gap while expanding into gas-phase measurements, we deployed 12 air quality sensor packages in California, USA during the 2022 wildfire season. These INSTEP (Inexpensive Network Sensor Technology Exploring Pollution) monitors measure carbon monoxide (CO), carbon dioxide (CO2), ozone (O3), nitrogen dioxide (NO2), and several hydrocarbons including methane (CH4) and formaldehyde (HCHO). Half of the monitors were co-located with remote sensing spectrometers: NASA Pandora and Total Column Carbon Observing Network (TCCON). The overlap in pollutants includes NO2, O3, and HCHO between the INSTEP monitors and the Pandora column measurements. TCCON covers column CO, CO2, and CH4, rounding out our comparison. Most of the monitors were distributed throughout the San Francisco Bay area, and an additional three were located within 100 km of Los Angeles. The sites ranged in geographic and population characteristics, including desert, mountainous, coastal, and urban locations. Since varying environmental conditions such as temperature and pressure are known to challenge sensor performance, we will apply newer sensor “calibration” techniques meant to combat this. We will normalize our sensor signals by z-scoring them prior to applying a single calibration model in the form of multivariate linear regression or an artificial neural network. While this technique has been validated for the hydrocarbon and ozone sensor types (metal oxide), it has not yet been tested on electrochemical and non-dispersive infrared sensors, which are also used in the INSTEP monitors. This will serve as a test to see if this normalization technique – or another – is most effective in accounting for environmental differences among sensors. Related data analysis efforts have found success with a variety of geospatial analysis techniques, including weighted network models in which high-quality instruments are given higher weights than their low-cost counterparts. Our preliminary analysis will focus on kriging, which uses a Gaussian algorithm to assign weights, providing estimated pollution levels at locations between monitors. Smoke trajectory and evolution will also be considered using both measurement types. We also aim to baseline subtract our emission estimates from each region to determine which portion of emissions are regional and local, further characterizing burn differences in northern and southern California fires. Future directions include using INSTEP jointly with TEMPO satellite data, and mobile deployments on aircraft and uncrewed aerial vehicles (UAV).

Low-cost sensors↗

Sampling of Space-Based Observations of XCO2 Associated With Biomass Burning Events

- The OCO-2 Level 3 assimilation product fills in gaps where OCO-2 does not have observations. - Although the Level 3 clearly shows the increased CO2 over the Amazon associated with the September 2017 biomass burning event, the total amount of CO2 may be underestimated because of a sampling bias. - The TCCON observations match fairly well with the OCO-2 observations and can be used to estimate the accuracy of the OCO observations and the temporal sampling bias. However, since it is a remotely sensed product, it may suffer the same sampling biases associated with biomass burning events as OCO-2. - Although aircraft data show that there is a significant variation in xCO2 with altitude, they suggest the OCO-2 fused products created by kriging may be better than the assimilated product at detecting large deviations from the average state on daily and regional scales.

Thomas J. Hearty, III↗

Efficient Global Optimization with Gradient Finish for Design Under Uncertainty

The Efficient Global Optimization (EGO) algorithm is extended to include a gradient-descent-based finish upon reaching a threshold value of the expected improvement function. Emphasis is placed on efficient evaluation of local gradients using Kriging models during the gradient-based finish to enable application to design under uncertainty (DUU) problems. The modified algorithm is applied to both the well-known Rosenbrock function and a more challenging hypersonic inlet design under uncertainty problem. Results demonstrate improvement in locating the global optimum compared to the classical implementation of EGO, as well as a reduced number of true function evaluations compared to pure gradient-based algorithms. For the Rosenbrock function, a global optimum is returned using an average of 12% fewer function calls than a gradient-based optimizer with comparable tolerances. For the design under uncertainty problem, a global optimum is found using an average of 75% fewer function calls than a gradient-based optimizer. Global Pareto fronts of multiobjective DUU problems are obtained at little additional cost after a single optimization is complete.

Optimization↗

Sampling Biases of Space-Based Observations of XCO2 Associated With South American Biomass Burning Events During 2017-2020

- We examine time series of OCO-2 xCO 2 observations coincident with other remotely sensed and in-situ observations to better understand external sampling effects. - The OCO-2 Level 3 assimilation product fills in gaps where OCO-2 does not have observations. - Although the Level 3 clearly shows the increased CO 2 over the Amazon associated with the September 2017 biomass burning event, the total amount of CO 2 may be underestimated because of a sampling bias. - Since TCCON is also a remotely sensed product, it may suffer the same sampling biases associated with biomass burning events as OCO-2. - Aircraft data show that there is significant variation in xCO 2 with altitude. - The OCO-2 fused products created by kriging may be better than the assimilated product at detecting deviations from the average state.

OCO↗

PRODEM: An Annual Series of Summer DEMs (2019 through 2022) of the Marginal Areas of the Greenland Ice Sheet

Surface topography across the marginal zone of the Greenland Ice Sheet is constantly evolving in response to changing weather, season, climate, and ice dynamics. However, current digital elevation models (DEMs) for the ice sheet are usually based on data from a multi-year period, thus obscuring these changes over time. Here we present four 500 m resolution summer DEMs (PRODEMs) of the Greenland Ice Sheet marginal zone for 2019 through 2022. The PRODEMs cover the marginal zone from the ice edge to 50 km inland, hence capturing all Greenland outlet glaciers. Each PRODEM is based on data fusion of CryoSat-2 radar altimetry and ICESat-2 laser altimetry using regionally varying kriging of elevation anomalies relative to ArcticDEM. The PRODEMs are validated using leave-one-out cross-validation, and PRODEM19 is further validated against an external data set, showcasing their ability to correctly represent surface elevations within the associated spatially varying prediction uncertainties. We observe a general lowering of surface elevations during the 4-year PRODEM period, but the spatial pattern of change is highly complex and with annual changes superimposed. The PRODEMs enable detailed studies of the marginal ice sheet elevation changes. With their high spatio-temporal resolution, the PRODEMs will be of value to a wide range of researchers and users studying ice sheet dynamics and monitoring how the ice sheet responds to changing environmental conditions. PRODEMs from summer 2019 through 2022 are available at https://doi.org/10.22008/FK2/52WWHG (Winstrup, 2024), and we plan to annually update the product henceforth.

Digital Elevation Models↗

Development of Versatile Compressor Modeling using Approximation Techniques for Alternative Refrigerants Evaluation

Refrigerants are the life-blood of vapor compression systems that are widely used in Heating, Ventilation, Air-Conditioning, and Refrigeration (HVAC&R) applications. The HVAC&R community is currently transitioning from main-stream refrigerants that have high Global Warming Potential (GWP) to alternative lower-GWP refrigerants. During this transition, it is important to account for the life cycle climate performance of alternative refrigerants since their performance will be different than that of higher-GWP refrigerants. This requires the evaluation of the system performance with the new refrigerants. Unfortunately, it is extremely difficult to predict the realistic performance of new alternative refrigerants without experimental validation. One of the main challenges in this regard is modeling the compressor performance with high fidelity due to the complex interaction of operating parameters, geometry, boundary conditions, and fluid properties. High fidelity compressor models are computationally expensive and require significant pre-processing to evaluate the performance of alternative refrigerants. This paper presents a new approach to modeling compressor performance when alternative refrigerants are used. The new modeling concept relies on using existing compressor performance to create an approximate model that captures the dependence of compressor performance on key operating parameters and fluid properties. The model can be built using a myriad of approximation techniques. This paper focuses on Kriging-based techniques to develop higher fidelity approximate compressor models. Baseline and at least one alternative refrigerant performance data are used to build the model. The model accuracy was evaluated by comparing the model results with compressor performance data using other refrigerants. Preliminary results show that the approximate model can predict the compressor mass flow rate and power consumption within 5%.

Approximation↗

MULTI-SCALE MODELING AND APPROXIMATION ASSISTED OPTIMIZATION OF BARE TUBE HEAT EXCHANGERS

Air-to-refrigerant heat exchangers are very common in air-conditioning, heat pump and refrigeration applications. In these heat exchangers, there is a great benefit in terms of size, weight, refrigerant charge and heat transfer coefficient, by moving from conventional channel sizes (~ 9mm) to smaller channel sizes (< 5mm). This work investigates new designs for air-to-refrigerant heat exchangers with tube outer diameter ranging from 0.5 to 2.0mm. The goal of this research is to develop and optimize the design of these heat exchangers and compare their performance with existing state of the art designs. The air-side performance of various tube bundle configurations are analyzed using a Parallel Parameterized CFD (PPCFD) technique. PPCFD allows for fast-parametric CFD analyses of various geometries with topology change. Approximation techniques drastically reduce the number of CFD evaluations required during optimization. Maximum Entropy Design method is used for sampling and Kriging method is used for metamodeling. Metamodels are developed for the air-side heat transfer coefficients and pressure drop as a function of tube-bundle dimensions and air velocity. The metamodels are then integrated with an air-to-refrigerant heat exchanger design code. This integration allows a multi-scale analysis of air-side performance heat exchangers including air-to-refrigerant heat transfer and phase change. Overall optimization is carried out using a multi-objective genetic algorithm. The optimal designs found can exhibit 50 percent size reduction, 75 percent decrease in air side pressure drop and doubled air heat transfer coefficients compared to a high performance compact micro channel heat exchanger with same capacity and flow rates.

CFD↗

On the benefit of using both hand calculation methods and computer experiments algorithms to solve complex nuclear criticality safety problems - 14299

A companion paper shows the 'anomalies', and even the safety issues, that may be raised by the very usual problem of making safe a fissile mass inside a controlled geometry. These anomalies mainly come from the solving of such problems which belong to the mathematical family of 'robust inversion' including two inversion parameters (the fissile mass and the geometry) and one nuisance parameter (the moderation ratio of the fissile material). Since a 'brute force' method of resolution requires a lot of calculations and since a 'smart' method of resolution requires a lot of engineer time, it is tempting to solve these problems by a 'custom made' approach rather than to achieve a 'whole' resolution. The present paper shows how one can use a hand calculation method (derivation of the one-speed diffusion theory) to perform the whole resolution of such kind of problems with a minimum calculations and human efforts. The powerful capability of predicting the 'anomalies' highlighted in the companion paper is shown. Then the limitations of such a method are discussed: limitation in the geometry of the problem and bias in absolute values prediction. Finally, new algorithms of computer experiments design, dedicated to solve robust inversion problems and based on Kriging surrogates, are tested to solve these problems with limiting calculations and human efforts. The performance of these new algorithms may supplement the one-speed diffusion theory limitations. (authors)

42 ENGINEERING↗