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,153 records · Page 64

Interactions Between the Amazonian Rainforest and Cumuli Clouds: A Large‐Eddy Simulation, High‐Resolution ECMWF, and Observational Intercomparison Study

Abstract The explicit coupling at meter and second scales of vegetation's responses to the atmospheric‐boundary layer dynamics drives a dynamic heterogeneity that influences canopy‐top fluxes and cloud formation. Focusing on a representative day during the Amazonian dry season, we investigate the diurnal cycle of energy, moisture and carbon dioxide at the canopy top, and the transition from clear to cloudy conditions. To this end, we compare results from a large‐eddy simulation technique, a high‐resolution global weather model, and a complete observational data set collected during the GoAmazon14/15 campaign. The overall model‐observation comparisons of radiation and canopy‐top fluxes, turbulence, and cloud dynamics are very satisfactory, with all the modeled variables lying within the standard deviation of the monthly aggregated observations. Our analysis indicates that the timing of the change in the daylight carbon exchange, from a sink to a source, remains uncertain and is probably related to the stomata closure caused by the increase in vapor pressure deficit during the afternoon. We demonstrate quantitatively that heat and moisture transport from the subcloud layer into the cloud layer are misrepresented by the global model, yielding low values of specific humidity and thermal instability above the cloud base. Finally, the numerical simulations and observational data are adequate settings for benchmarking more comprehensive studies of plant responses, microphysics, and radiation.

Vilà‐Guerau de Arellano, J.↗

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗

Bloom Modeling and Prediction of the Harmful Algae Alexandrium in Bellingham Bay, WA

Bellingham, Washington is located near the American-Canadian border in the northwestern part of Washington state, and is bordered on its west side by the Lummi reservation. Between Lummi and Bellingham lies Bellingham bay, which has a history of harmful algae related closures dating back to 1978. The subject of this work is a genus of dinoflagellates: Alexandrium, within which many species have the capacity to produce a suite of toxins known as saxitoxin. These toxins bioaccumulate in bivalves, which in turn cause paralytic shellfish poisoning in marine consumers (mammals, birds, and fish), including humans. Symptoms in humans can range from tingling and numbness to difficulty or inability to breathe, resulting in death. Because of the longstanding history of shellfish gathering among Salish tribes and the fact that Washington is the leading U.S. producer of farmed bivalves, harmful algae blooms impact both native and non-native peoples living in the Salish Sea area negatively. The objective of this work is to identify factors which influence Alexandrium blooms in Bellingham Bay, as well as predict Alexandrium blooms in the future. The method of doing so involves two processes: an initial statistical modeling phase to find in situ and remote sensing observations correlated to bloom density (including, but not limited to: water temperature, chlorophyll-a, salinity, color dissolved organic matter, and discharge rate of local rivers), followed by use of that data as a training set for a recursive neural network. This predictive capacity may inform future closures, help ensure the safety of shellfish consumers, and act as a baseline for future modeling efforts in the region.

Harmful↗

Constraining Aerosol Optical Models Using Ground-Based, Collocated Particle Size and Mass Measurements in Variable Air Mass Regimes During the 7-SEAS/Dongsha Experiment

During the spring of 2010, NASA Goddard's COMMIT ground-based mobile laboratory was stationed on Dongsha Island off the southwest coast of Taiwan, in preparation for the upcoming 2012 7-SEAS field campaign. The measurement period offered a unique opportunity for conducting detailed investigations of the optical properties of aerosols associated with different air mass regimes including background maritime and those contaminated by anthropogenic air pollution and mineral dust. What appears to be the first time for this region, a shortwave optical closure experiment for both scattering and absorption was attempted over a 12-day period during which aerosols exhibited the most change. Constraints to the optical model included combined SMPS and APS number concentration data for a continuum of fine and coarse-mode particle sizes up to PM2.5. We also take advantage of an IMPROVE chemical sampler to help constrain aerosol composition and mass partitioning of key elemental species including sea-salt, particulate organic matter, soil, non sea-salt sulphate, nitrate, and elemental carbon. Our results demonstrate that the observed aerosol scattering and absorption for these diverse air masses are reasonably captured by the model, where peak aerosol events and transitions between key aerosols types are evident. Signatures of heavy polluted aerosol composed mostly of ammonium and non sea-salt sulphate mixed with some dust with transitions to background sea-salt conditions are apparent in the absorption data, which is particularly reassuring owing to the large variability in the imaginary component of the refractive indices. Extinctive features at significantly smaller time scales than the one-day sample period of IMPROVE are more difficult to reproduce, as this requires further knowledge concerning the source apportionment of major chemical components in the model. Consistency between the measured and modeled optical parameters serves as an important link for advancing remote sensing and climate research studies in dynamic aerosol-rich environments like Dongsha.

Bell, Shaun W.↗

Verification and Validation of a Conceptual Model of the Auto-Rigging Payload Handling and Off-Loading System Using LEGO Technic System and Three-Dimensional Printed Parts

A conceptual model (CM) can be used to validate a concept in modeling and simulation life cycles. During the 2020-2022 Coronavirus disease 2019 (COVID-19) pandemic, for employee safety NASA implemented center closures and mandatory telework for the entire workforce. During this challenging time, engineers and researchers at NASA Langley Research Center (LaRC) looked for safe and innovative approaches and methods to continue the development of CMs for various projects. Engineers and researchers at LaRC researched Auto-Rigging Payload Handling and Off-Loading System (ARPHOLS) for payload handing and off-loading a system on an inclined lunar lander deck. In this paper, the development of a CM and the verification and validation of a conceptual idea for ARPHOLS using a LEGO Technic system and three-dimensional printed parts is presented.

Design↗

Scrutinizing entrainment and mass flux closures in shallow cumulus parameterizations using cloud-radar observations and large-eddy simulation (Final Report)

Representing warm, shallow clouds is a challenge for Global Climate Models, which leads to uncertainty in model projections of future climate. Our study was particularly concerned with how climate models represent the effects of shallow, non-precipitating cumulus clouds. Climate models must make certain assumptions when representing low clouds such as how much dry air from outside the cloud is mixed into the warm, moist cloud updraft. This process is called “entrainment” and tends to weaken the cloud updraft and dry out the cloud, compared to a theoretical cloud with no mixing. This project was a joint observational and modeling study to derive estimates of entrainment from field observations and use high-resolution modeling to evaluate the assumptions used in climate model representations of shallow cumulus. Our research took advantage of a wide variety of Department of Energy measurements from fixed climate-monitoring locations over the Southern Great Plains of the U.S. and the in the Eastern North Atlantic, along with short-term field projects (MAGIC, CAP–MBL, ACE–ENA). During the project, the PI maintained regular collaborations with Department of Energy Scientists at Brookhaven and Argonne National Laboratories. The project funded two students, who both obtained their M.S. in Atmospheric Science and are currently in the Atmospheric Science Ph.D. program. Both were lead authors on research papers funded by this grant.

54 ENVIRONMENTAL SCIENCES↗

Two-Neutron Halo is Unveiled in F 29

Here, we report the measurement of reaction cross sections ( σ R ex ) of F 27 , 29 with a carbon target at RIKEN. The unexpectedly large σ R ex and derived matter radius identify F 29 as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the 2 p 3 / 2 orbital, thereby vanishing the shell closure associated with the neutron number N = 20 . The findings are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of F 27 but are challenged for F 29 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Observations of Birkeland currents

Recent measurements of precipitating energetic particles and vector magnetic fields from satellites and sounding rockets have verified the existence of geomagnetically-aligned electric currents at high latitudes in the ionosphere and magnetosphere. The spatial and temporal configuration of such currents, now commonly called Birkeland currents, has delineated their role in providing ionospheric closure of magnetospheric current systems, and gross features of these current systems may be understood in terms of theoretical models of magnetospheric convection. The association of Birkeland currents with auroral features on a very small scale suggests that auroral acceleration may result from the current flow.

Cloutier, P. A.↗

Velocity-coupled flow oscillations in a simulated solid-propellant rocket environment

A comprehensive numerical analysis has been carried out to study the unsteady flowfields in a simulated rocket-motor environment. The model is based on the time-dependent compressible Navier-Stokes equations with a two-equation turbulence closure scheme. Various important aspects of the coupling between acoustic oscillations and mean flowfields, including flow reversal, modification of transport properties, etc., are addressed. Results indicate that multi-dimensional effects play important roles in determining local flow structures and wave characteristics. In much of the domain, acoustic velocity nodal points are observed in the near-wall region. The classical one-dimensional theory fails to describe several important mechanisms associated with velocity-induced flow instabilities.

Yang, Vigor↗

Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions. Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz model as well as more realistic models of the means and atmosphere. A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter situations to allow for correct update of the ensemble members. The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to be quite puzzling in that results state estimates are worse than for their filter analogue. In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use the Lorenz model to test and compare the behavior of a variety of implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqui↗

A13K-0336: Airborne Multi-Wavelength High Spectral Resolution Lidar for Process Studies and Assessment of Future Satellite Remote Sensing Concepts

NASA Langley recently developed the world's first airborne multi-wavelength high spectral resolution lidar (HSRL). This lidar employs the HSRL technique at 355 and 532 nm to make independent, unambiguous retrievals of aerosol extinction and backscatter. It also employs the standard backscatter technique at 1064 nm and is polarization-sensitive at all three wavelengths. This instrument, dubbed HSRL-2 (the secondgeneration HSRL developed by NASA Langley), is a prototype for the lidar on NASA's planned Aerosols- Clouds-Ecosystems (ACE) mission. HSRL-2 completed its first science mission in July 2012, the Two-Column Aerosol Project (TCAP) conducted by the Department of Energy (DOE) in Hyannis, MA. TCAP presents an excellent opportunity to assess some of the remote sensing concepts planned for ACE: HSRL-2 was deployed on the Langley King Air aircraft with another ACE-relevant instrument, the NASA GISS Research Scanning Polarimeter (RSP), and flights were closely coordinated with the DOE's Gulfstream-1 aircraft, which deployed a variety of in situ aerosol and trace gas instruments and the new Spectrometer for Sky-Scanning, Sun-Tracking Atmospheric Research (4STAR). The DOE also deployed their Atmospheric Radiation Measurement Mobile Facility and their Mobile Aerosol Observing System at a ground site located on the northeastern coast of Cape Cod for this mission. In this presentation we focus on the capabilities, data products, and applications of the new HSRL-2 instrument. Data products include aerosol extinction, backscatter, depolarization, and optical depth; aerosol type identification; mixed layer depth; and rangeresolved aerosol microphysical parameters (e.g., effective radius, index of refraction, single scatter albedo, and concentration). Applications include radiative closure studies, studies of aerosol direct and indirect effects, investigations of aerosol-cloud interactions, assessment of chemical transport models, air quality studies, present (e.g., CALIPSO) and future (e.g., EarthCARE) satellite calibration/validation, and development/assessment of advanced retrieval techniques for future satellite applications (e.g., lidar+polarimeter retrievals of aerosol and cloud properties). We will also discuss the relevance of HSRL-2 measurement capabilities to the ACE remote sensing concept.

Hostetler, Chris A.↗

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

Aft Body Closure: Predicted Strut Effects at M=2.4

This paper reports the predicted M = 2.4 strut-interference effects on a closed aftbody with empennage for the TCA baseline model. The strut mounting technique was needed in order to assess the impact of aft-end shaping, i.e. open for a sting or closed to better represent a flight vehicle. However,this technique can potentially lead to unanticipated effects that are measured on the aft body. Therefore, a set of computations were performed in order to examine the closed aft body with and without strut present, at both zero and non-zero angles of sideslip (AOS). The work was divided into a computational task performed by Javier A. Garriz, using an inviscid (Euler) solver, and a monitoring/reporting task done by John E. Lamar. All this work was performed during FY98 at the NASA Langley Research Center.

Lamar, John E.↗

Effect of the Addition of a Low Equivalent Stress Mechanism to the Analysis of Geomechanical Behavior of the SPR West Hackberry Site

Sandia National Laboratories has long used the Munson-Dawson (M-D) model to predict the geomechanical behavior of salt caverns used to store oil at the Strategic Petroleum Reserve (SPR). Salt creep causes storage caverns to deform inward, thus losing volume. This loss of volume affects the salt above and around the caverns, puts stresses and strains on borehole casings, and creates surface subsidence which affects surface infrastructure. Therefore, accurate evaluation of salt creep behavior drives decisions about cavern operations. Parameters for the M-D model are typically fit against laboratory creep tests, but nearly all historic creep tests have been performed at equivalent stresses of 8 MPa or higher. Creep rates at lower equivalent stresses are very slow, such that tests take months or years to run, and the tests are sensitive to small temperature perturbations (<0.1°C). A recent collaboration between US and German researchers, recently characterized the creep behavior at low equivalent (deviatoric) stresses (<8 MPa) of salt from the Waste Isolation Pilot Plant. In addition, the M-D model was recently extended to include a low stress creep “mechanism”. This paper details new simulations of SPR caverns that use this extended M-D model, called the M-D Viscoplastic model. The results show that the inclusion of low stress creep significantly alters the prediction of steady-state cavern closure behavior and indicate that low stress creep is the dominant displacement mechanism at the dome scale. The implications for evaluating cavern and well integrity are demonstrated by investigating three phenomena: the extent of stress changes around the cavern; the predicted vertical strains applied to wellbore casings; and the evaluation of oscillating stress changes around the cavern due to oil sale cycles and their potential effect on salt fatigue.

58 GEOSCIENCES↗

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗

Finite-element numerical modeling of atmospheric turbulent boundary layer

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Lee, H. N.↗

Anisotropy of thermal infrared exitance in sunflower canopies

Anisotropy of thermal infrared exitance above and within a relatively closed fully irrigated sunflower canopy is detailed. Azimuthal variation in thermal infrared exitance above canopies was weakly (statistically) related to solar position and was comparable to or larger than errors in satellite-based canopy estimates. Anisotropy within canopies was significantly lower and decreased with canopy closure and depth into the canopy. Measured azimuthal isotropy within canopies supports the use of this assumption in radiative transfer models. Significant differences in canopy temperature measurements were found depending upon whether the instruments were within or above the canopy. These differences could produce errors of 20-35 percent in latent energy estimates during periods of high evapotranspiration (ET) and greater errors in periods of restricted ET.

Tha Paw u, Kyaw↗