Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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 1,009 records · Page 56

GCM Simulations of Cirrus Clouds and Cloud Feedbacks

Cirrus clouds are a particularly uncertain component of general circulation model (GCM simulations of long-term climate change for a variety of reasons: (1) They encompass a wide range of optical thicknesses and altitudes, from thin tropopause cirrus to thick anvil cirrus that descend to the freezing level, and thus can exert both positive and negative forcing and feedback on the climate; (2) The dynamical processes that create them are poorly resolved in climate GCMs and different in the tropics and midlatitudes; (3) Predictions of their formation and microphysical properties depend on the accuracy of dynamical transports of small concentrations of water vapor to and within the upper troposphere; (4) The relative humidity conditions at which they form depends on the nature and concentration of nucleating particles and is poorly understood; (5) They are more difficult to observe than other cloud types, and hence their parameterization is more loosely constrained by available data. We will illustrate the potential sensitivity of the perturbed climate to uncertainties in cirrus cloud formulation. We will also examine the processes that form cirrus in climate models and discuss the accuracy with which climate GCMs represent these processes. We will also discuss ways in which GCM grid-scale parameterizations might be derived from cloud-scale observations. Finally, we will emphasize the types of global observations needed to constrain parameterizations of cirrus in climate GCMs.

DelGenio, Anthony D.↗

Measurement report: Closure Analysis of Aerosol–cloud Composition in Tropical Maritime Warm Convection

Cloud droplet chemical composition is a key observable property that can aid understanding of how aerosols and clouds interact. As part of the Clouds, Aerosols and Monsoon Processes – Philippines Experiment (CAMP2Ex), three case studies were analyzed involving collocated airborne sampling of relevant clear and cloudy air masses associated with maritime warm convection. Two of the cases represented a polluted marine background, with signatures of transported East Asian regional pollution, aged over water for several days, while the third case comprised a major smoke transport event from Kalimantan fires. Sea salt was a dominant component of cloud droplet composition, in spite of fine particulate enhancement from regional anthropogenic sources. Furthermore, the proportion of sea salt was enhanced relative to sulfate in rainwater and may indicate both a propensity for sea salt to aid warm rain production and an increased collection efficiency of large sea salt particles by rain in subsaturated environments. Amongst cases, as precipitation became more significant, so too did the variability in the sea salt to (non-sea salt) sulfate ratio. Across cases, nitrate and ammonium were fractionally greater in cloud water than fine-mode aerosol particles; however, a strong covariability in cloud water nitrate and sea salt was suggestive of prior uptake of nitrate on large salt particles. A mass-based closure analysis of non-sea salt sulfate compared the cloud water air-equivalent mass concentration to the concentration of aerosol particles serving as cloud condensation nuclei for droplet activation. While sulfate found in cloud was generally constrained by the sub-cloud aerosol concentration, there was significant intra-cloud variability that was attributed to entrainment – causing evaporation of sulfate-containing droplets –and losses due to precipitation. In addition, precipitation tended to promote mesoscale variability in the sub-cloud aerosol through a combination of removal, convective downdrafts, and dynamically driven convergence. Physical mechanisms exerted such strong control over the cloud water compositional budget that it was not possible to isolate any signature of chemical production/loss using in-cloud observations. The cloud-free environment surrounding the non-precipitating smoke case indicated sulfate enhancement compared to convective mixing quantified by a stable gas tracer; however, this was not observed in the cloud water (either through use of ratios or the mass closure), perhaps implying that the warm convective cloud timescale was too short for chemical production to be a leading-order budgetary term and because precursors had already been predominantly exhausted. Closure of other species was truncated by incomplete characterization of coarse aerosol (e.g., it was found that only 10 %–50% of sea salt mass found in cloud was captured during clear-air sampling) and unmeasured gasphase abundances affecting closure of semi-volatile aerosol species (e.g., ammonium, nitrate and organic) and soluble volatile organic compound contributions to total organic carbon in cloud water.

Ewan Crosbie↗

Analog Environments for Venus Aerosol Instrument and Mission Concept Testing

The most common hypothesis discussed for life on Venus is a habitat in the dense, persistent clouds, which are primarily sulfuric acid with a secondary water component. Although Earth lacks a complete Venus analog environment, several partial analogs exist which can inform our understanding of the requirements for cloud habitability as well as serving as potential testbeds for future Venus instruments or missions. Earth’s tropospheric clouds have larger particles, more water, and shorter lifetimes than Venus’s; however, they have a partial overlap in pressure and temperature, and are our only example of an airborne habitat.They can inform how particle residence time and microbial generation time constrain habitable atmospheric regions. On the biological side, this includes available nutrients and energy, environmental stressors such as radiation, osmotic potential and acidity, and time spent in protective but inactive forms (endospores, cysts, akinetes); on the physical side, gravity, air density and viscosity, thermal lofting, gravity waves, scavenging due to precipitation, and other airflow dynamics, as well as effective particle radius and electrical charge. They can also enable sampling strategy tests for airborne life detection, ranging from10 2 to10 8 viable cells/mL (10 -3 to 10 -9 cells/particle). Conversely, Earth's stratospheric sulfate layer lacks an active microbial presence, and is at a colder and less dense altitude but has partially analogous values in terms of H2SO4 and water concentration, particle size, and number density. These aerosols are accessible to aircraft and high-altitude balloons. Even more accessible potential analogs include sea spray, which overlaps with the smaller (haze) particle sizes and densities at Venus and contains trace organic constituents; marine fog, which at higher densities can match Venus’s larger (Mode 2) particle properties and contains a range of trace compounds including ammonia, phosphates, sulfates, and nitrates; and spray from acidic hydrothermal systems, which can cover a wide variety of sizes, densities, and chemical compositions. These analogs could improve the development of Venus cloud sampling instruments or sondes, including field tests of material compatibility, particle capture efficiency, and analysis sensitivity.

Analog↗

A machine learning approach targeting parameter estimation for plant functional type coexistence modeling using ELM-FATES (v2.0)

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

54 ENVIRONMENTAL SCIENCES↗

Conceptual design of highly-constrained splitters for the FFA@CEBAF energy upgrade study

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab is investigating a significant energy upgrade utilizing Fixed-Field Alternating-gradient (FFA) recirculating arcs. This upgrade requires the design of complex horizontal beam splitters to manage up to six concurrent beam passes. This paper presents the conceptual design of these splitters, which are subject to severe physical constraints imposed by the existing accelerator tunnel and multifaceted beam dynamics requirements for matching into the permanent-magnet FFA arcs. The design methodology, centered on multi-pass simulations in the Bmad toolkit, is detailed from the initial geometric layout through the advanced optics matching. Key results include a robust geometric arrangement that fits within the spatial boundaries and the development of multiple, flexible optics matching solutions. Furthermore, the design integrates a viable scheme for extracting high-energy beams for the experimental halls, a critical operational requirement. This work establishes a comprehensive and viable conceptual design, forming a baseline for future engineering and performance optimization studies.

Bodenstein, R.M. [Thomas Jefferson National Accele↗

Spin and lattice dynamics of the two-dimensional van der Waals ferromagnet CrI3

Chromium tri-iodide (CrI 3 ) is a prototypical ferromagnetic van der Waals insulator with its genuinely two-dimensional (2D) long-range magnetic order below 45 K demonstrated recently. The underlying magnetic anisotropy has not been completely understood while both the Dzyaloshinskii-Moriya (DM) interaction and the Kitaev—Γ type interaction have been considered as the relevant magnetic Hamiltonian. In addition, the relation between the crystal structure and the magnetic order needs to be further elucidated concerning their possible coupling in few-layer samples and in the topmost surface layers of bulk samples. Here, we investigate these issues via temperature- and magnetic field-dependent terahertz spectroscopy on bulk CrI 3 single crystals, focusing on the dynamics of ferromagnetic resonances (FMRs) and optical phonons in the terahertz (THz) region from 4 to 120 cm -1 (from 0.5 to 15 meV). We narrow down the possible ranges of the interaction parameters such as the off-diagonal symmetric terms and the single-ion anisotropy. The accurate values of these parameters significantly constrain the magnitude of possible Kitaev—Γ exchange interaction and the topological magnon gap. Moreover, the structure-magnetism relationship was critically analyzed based on the temperature- and field-dependences of two E u in-plane optical phonon modes, which shows that the commonly believed structural phase diagram of CrI 3 , derived from surface-preferential data, has to be seriously modified.

36 MATERIALS SCIENCE↗

A Stochastic Fractional Dynamics Model of Space-time Variability of Rain

Rainfall varies in space and time in a highly irregular manner and is described naturally in terms of a stochastic process. A characteristic feature of rainfall statistics is that they depend strongly on the space-time scales over which rain data are averaged. A spectral model of precipitation has been developed based on a stochastic differential equation of fractional order for the point rain rate, that allows a concise description of the second moment statistics of rain at any prescribed space-time averaging scale. The model is thus capable of providing a unified description of the statistics of both radar and rain gauge data. The underlying dynamical equation can be expressed in terms of space-time derivatives of fractional orders that are adjusted together with other model parameters to fit the data. The form of the resulting spectrum gives the model adequate flexibility to capture the subtle interplay between the spatial and temporal scales of variability of rain but strongly constrains the predicted statistical behavior as a function of the averaging length and times scales. We test the model with radar and gauge data collected contemporaneously at the NASA TRMM ground validation sites located near Melbourne, Florida and in Kwajalein Atoll, Marshall Islands in the tropical Pacific. We estimate the parameters by tuning them to the second moment statistics of radar data. The model predictions are then found to fit the second moment statistics of the gauge data reasonably well without any further adjustment.

TRMM↗

On the Reconstruction of Palaeo-Ice Sheets: Recent Advances and Future Challenges

Reconstructing the growth and decay of palaeo-ice sheets is critical to understanding mechanisms of global climate change and associated sea-level fluctuations in the past, present and future. The significance of palaeo-ice sheets is further underlined by the broad range of disciplines concerned with reconstructing their behaviour, many of which have undergone a rapid expansion since the 1980s. In particular, there has been a major increase in the size and qualitative diversity of empirical data used to reconstruct and date ice sheets, and major improvements in our ability to simulate their dynamics in numerical ice sheet models. These developments have made it increasingly necessary to forge interdisciplinary links between sub-disciplines and to link numerical modelling with observations and dating of proxy records. The aim of this paper is to evaluate recent developments in the methods used to reconstruct ice sheets and outline some key challenges that remain, with an emphasis on how future work might integrate terrestrial and marine evidence together with numerical modelling. Our focus is on pan-ice sheet reconstructions of the last deglaciation, but regional case studies are used to illustrate methodological achievements, challenges and opportunities. Whilst various disciplines have made important progress in our understanding of ice-sheet dynamics, it is clear that data-model integration remains under-used, and that uncertainties remain poorly quantified in both empirically-based and numerical ice-sheet reconstructions. The representation of past climate will continue to be the largest source of uncertainty for numerical modelling. As such, palaeo-observations are critical to constrain and validate modelling. State-of-the-art numerical models will continue to improve both in model resolution and in the breadth of inclusion of relevant processes, thereby enabling more accurate and more direct comparison with the increasing range of palaeo-observations. Thus, the capability is developing to use all relevant palaeo-records to more strongly constrain deglacial (and to a lesser extent pre-LGM) ice sheet evolution. In working towards that goal, the accurate representation of uncertainties is required for both constraint data and model outputs. Close cooperation between modelling and data-gathering communities is essential to ensure this capability is realised and continues to progress.

Ice sheet reconstruction↗

Reduced Complexity Model Intercomparison Project Phase 2: Synthesizing Earth System Knowledge for Probabilistic Climate Projections

Over the last decades, climate science has evolved rapidly across multiple expert domains. Our best tools to capture state-of-the-art knowledge in an internally self-consistent modelling framework are the increasingly complex fully coupled Earth System Models (ESMs). However, computational limitations and the structural rigidity of ESMs mean that the full range of uncertainties across multiple domains are difficult to capture with ESMs alone. The tools of choice are instead more computationally efficient reduced complexity models (RCMs), which are structurally flexible and can span the response dynamics across a range of domain-specific models and ESM experiments. Here we present Phase 2 of the Reduced Complexity Model Intercomparison Project (RCMIP Phase 2), the first comprehensive intercomparison of RCMs that are probabilistically calibrated with key benchmark ranges from specialized research communities. Unsurprisingly, but crucially, we find that models which have been constrained to reflect the key benchmarks better reflect the key benchmarks. Under the low-emissions SSP1-1.9 scenario, across the RCMs, median peak warming projections range from 1.3 to 1.7°C (relative to 1850-1900, using an observationally-based historical warming estimate of 0.8°C between 1850-1900 and 1995-2014). Further developing methodologies to constrain these projection uncertainties seems paramount given the international community's goal to contain warming to below 1.5°C above pre-industrial in the long-term. Our findings suggest that users of RCMs should carefully evaluate their RCM, specifically its skill against key benchmarks and consider the need to include projections benchmarks either from ESM results or other assessments to reduce divergence in future projections.

Climate↗

Unravelling the orbits of cluster galaxy populations according to their dominant gas ionization source

ABSTRACT We investigate the kinematical and dynamical properties of cluster galaxy populations classified according to their dominant source of gas ionization, namely: star-forming (SF) galaxies, optical active galactic nuclei (AGNs), mixed SF plus AGN ionization (transition objects, T), and quiescent (Q) galaxies. We stack 8892 member galaxies from 336 relaxed galaxy clusters to build an ensemble cluster and estimate the observed projected profiles of numerical density and velocity dispersion, $\sigma _P(R)$, of each galaxy population. The MAMPOSSt code and the Jeans equations inversion technique are used to constrain the velocity anisotropy profiles of the galaxy populations in both parametric and non-parametric ways. We find that Q (SF) galaxies display the lowest (highest) typical cluster-centric distances and velocity dispersion values. Transition galaxies are more concentrated and tend to exhibit lower velocity dispersion values than SF galaxies. Galaxies that host an optical AGN are as concentrated as Q galaxies but display velocity dispersion values similar to those of the SF population. MAMPOSSt is able to find equilibrium solutions that successfully recover the observed $\sigma _P(R)$ profile only for the Q, T, and AGN populations. We find that the orbits of all populations are consistent with isotropy in the inner regions, becoming increasingly radial with the distance from the cluster centre. These results suggest that Q galaxies are in equilibrium within their clusters, while SF galaxies have more recently arrived in the cluster environment. Finally, the T and AGN populations appear to be in an intermediate dynamical state between those of the SF and Q populations.

Valk, Greique A. (ORCID:0009000827731299)↗

Constraining Primordial Black Holes with Dwarf Galaxy Heating

Black holes formed in the early universe, prior to the formation of stars, can exist as dark matter and also contribute to the black hole merger events observed in gravitational waves. We set a new limit on the abundance of primordial black holes (PBHs) by considering interactions of PBHs with the interstellar medium, which result in the heating of gas. We examine generic heating mechanisms, including emission from the accretion disk, dynamical friction, and disk outflows. Using the data from the Leo T dwarf galaxy, we set a new cosmology-independent limit on the abundance of PBHs in the mass range ${ \mathcal O }(1){M}_{\odot }\mbox{-}{10}^{7}{M}_{\odot }$, relevant for the recently detected gravitational-wave signals from intermediate-mass BHs.

79 ASTRONOMY AND ASTROPHYSICS↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

CryoFold: Determining protein structures and data-guided ensembles from cryo-EM density maps

Cryoelectron microscopy requires molecular modeling for refinement of structures. Ensemble models arrive at low free-energy molecular structures, but are computationally expensive and limited to resolving only small proteins. Here, we introduce CryoFold, a pipeline of molecular dynamics simulations that determines ensembles of protein structures by integrating density data of varying sparsity at 3–5 Å resolution with sequence information and coarse-grained topological knowledge of the protein folds. We present six examples, folding proteins between 72 and 2,000 residues, including large membrane and multi-domain systems, and results from two Electron Microscopy Data Bank (EMDB) competitions. Driven by data from a single state, CryoFold discovers ensembles of common low-energy models together with rare low-probability structures that capture the equilibrium distribution of proteins constrained by the density maps. Many of these conformations are experimentally validated and functionally relevant. We arrive at a set of best practices for data-guided protein folding that are controlled using a Python graphical user interface (GUI).

59 BASIC BIOLOGICAL SCIENCES↗

The Open-Ocean Sensible Heat Flux and Its Significance for Arctic Boundary Layer Mixing During Early Fall

The increasing ice-free area during late summer has transformed the Arctic to a climate system with more dynamic boundary layer (BL) clouds and seasonal sea ice growth. The open-ocean sensible heat flux, a crucial mechanism of excessive ocean heat loss to the atmosphere during the fall freeze season, is speculated to play an important role in the recently observed cloud cover increase and BL instability. However, lack of observations and understanding of the resilience of the proposed mechanisms, especially in relation to meteorological and interannual variability, has left a poorly constrained BL parameterization scheme in Arctic climate models. In this study, we use multiyear Japanese cruise-ship observations from RV Mirai over the open Arctic Ocean to characterize the surface sensible heat flux (SSHF) during early fall and investigate its contribution to BL turbulence. It is found that mixing by SSHF is favored during episodes of high surface wind speed and is also influenced by the prevailing cloud regime. The deepest BLs and maximum ocean-atmosphere temperature difference are observed during cold air advection (associated with the stratocumulus regime), yet, contrary to previous speculation, the efficiency of sensible heat exchange is low. On the other hand, the SSHF contributes significantly to BL mixing during the uplift (low pressure) followed by the highly stable (stratus) regime. Overall, it can explain 10 of the open ocean BL height variability, whereas cloud-driven (moisture and radiative) mechanisms appear to be the other dominant source of convective turbulence. Nevertheless, there is strong interannual variability in the relationship between the SSHF and the BL height which can be intensified by the changing occurrence of Arctic climate patterns, such as positive surface wind speed anomalies and more frequent conditions of uplift. This study highlights the need for comprehensive BL observations like the RV Mirai for better understanding and predicting the dynamic nature of the Arctic climate.

boundary layer↗

Cation Disorder of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ in Super‐Earth Mantles

Understanding the mineralogy of exoplanets is essential for unraveling their interior structures, dynamics, and evolution. For large super-Earths, the post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$, one of the major mantle phases, may undergo the order-disorder transition (ODT) at high temperatures. However, the ODT phase boundary of ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ has not been rigorously constrained. Additionally, fundamental thermodynamic properties of the disordered ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ remain poorly investigated. Here, we develop a unified machine learning potential (MLP) for ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ of ab initio accuracy under super-Earth mantle conditions. With the efficient MLP, we extensively calculate the free energy of post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ via the thermodynamic integration method. The results are used to constrain the ODT phase boundary. Furthermore, we report the P-V-T equation of state and Grüneisen parameters for post-post spinel ${\text{Mg}}_{\mathbf{2}}{\text{SiO}}_{\mathbf{4}}$ across various degrees of disorder. These thermodynamic properties are further applied to update the adiabatic thermal profiles and the mass-radius relation of super-Earths.

36 MATERIALS SCIENCE↗

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Final Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aerodynamic design optimization using sensitivity analysis and computational fluid dynamics

A new and efficient method is presented for aerodynamic design optimization, which is based on a computational fluid dynamics (CFD)-sensitivity analysis algorithm. The method is applied to design a scramjet-afterbody configuration for an optimized axial thrust. The Euler equations are solved for the inviscid analysis of the flow, which in turn provides the objective function and the constraints. The CFD analysis is then coupled with the optimization procedure that uses a constrained minimization method. The sensitivity coefficients, i.e. gradients of the objective function and the constraints, needed for the optimization are obtained using a quasi-analytical method rather than the traditional brute force method of finite difference approximations. During the one-dimensional search of the optimization procedure, an approximate flow analysis (predicted flow) based on a first-order Taylor series expansion is used to reduce the computational cost. Finally, the sensitivity of the optimum objective function to various design parameters, which are kept constant during the optimization, is computed to predict new optimum solutions. The flow analysis of the demonstrative example are compared with the experimental data. It is shown that the method is more efficient than the traditional methods.

Baysal, Oktay↗