Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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,603 records · Page 89

Effect of Stability on Mixing in Open Canopies

In open canopies, the within-canopy flux from the ground surface and understory can account for a significant fraction of the total flux above the canopy. This study incorporates the important influence of within-canopy stability on turbulent mixing and subcanopy fluxes into a first-order closure scheme. Toward this goal, we analyze within-canopy eddy-correlation data from the old aspen site in the Boreal Ecosystem - Atmosphere Study (BOREAS) and a mature ponderosa pine site in Central Oregon, USA. A formulation of within-canopy transport is framed in terms of a stability- dependent mixing length, which approaches Monin-Obukhov similarity theory above the canopy roughness sublayer. The new simple formulation is an improvement upon the usual neglect of the influence of within-canopy stability in simple models. However, frequent well-defined cold air drainage within the pine subcanopy inversion reduces the utility of simple models for nocturnal transport. Other shortcomings of the formulation are discussed.

Lee, Young-Hee↗

Magnetic state switching in FeGa microstructures

This work demonstrates that magnetoelectric composite heterostructures can be designed at the length scale of 10µms that can be switched from a magnetized state to a vortex state, effectively switching the magnetization off, using electric field induced strain. This was accomplished using thin film magnetoelectric heterostructures of Fe 81.4 Ga 18.6 on a single crystal (011) [Pb(Mg 1/3 Nb 2/3 )O 3 ] 0.68 -[PbTiO 3 ] 0.32 (PMN-32PT) ferroelectric substrate. Here, the heterostructures were tripped from a multi-domain magnetized state to a flux closure vortex state using voltage induced strain in a piezoelectric substrate. FeGa heterostructures were deposited on a Si-substrate for superconducting quantum interference device magnetometry characterization of the magnetic properties. The magnetoelectric coupling of a FeGa continuous film on PMN-32PT was characterized using a magneto optical Kerr effect magnetometer with bi-axial strain gauges, and magnetic multi-domain heterostructures were imaged using x-ray magnetic circular dichroism—photoemission electron microscopy during the transition to the vortex state. The domain structures were modelled using MuMax 3 , a micromagnetics code, and compared with observations. The results provide considerable insight into designing magnetoelectric heterostructures that can be switched from an ‘on’ state to an ‘off’ state using electric field induced strain.

36 MATERIALS SCIENCE↗

The Sealing Mechanisms of a Fracture in Opalinus Clay as Revealed by in situ Synchrotron X-Ray Micro-Tomography

The detailed mechanisms of the sealing of a single fracture, from hydration to almost complete closure by increase of confining pressure, as monitored from in situ synchrotron X-ray microtomography during the flow of carbonated water, is here shown for the first time. Different mechanisms play the key role at different stages in the evolution of the fracture. Hydration mechanically weakens the surfaces of the fracture and induces a first closure due to microcracking at the asperity contacts, increasing their size and creating choke points. Increase in confining stress promptly hydraulically seals the fracture by closing the main choke point, with a relative small deformation of the sample. Finally, the more pervasive mechanical deformation observed at higher stresses almost completely seals the whole fracture. The evolution of the sample has been described and quantified using 4D image processing, focusing on the evolution of aperture and digital volume correlation. Hydraulic properties of the sample at different stages have been modeled via Stokes flow simulation, and the results compared to the morphometric analysis, finding positive correlations with the average fracture aperture variation along the flowpath in function of confining pressure. Opalinus Clay is found to be a rock markedly prone to sealing in case of flow with carbonated water; this behavior is the result of the large fraction of clays and of its microstructure, lacking both cementing phases and large stiff particles. CO2 in this sample has no evident role in the evolution of the fracture; chemically induced weathering on the surface has not been detected, in contrast with the behavior observed in samples with carbonates as cementing phase.

58 GEOSCIENCES↗

Auroral Current and Electrodynamics Structure (ACES) Observations of Ionospheric Feedback in the Alfven Resonator

In 2009, the Auroral Current and Electrodynamics Structure (ACES) High and Low sounding rockets were launched from the Poker Flat Rocket Range (PFRR) in Alaska, with the science objective of gathering in-situ data to quantify current closure in a discrete auroral arc. As ACES High crossed through the return current of an arc (that was monitored using an all sky camera from the ground at Fort Yukon), its instruments recorded clear Alfv nic signatures both poleward and equatorward of the return current region, but not within the main region of the return current itself. These data provide an excellent opportunity to study ionospheric feedback and how it interacts with the Alfv n resonator. We compare the observations with predictions and new results from a model of ionospheric feedback in the ionospheric Alfv n resonator (IAR) and report the significance and impact of these new data for the Magnetosphere-Ionosphere Coupling in the Alfv n Resonator (MICA) rocket mission to launch from PFRR this winter. MICA s primary science objectives specifically focus on better understanding the small-scale structure that the model predicts should exist within the return current region.

Cohen, Ian J.↗

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A kinetic approach to compressible turbulence

A kinetic theory of strong turbulence is developed by decomposing fluctuation of distribution function into scaled components, or ranks, to represent three transport processes: macroscopic evolution, transport property and relaxation. Closure is obtained by formulating a relaxation process for the approach of eddy diffusivity to equilibrium. The kinetic equation is solved and three transport functions of turbulence in adiabatic gas are derived: two transfer functions describe turbulent mode transfer, and a coupling function describes the build-up of kinetic energy at the expense of internal energy. The energy balance in the inertial subrange of homogeneous and isotropic turbulence in an adiabatic gas yields a k to the -3rd law, in agreement with numerical computations from one dimensional model of compressible turbulence.

Tchen, C. M.↗

Mean climate and transience in the tropics of the UGAMP GCM: Sensitivity to convective parametrization

The sensitivity of the UK Universities Global Atmospheric Modelling Programme (UGAMP) General Circulation Model (UGCM) to two very different approaches to convective parametrization is described. Comparison is made between a Kuo scheme, which is constrained by large-scale moisture convergence, and a convective-adjustment scheme, which relaxes to observed thermodynamic states. Results from 360-day integrations with perpetual January conditions are used to describe the model's tropical time-mean climate and its variability. Both convection schemes give reasonable simulations of the time-mean climate, but the representation of the main modes of tropical variability is markedly different. The Kuo scheme has much weaker variance, confined to synoptic frequencies near 4 days, and a poor simulation of intraseasonal variability. In contrast, the convective-adjustment scheme has much more transient activity at all time-scales. The various aspects of the two schemes which might explain this difference are discussed. The particular closure on moisture convergence used in this version of the Kuo scheme is identified as being inappropriate.

Slingo, Julia↗

Assessing Climate Change Impacts on the Stability of Small Tidal Inlets: Part 2- Data Rich Environments

Climate change (CC) is likely to affect the thousands of bar-built or barrier estuaries (here referred to as Small tidal inlets - STIs) around the world. Any such CC impacts on the stability of STIs, which governs the dynamics of STIs as well as that of the inlet-adjacent coastline, can result in significant socio-economic consequences due to the heavy human utilisation of these systems and their surrounds. This article demonstrates the application of a process based snap-shot modelling approach, using the coastal morphodynamic model Delft3D, to 3 case study sites representing the 3 main STI types; Permanently open, locationally stable inlets (Type 1), Permanently open, alongshore migrating inlets (Type 2) and Seasonally/Intermittently open, locationally stable inlets (Type 3). The 3 case study sites (Negombo lagoon - Type 1, Kalutara lagoon - Type 2, and Maha Oya river - Type 3) are all located along the southwest coast of Sri Lanka. After successful hydrodynamic and morphodynamic model validation at the 3 case study sites, CC impact assessment are undertaken for a high end greenhouse gas emission scenario. Future CC modified wave and riverflow conditions are derived from a regional scale application of spectral wave models (WaveWatch III and SWAN) and catchment scale applications of a hydrologic model (CLSM) respectively, both of which are forced with IPCC Global Climate Model output dynamically downscaled to approximately 50 km resolution over the study area with the stretched grid Conformal Cubic Atmospheric Model CCAM. Results show that while all 3 case study STIs will experience significant CC driven variations in their level of stability, none of them will change Type by the year 2100. Specifically, the level of stability of the Type 1 inlet will decrease from 'Good' to 'Fair to poor' by 2100, while the level of (locational) stability of the Type 2 inlet will also decrease with a doubling of the annual migration distance. Conversely, the stability of the Type 3 inlet will increase, with the time till inlet closure increasing by approximately 75%. The main contributor to the overall CC effect on the stability of all 3 STIs is CC driven variations in wave conditions and resulting changes in longshore sediment transport, not Sea level rise as commonly believed.

IPC↗

Predicting tungsten erosion and leakage properties for the new V-shaped small angle slot divertor in DIII-D

Impurity transport modeling of the new tungsten (W)-coated, V-shaped small angle slot (SAS) divertor in the DIII-D tokamak was conducted using the SOLPS-ITER plasma edge code package and the DIVIMP impurity tracking code. The inboard baffle of the current SAS divertor will be shifted closer to the outboard baffle, creating a V-corner at the slot vertex. In addition, the outboard baffle will be coated with 10–15 μm of W for experiments studying high-Z sourcing and leakage in a closed divertor. Modeling of the 'SAS-VW' divertor predicts that these changes to the inner baffle will reduce W gross erosion by 40× relative to the existing SAS divertor when the outer strike point (OSP) is at the V-corner and the ion B × ∇B drift is towards the divertor, driven primarily by significant cooling near the slot vertex. Most W erosion in SAS-VW is expected to occur near the slot entrance, which may pose a higher risk to core contamination than W eroded deeper in the slot. Furthermore, adding a new sheath-based prompt redeposition model increases the sensitivity of redeposition estimates to near-target plasma conditions and may provide more accurate predictions of net erosion. Moving the OSP outboard from the slot vertex ~4 cm onto the W-coated region yields a 40× increase in the gross erosion rate and a 50% decrease in the core leakage fraction. Thus slight variations in strike point location may counteract the potential benefits of the tightly-baffled V slot on minimizing erosion. This impurity transport modeling provides useful guidance for future experiments on the SAS-VW divertor focused on high-Z erosion/redeposition, scrape-off layer transport, and core leakage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic response of a freely rotating butterfly valve in the advanced test reactor – dynamic fluid-body interaction modeling

To regulate primary coolant flow in the Advanced Test Reactor (ATR), a butterfly valve was installed between the primary coolant pumps and the reactor core. If the mechanical connection between the valve's disk and its shaft ever fails, the disk will rotate freely. Rapid disk rotation may induce pressure surges that could damage upstream pipes. In the present work, the rotational trajectory and pressure evolution during a disk free-rotation scenario were analyzed via the dynamic fluid-body interaction (DFBI) approach in STAR-CCM+, with the movement of a solid being driven by the forces and moment/torque imposed by its surrounding fluid. Assuming a large initial opening angle, the disk accelerates slowly, but swiftly passes the closed position. As a result of the sudden valve closure, a pressure surge occurs in the upstream pipes, exceeding their maximum allowable pressure. Furthermore, the disk does not stabilize at the closed position but continually oscillates due to the unsteady nature of the coolant flow. Because of the significant and continuous water hammer effect, a fix to the butterfly valve is being implemented to prevent rapid valve closure due to failure at the valve's disk/shaft mechanical connection.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Feasibility and the Benefits of the Advanced Nuclear Fuel Pellet Designs with Radially Varying Fuel Zoning and Burnable Poison Concentration (Final Summary Report)

As part of the work supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE) Gateway for Accelerated Innovation in Nuclear (GAIN) FY 2021 Voucher, Exelon Generation, now Constellation, and Oak Ridge National Laboratory (ORNL) entered into a cooperative research and development agreement (CRADA) to evaluate and assess the feasibility and impact of various conceptual advanced nuclear fuel pellet designs (ANFPDs). The objective of this project was to perform modeling and simulation and analyses using the advanced modeling and simulation capabilities of VERA/BISON, developed by DOE, to determine the viability and benefits of numerous advanced nuclear fuel pellet design concepts in terms of fuel cycle costs, operational safety, and margin improvement. Whereas the detailed coupled neutronic and thermal hydraulic analyses performed using VERA provided in-depth knowledge in terms of fuel cycle performance, the detailed VERA results were used in subsequent fuel performance analyses using BISON. These subsequent analyses focused on several key fuel performance criteria, such as peak fuel centerline temperature (FCT), fission gas release (FGR), gap closure, plenum pressure, and cladding hoop stress. These key fuel performance criteria were analyzed for some of the conceptual fuel designs and were compared with the results obtained for UO 2 fuel. The results provided detailed information to enable better understanding of the performance of the fuel types analyzed. Understanding the advantage of loading these conceptual fuel designs into the core is important not only to Constellation but also to the entire light-water reactor (LWR) fleet in the United States.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Theoretical analysis of aircraft afterbody flow

The strongly interactive flow field about aircraft afterbodies is investigated using computational techniques by which the thin-shear-layer formulation of the compressible, Reynolds-averaged Navier-Stokes equations is solved. A time-dependent implicit numerical algorithm is used to obtain solutions for a variety of afterbody and nozzle geometries, within the class of bodies of revolution, for both subsonic and supersonic external flow, and for sonic and supersonic underexpanded jets. Only centered nozzles with either a sharp lip or a blunt base are considered. In all cases, computed results are compared with experimental data taken at flight Reynolds numbers for like-flow conditions. Turbulence closure is realized using algebraic eddy-viscosity concepts. A new and unique adaptive-grid technique is used to resolve flow regimes with large gradients and to improve the accuracy and efficiency of the computational scheme. Special singular point boundary conditions are used for similar purposes, and are especially effective for highly under-expanded jets. For all cases considered, except one with a very large base-to-nozzle-exit-diameter ratio, the agreement with experimental measurements is excellent. For geometries with large base regions, enhancements in the turbulence transport model are necessary to support improvements in the flow-field simulation.

Deiwert, G. S.↗

Simulation of Multistage Turbine Flows

The numerical simulation of turbine flows serves to enhance the understanding of the flow phenomena within multistage turbomachinery components. The direct benefit of this activity is improved modeling capability, which can be used to improve component efficiency and durability. A hierarchy of equations was formulated to assess the difficulty in analyzing the flow field within multistage turbomachinery components. The Navier-Stokes equations provides the most complete description. The simplest description is given by a set of equations that govern the quasi-one-dimensional flow. The number of unknowns to be solved for increases monotonically above the number of equations. The development of the additional set of equations needed to mathematically close the system of equations forms the closure problem associated with that level of description. For the Navier-Stokes equation there is no closure problem. For the quasi-one-dimensional equation set random flow fluctuations, unsteady fluctuations, nonaxisymmetric flow variations, and hub-to-shroud variations on the quasi-one-dimensional flow must be accounted for.

Celestina, M. L.↗

Dependence of CCN closure relationship with organic fraction from two airborne field campaigns over mid-latitude land and ocean

This study investigates the relationship between measured and calculated cloud condensation nuclei (CCN) number concentration and its dependence with organic fraction utilizing aircraft observations from The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA, 2017–2018) and The Holistic Interactions of Shallow Clouds, Aerosols, and Land Ecosystems (HI-SCALE, 2016) campaigns, which represent midlatitude marine and continental environments, respectively. For the ACE-ENA marine region, aerosol and CCN concentrations were significantly higher in summer than in winter, whereas at continental site for HI-SCALE, aerosol and CCN concentrations showed no pronounced differences between spring and autumn. Using aerosol chemical composition and number size distribution data, CCN concentrations at various supersaturations are calculated based on Köhler theory and then compared with observations from CCN counter. The results show that CCN closure performs well at both sites with a slight overestimation, with mean closure ratio (CR) of 1.13 and 1.17, respectively. Further investigation reveals that CR at lower supersaturation perform better than that at higher supersaturation. The dependence of CR on organic mass fraction (MForg) varies by environment: for marine aerosols, CR decreases with increasing organic fraction at lower supersaturations, whereas continental aerosols exhibit a consistent overestimation, with CR decreasing as organic fraction increases at higher supersaturations. This study provides key insights into CCN characteristics over midlatitude marine and continental environments, emphasizing the necessity of incorporating size-resolved chemical composition and mixing states into future model parameterizations, and contributing to a better understanding of aerosol–cloud interactions.

ACE-ENA field campaign↗

High-voltage sheaths and charge neutralization in space power systems

The authors examine the electrodynamics of charged platforms in the ionosphere with a variety of analytical and numerical models. These models have been specifically designed to study the tethered satellite system (TSS-1) due to be launched in early 1992. One of the objectives of TSS-1 is to determine the potential of tethers for electrical power generation from orbital motion across the earth's magnetic field. The author identifies and explores important aspects of the interaction between the ambient ionospheric plasma and moving charged orbital platforms: (1) the formation of energetic particles in the wake of a nominally neutral satellite, (2) transient current collection by a highly charged platform, (3) the current closure paths in the ionosphere between multiple polarized platforms, and (4) the conditions for rapid neutralization by enhanced plasma formation in the presence of effluent gases.

Satyanarayana, P.↗

Large Eddy Simulation Of Gravitational Effects In Transitional And Turbulent Gas-Jet Diffusion Flames

The influence of gravity on the spatial and the compositional structures of transitional and turbulent hydrocarbon diffusion flames are studies via large eddy simulation (LES) and direct numerical simulation (DNS) of round and planar jets. The subgrid-scale (SGS) closures in LES are based on the filtered mass density function (FMDF) methodology. The FMDF represents the joint probability density function (PDF) of the SGS scalars, and is obtained by solving its transport equation. The fundamental advantage of LES/FMDF is that it accounts for the effects of chemical reaction and buoyancy exactly. The methodology is employed for capturing some of the fundamental influences of gravity in equilibrium flames via realistic chemical kinetic schemes. Some preliminary investigation of the gravity effects in non-equilibrium flames is also conducted, but with idealized chemical kinetics models.

Jaberi, Farhad A.↗

ST5 Observations of the Imbalance of Region 1 and 2 Field-Aligned Currents and Its Implication to the Cross-Polar Cap Pedersen Currents

In this study, we use the in-situ magnetic field observations from Space Technology 5 mission to quantify the imbalance of Region 1 (R1) and Region 2 (R2) currents. During the three-month duration of the ST5 mission, geomagnetic conditions range from quiet to moderately active. We find that the R1 current intensity is consistently stronger than the R2 current intensity both for the dawnside and the duskside large-scale field-aligned current system. The net currents flowing into (out of) the ionosphere in the dawnside (duskside) are in the order of 5% of the total R1 currents. We also find that the net currents flowing into or out of the ionosphere are controlled by the solar wind-magnetosphere interaction in the same way as the field-aligned currents themselves are. Since the net currents due to the imbalance of the R1 and R2 currents require that their closure currents flow across the polar cap from dawn to dusk as Pedersen currents, our results indicate that the total amount of the cross-polar cap Pedersen currents is in the order of ~ 0.1 MA. This study, although with a very limited dataset, is one of the first attempts to quantify the cross-polar cap Pedersen currents. Given the importance of the Joule heating due to Pedersen currents to the high-latitude ionospheric electrodynamics, quantifying the cross-polar cap Pedersen currents and associated Joule heating is needed for developing models of the magnetosphere-ionosphere coupling.

Le, Guan↗