Search NASA⌕ Search

SEARCH · Search NASA

Results for “convection permitting 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.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Spherical Diffusion Flames: Structure and Dynamics

The spherical geometry is the most suitable one to study fundamental issues concerning the structure and the dynamics of diffusion flames. From a theoretical point of view, it is the only geometry that permits the existence of a truly one dimensional stationary diffusion flame. A stationary planar diffusion flame with the fuel supplied upstream, at x approaches -infinity say, and the oxidant downstream, at x approaches +infinity is not possible. For a steady diffusion flame to exist, one must have nonzero fluxes of fuel and oxidant towards the flame. However, in the unlimited region behind the planar flame the only bounded solutions to the reaction-free convective-diffusive operator are constants. Hence the oxidant concentration behind the flame remains constant and there is no mechanism to generate the necessary flux towards the flame. A one-dimensional problem can be formulated if the reactants are supplied at finite locations; but the boundary conditions in this case introduce unnecessary complications and do not appropriately model the physical reality. Indeed, a planar diffusion flame can be established in the stagnation-point flow of two opposed jets but the flame in this case is stretched and the flow is essentially two-dimensional. The only stationary one-dimensional diffusion flame in an unlimited environment is therefore the spherical flame.

Matalon, Moshe↗
NASA NTRS2013

The NASA-Goddard Multi-Scale Modeling Framework - Land Information System: Global Land/atmosphere Interaction with Resolved Convection

The present generation of general circulation models (GCM) use parameterized cumulus schemes and run at hydrostatic grid resolutions. To improve the representation of cloud-scale moist processes and landeatmosphere interactions, a global, Multi-scale Modeling Framework (MMF) coupled to the Land Information System (LIS) has been developed at NASA-Goddard Space Flight Center. The MMFeLIS has three components, a finite-volume (fv) GCM (Goddard Earth Observing System Ver. 4, GEOS-4), a 2D cloud-resolving model (Goddard Cumulus Ensemble, GCE), and the LIS, representing the large-scale atmospheric circulation, cloud processes, and land surface processes, respectively. The non-hydrostatic GCE model replaces the single-column cumulus parameterization of fvGCM. The model grid is composed of an array of fvGCM gridcells each with a series of embedded GCE models. A horizontal coupling strategy, GCE4fvGCM4Coupler4LIS, offered significant computational efficiency, with the scalability and I/O capabilities of LIS permitting landeatmosphere interactions at cloud-scale. Global simulations of 2007e2008 and comparisons to observations and reanalysis products were conducted. Using two different versions of the same land surface model but the same initial conditions, divergence in regional, synoptic-scale surface pressure patterns emerged within two weeks. The sensitivity of largescale circulations to land surface model physics revealed significant functional value to using a scalable, multi-model land surface modeling system in global weather and climate prediction.

cumulus schemes↗
DOE OSTI2021

Growth of Mesoscale Convective Systems in Observations and a Seasonal Convection-Permitting Simulation over Argentina

Abstract A 6.5-month, convection-permitting simulation is conducted over Argentina covering the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) field campaign and is compared with observations to evaluate mesoscale convective system (MCS) growth prediction. Observed and simulated MCSs are consistently identified, tracked, and separated into growth, mature, and decay stages using top-of-the-atmosphere infrared brightness temperature and surface rainfall. Simulated MCS number, lifetime, seasonal and diurnal cycles, and various cloud-shield characteristics including growth rate are similar to those observed. However, the simulation produces smaller rainfall areas, greater proportions of heavy rainfall, and faster system propagations. Rainfall area is significantly underestimated for long-lived MCSs but not for shorter-lived MCSs, and rain rates are always overestimated. These differences result from a combination of model and satellite retrieval biases, in which simulated MCS rain rates are shifted from light to heavy, while satellite-retrieved rainfall is too frequent relative to rain gauge estimates. However, the simulation reproduces satellite-retrieved MCS cloud-shield evolution well, supporting its usage to examine environmental controls on MCS growth. MCS initiation locations are associated with removal of convective inhibition more than maximized low-level moisture convergence or instability. Rapid growth is associated with a stronger upper-level jet (ULJ) and a deeper northwestern Argentinean low that causes a stronger northerly low-level jet (LLJ), increasing heat and moisture fluxes, low-level vertical wind shear, baroclinicity, and instability. Sustained growth corresponds to similar LLJ, baroclinicity, and instability conditions but is less sensitive to the ULJ, large-scale vertical motion, or low-level shear. Growth sustenance controls MCS maximum extent more than growth rate.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Analytical and numerical studies of the thermocapillary flow in a uniformly floating zone

The microgravity environment of an orbiting vehicle permits crystal growth experiments in the presence of greatly reduced buoyant convection in the liquid melt. Crystals grown in ground-based laboratories do not achieve their potential properties because of dopant variations caused by flow in the melt. The floating zone crystal growing system is widely used to produce crystals of silicon and other materials. However, in this system the temperature gradient on the free sidewall surface of the melt is the source of a thermocapillary flow which does not disappear in the low-gravity environment. The idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow to the melt sidewall leaving the interior of the melt passive is examined. A cylinder of fluid with an axial temperature gradient imposed on the cylindrical sidewall is considered. A half zone and the linearized, axisymmetric flow in the absence of crystal growth is examined. Rotation is found to confine the linear thermocapillary flow. A simplified model is extended to a full zone and both linear and nonlinear thermocapillary flows are studied theoretically. Analytical and numerical methods are used for the linear flows and numerical methods for the nonlinear flows. It was found that the linear flows in the full zone have more complicated and thicker boundary layer structures than in the half zone, and that these flows are also confined by the rotation. However, for the simplified model considered and for realistic values for silicon, the thermocapillary flow is not linear. The fully nonlinear flow is strong and unsteady (a weak oscillation is present) and it penetrates the interior. Some non-rotating flow results are also presented. Since silicon as a large value of thermal conductivity, one would expect the temperature fields to be determined by conduction alone. This is true for the linear and weakly nonlinear flows, but for the stronger nonlinear flow the results show that temperature advection is also important. Uniform rotation may still be a means of confining the flow and the results obtained define the procedure to be used to examine this hypothesis.

Fowlis, W. W.↗
NASA NTRS1992

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗
DOE OSTI2024

Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model

Abstract. Emissions of anthropogenic aerosol and their precursors are often prescribed in global aerosol models. Most of these emissions are spatially heterogeneous at model grid scales. When remapped from low-resolution data, the spatial heterogeneity in emissions can be lost, leading to large errors in the simulation. It can also cause the conservation problem if non-conservative remapping is used. The default anthropogenic emission treatment in the Energy Exascale Earth System Model (E3SM) is subject to both problems. In this study, we introduce a revised emission treatment for the E3SM Atmosphere Model (EAM) that ensures conservation of mass fluxes and preserves the original emission heterogeneity at the model-resolved grid scale. We assess the error estimates associated with the default emission treatment and the impact of improved heterogeneity and mass conservation in both globally uniform standard-resolution (∼ 165 km) and regionally refined high-resolution (∼ 42 km) simulations. The default treatment incurs significant errors near the surface, particularly over sharp emission gradient zones. Much larger errors are observed in high-resolution simulations. It substantially underestimates the aerosol burden, surface concentration, and aerosol sources over highly polluted regions, while it overestimates these quantities over less-polluted adjacent areas. Large errors can persist at higher elevation for daily mean estimates, which can affect aerosol extinction profiles and aerosol optical depth (AOD). We find that the revised treatment significantly improves the accuracy of the aerosol emissions from surface and elevated sources near sharp spatial gradient regions, with significant improvement in the spatial heterogeneity and variability of simulated surface concentration in high-resolution simulations. In the next-generation E3SM running at convection-permitting scales where the resolved spatial heterogeneity is significantly increased, the revised emission treatment is expected to better represent the aerosol emissions as well as their lifecycle and impacts on climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Exploring Western North Pacific Tropical Cyclone Activity in the High‐Resolution Community Atmosphere Model

Abstract High‐resolution climate models (∼28 km grid spacing) can permit realistic simulations of tropical cyclones (TCs), thus enabling their investigation in relation to the climate system. On the global scale, previous works have demonstrated that the Community Atmosphere Model (CAM) version 5 presents a reasonable TC climatology under prescribed present‐day (1980–2005) forcing. However, for the Western North Pacific (WNP) region, known biases in simulated TC genesis frequency and location under‐represent the basin's dominant share in observations. This study addresses these model biases in WNP by evaluating WNP TCs in a decadal simulation, and exploring potential improvements through nudging experiments. Among the major environmental controls of TC genesis, the lack of mid‐level moisture is identified as the leading cause of the deficit in simulated WNP TC genesis over the Pacific Warm Pool. Subsequent seasonal experiments explore the effect of constraining the large‐scale environment on TC development by nudging WNP temperature field toward reanalysis at various strengths. Temperature nudging elicits a significant response in TC genesis and intensity development, as well as in moisture and convection over the Warm Pool. These responses are sensitive to the choice of nudging timescale. Overall, the nudging experiments demonstrate that improvements in the large‐scale environment can lead to improvements in simulated TCs, suggesting future model developments in relation to model physics. In this way, the potential improvements in model fidelity will contribute to the understanding of how the mean state of current or future climates may give rise to extremes such as TCs.

58 GEOSCIENCES↗
NASA NTRS2010

Genesis of Pre-Hurricane Felix (2007): The Role of the Easterly Wave Critical Layer - Part 1

The formation of pre Hurricane Felix (2007) in a tropical easterly wave is examined in a two-part study using the Weather Research and Forecasting (WRF) model with a high-resolution nested grid configuration that permits the representation of cloud system processes. The simulation commences during the wave stage of the precursor African easterly-wave disturbance. Here the simulated and observed developments are compared, while in Part II of the study various large-scale analyses, physical parameterizations, and initialization times are explored to document model sensitivities. In this first part the authors focus on the wave/vortex morphology, its interaction with the adjacent intertropical convergence zone complex, and the vorticity balance in the neighborhood of the developing storm. Analysis of the model simulation points to a bottom-up development process within the wave critical layer and supports the three new hypotheses of tropical cyclone formation proposed recently by Dunkerton, Montgomery, and Wang. It is shown also that low-level convergence associated with the ITCZ helps to enhance the wave signal and extend the "wave pouch" from the jet level to the top of the atmospheric boundary layer. The region of a quasi-closed Lagrangian circulation within the wave pouch provides a focal point for diabatic merger of convective vortices and their vortical remnants. The wave pouch serves also to protect the moist air inside from dry air intrusion, providing a favorable environment for sustained deep convection. Consistent with the authors' earlier findings, the tropical storm forms near the center of the wave pouch via system-scale convergence in the lower troposphere and vorticity aggregation. Components of the vorticity balance are shown to be scale dependent, with the immediate effects of cloud processes confined more closely to the storm center than the overturning Eliassen circulation induced by diabatic heating, the influence of which extends to larger radii.

Wang, Zhuo↗
NASA NTRS1984

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗
NASA NTRS1992

Materials science research in microgravity

There are several important attributes of an extended duration microgravity environment that offer a new dimension in the control of the microstructure, processing, and properties of materials. First, when gravitational effects are minimized, buoyancy driven convection flows are also minimized. The flows due to density differences, brought about either by composition or temperature gradients will then be reduced or eliminated to permit a more precise control of the temperature and the composition of a melt which is critical in achieving high quality crystal growth of electronic materials or alloy structures. Secondly, body force effects such as sedimentation, hydrostatic pressure, and deformation are similarly reduced. These effects may interfere with attempts to produce uniformly dispersed or aligned second phases during melt solidification. Thirdly, operating in a microgravity environment will facilitate the containerless processing of melts to eliminate the limitations of containment for reactive melts. The noncontacting forces such as those developed from electromagnet, electrostatic, or acoustic fields can be used to position samples. With this mode of operation, contamination can be minimized to enable the study of reactive melts and to eliminate extraneous crystal nucleation so that novel crystalline structures and new glass compositions may be produced. In order to take advantage of the microgravity environment for materials research, it has become clear that reliable processing models based on a sound ground based experimental experience and an established thermophysical property data base are essential.

Perepezko, John H.↗
NASA NTRS1994

Baroclinic instability in the interiors of the giant planets: A cooling history of Uranus?

We propose a quasigeostrophic, baroclinic model for heat transport within the interior of a stably stratified Jovian planet, based on motion in thin cylindrical annuli. Density decreases from the center outward and is zero at the surface of the planet. In the homogeneous case (no core), we find instability for the poles hotter than the equator, but not for the reverse. If the motion is bounded by an impenetrable core, instability occurs for both cases. Much of the behavior can be explained by analogy to conventional baroclinic instability theory. Motivated by our results, we explore a possible connection between the highly inclined rotation axis of Uranus and its anomalously low surface heat flux. We assume that the planets formed hot. Our conjecture is that heat was efficiently convected outwards by baroclinic instability in Uranus (with the poles hotter than the equator), but not in the other three Jovian planets. The surface temperature was higher for the stably stratified case (Uranus), leading to a higher rate of infrared emission and faster cooling. Therefore, we propose that Uranus lost its internal heat sooner than Neptune because baroclinic motions, permitted by its inclination to the sun, were able to extract its internal heat while the surface was still warm.

Holme, Richard↗
NASA NTRS2018

Aerosol Midlatitude Cyclone Indirect Effects in Observations and High-Resolution Simulations

Aerosol-cloud interactions are a major source of uncertainty in inferring the climate sensitivity from the observational record of temperature. The adjustment of clouds to aerosol is a poorly constrained aspect of these aerosol-cloud interactions. Here, we examine the response of midlatitude cyclone cloud properties to a change in cloud droplet number concentration (CDNC). Idealized experiments in high-resolution, convection-permitting global aquaplanet simulations with constant CDNC are compared to 13 years of remote-sensing observations. Observations and idealized aquaplanet simulations agree that increased warm conveyor belt (WCB) moisture flux into cyclones is consistent with higher cyclone liquid water path (CLWP). When CDNC is increased a larger LWP is needed to give the same rain rate. The LWP adjusts to allow the rain rate to be equal to the moisture flux into the cyclone along the WCB. This results in an increased CLWP for higher CDNC at a fixed WCB moisture flux in both observations and simulations. If observed cyclones in the top and bottom tercile of CDNC are contrasted it is found that they have not only higher CLWP but also cloud cover and albedo. The difference in cyclone albedo between the cyclones in the top and bottom third of CDNC is observed by CERES to be between 0.018 and 0.032, which is consistent with a 4.6-8.3 Wm(exp -2) in-cyclone enhancement in upwelling shortwave when scaled by annual-mean insolation. Based on a regression model to observed cyclone properties, roughly 60% of the observed variability in CLWP can be explained by CDNC and WCB moisture flux.

Mccoy, Daniel T.↗
NASA NTRS1998

LASE measurements of convective boundary layer development during SGP97

The Southern Great Plains 1997 (SGP97) field experiment was conducted in Oklahoma during June-July 1997 to validate the models used for computing remote soil moisture using measurements by microwave radiometers. One of the objectives of SGP97 was to examine the effect of soil moisture on the evolution of the Atmospheric Boundary Layer (ABL) and clouds over the Southern Great Plains (SGP) during the warm season. The LASE (Lidar Atmospheric Sensing Experiment) airborne DIAL (Differential Absorption Lidar) system, which was flown autonomously on the NASA ER-2 aircraft during previous missions, was reconfigured to fly on the NASA P3 research aircraft. During SGP97 LASE was used to study the morning evolution of the ABL, particularly as manifested in the development of the convective boundary layer, and to study the influence of soil moisture variations on the development of ABL. The ABL development is strongly influenced by the surface energy budget, which is in turn influenced by soil moisture, mesoscale meteorology, clouds, and solar insolation. LASE data acquired during this mission are being used to study the ABL water vapor budget, the development of the ABL, spatial and temporal variabilities in the ABL, and the meteorological factors that influence the ABL development. This field experiment also permitted comparisons of LASE water vapor measurements with water vapor profiles acquired by radiosondes launched at the DOE (Department of Energy) Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site and at NASA/Wallops Flight Facility, as well as with measurements from other SGP97 aircraft.

Ismail, Syed↗
DOE OSTI2021

Five percent measurement of the gravitational constant in the Large Magellanic Cloud

We perform a novel test of General Relativity by measuring the gravitational constant in the Large Magellanic Cloud (LMC). The LMC contains six well-studied Cepheid variable stars in detached eclipsing binaries. Radial velocity and photometric observations enable a complete orbital solution, and precise measurements of the Cepheids’ periods permit detailed stellar modelling. Both are sensitive to the strength of gravity, the former via Kepler’s third law and the latter through the gravitational free-fall time. We jointly fit the observables for stellar parameters and the gravitational constant. Performing a full Markov Chain Monte Carlo analysis of the parameter space including all relevant nuisance parameters, we constrain the gravitational constant in the Large Magellanic Cloud relative to the Solar System to be G LMC /G SS = $0.93{^{+0.05}_{-0.04}}$. We discuss the implications of this 5% measurement of Newton’s constant in another galaxy for dark energy and modified gravity theories. This result excludes one Cepheid, CEP-1812, which is an outlier and needs further study: it is either a highly unusual system to which our model does not apply, or it prefers G LMC < G SS at 2.6σ. We also obtain new bounds on critical parameters that appear in semi-analytic descriptions of stellar processes. In particular, we measure the mixing length parameter to be $α = 0.90{^{+0.36}_{-0.26}}$ (when assumed to be constant across our sample), and obtain constraints on the parameters describing turbulent dissipation and convective flux.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Frozen-Flux Modelling for Epochs 1915 and 1980

The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.

OBrien, Michael S.↗
DOE OSTI2022

Sensitivity of Cloud-Radiative Effects to Cloud Fraction Parametrizations in Tropical, Mid-Latitude and Arctic Kilometre-Scale Simulations

The Regional Atmosphere (RA) configuration of the Met Office Unified Model currently requires different cloud fraction parametrizations (CFP) for tropical and midlatitude simulations. To explore the scope for unification of these two RA configurations, this paper presents a detailed evaluation of simulations over tropical, midlatitude and arctic domains, with two different diagnostic CFPs, a prognostic CFP, and no CFP at all. Furthermore, a novel, hybrid approach was used that treats liquid cloud diagnostically and ice cloud prognostically. Using observations from three U.S. Department of Energy Atmospheric Radiation Measurement supersites, it is shown that none of these CFPs stands out as superior over all domains. Over the frequently overcast Arctic, the all-or-nothing approach best captures the cloud radiative properties. Conversely, CFPs are of benefit in regions with frequent partial cloudiness, such as the midlatitudes and the tropics. However, their improved cloud radiative properties often hide an error compensation. All models underestimate overcast, low-base cloud with small water paths in convective environments. In addition, mid-latitude overcast, low-base, optically thick clouds in the morning, possibly associated with overnight convection, are frequently too broken. Diagnostic schemes compensate for these errors by producing spurious, scattered afternoon cloud, which could be due to a correct cloud response to too eager convective initiation. Winter clouds over the midlatitudes are improved when liquid cloud is represented diagnostically with a bimodal saturation-departure PDF, without error compensation. While it is difficult to unify the RA across the globe around a single CFP scheme, the newly proposed Hybrid scheme performs reasonably well for cloud cover across all regions. It also exhibits SW biases that are smaller than most other configurations and is less affected by excessive liquid water paths and compensating errors than fully diagnostic schemes. Surface precipitation is fairly insensitive to the CFP in the simulations shown here.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Optimized Profile Retrievals of Aerosol Microphysical Properties from Simulated Spaceborne Multiwavelength Lidar

This work is an expanded study of one previously published onretrievals of aerosol microphysical properties from space-borne multiwavelengthlidarmeasurements. The earlier studiesand this one weredone in the framework of the NASA Aerosol-Clouds-Ecosystems (now the Aerosol Clouds Convection and Precipitation) NASA mission. The focus here is on the capabilities of a simulated spacebornemultiwavelengthlidar system for retrieving aerosol complex refractive index (m = mr+ imi) and spectral single scattering albedo (SSA(λ)), although other bulk parameters such as effective (reff) radius and particle volume (V) and surface (S) concentrations are also studied. The novelty presented here is the use of recently published, case dependent optimized-constraints on the microphysical retrievals using three backscattering coefficients (β) at 355, 532 and 1064 nm and two extinction coefficients (α) at 355 and 532 nm, typically known as the stand-alone 3β+2α lidar inversion. Case-dependent optimized-constraints (CDOC) limit the ranges of refractive index, both real (mr) and imaginary (mi) parts, and of radii that are permitted in the retrievals. Such constraints are selected directly from the 3β+2α41measurements through an analysis of the relationship between spectral dependence of aerosol extinction-to-backscatter ratios (LR) and the Ångström exponent of extinction. The analyses presented here for different sets ofsize distributions and refractive indices reveal that the direct determination of CDOCareonly feasible for cases where the uncertaintiesin the input optical data areless than 15 %.Forthe same simulated spacebornesystem and yield than in Whiteman et al., (2018), we demonstrated that the use of CDOC as essential for the retrievals of refractive index and also largely improved retrieval of bulk parameters. A discussion of the global representativeness of CDOC is presented using simulated lidar data from a 24-hour satellite track using GEOS model output to initialize the lidar simulator.We found that CDOCare representative of many aerosol mixtures in spite of some outliers (e.g. highly hydrated particles) associatedwith the assumptions of bimodal size distributions and of the same refractive index for fine and coarse modes. Moreover, sensitivity tests performed using synthetic data reveal that retrievals of imaginary refractive index (mi) and SSA are extremely sensitive to β(355).

NASA Aerosol-Clouds-Ecosystems↗
DOE OSTI2025

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗