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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the significantly enhanced (u,v) coverage in the 2021 Event Horizon Telescope (EHT) observations. The integration of the 12−m Kitt Peak Telescope (USA) and NOEMA (France) stations into the array introduces two critical intermediate-length baselines to SMT (USA) and IRAM 30−m (Spain), providing sensitivity to emission structures at spatial scales of ∼250 μas and ∼2500 μas (∼ 0.02 pc and ∼ 0.02 pc). Without these new baselines, previous EHT observations of the source in 2017 and 2018 lacked the capability to constrain emission on large scales, where a “missing flux” of order ∼1 Jy is expected to reside. To probe these scales, we analyzed closure phases–robust against station-based gain calibration errors–and model the jet base emission using a simple Gaussian component offset from the compact ring emission at spatial separations > 100 μas. Our analysis revealed a Gaussian feature centered at (ΔRA ≈ 320 μ as, ΔDec. ≈ 60 μ as), projected separation of ≈ 5500 AU, with an estimated flux density of only ∼60 mJy, implying that most of the missing flux identified in previous EHT studies had to originate from different, larger scales. Brighter emission at the relevant spatial scales is firmly ruled out, and the data do not favor more complex models. This component aligns with the inferred position of the large-scale jet and is therefore physically consistent with the emission of the jet base. While our findings point to detectable jet base emission at 230 GHz, the limited coverage provided by only two intermediate baselines limits our ability to robustly reconstruct its morphology. Consequently, we treated the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate baseline coverage will be essential to constrain the structure and nature of this component with higher precision.

accretion, accretion disks↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Performance Evaluation of an Occupancy-Based HVAC Control System in an Office Building

As new algorithms incorporate occupancy count information into more sophisticated HVAC control, these technologies offer great potential for reductions in energy costs while enhancing flexibility. This study presents results from a two-year field evaluation of an occupancy-based HVAC control system installed in an office building. Two wings on each of the building’s 2–11 floors were equipped with occupancy counters to learn occupancy patterns. In combination with proprietary machine learning algorithms and thermal modeling, the occupancy data were leveraged to implement optimized start, early closure, and adjustments to fan operation at the air handling unit (AHU) level. This study conducted a holistic evaluation of technical performance, cost-effectiveness analysis, and user satisfaction. Results show the platform reduced weekday AHU run times by 2 h and 35 min per AHU per day during the pandemic time period. Simulation shows that 6.1% annual whole-building savings can be achieved when the building is fully occupied. The results are compared with prior studies, and potential drivers are discussed for future opportunities. The assessment results shed light on the expected in-the-field performance for researchers and industry stakeholders and enabled practical considerations as the technology strives to move beyond research-grade pilot trials into product-grade deployment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1977

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗
NASA NTRS1980

QSO evolution and the intergalactic medium

A model is developed of a cold, photoionized intergalactic medium in which only discrete, evolving sources (quasars) produce the observed extragalactic hard X-ray background. Constraints on the model are imposed by observed optical depth limits to various redshifts, far UV and visible flux measurements, X-ray background and source counts in the region 2-6 keV, and it makes use of results of the Palomar Bright Quasar Survey and the Braccesi Quasar survey to obtain expressions for the quasar evolution function, absolute luminosity, present number density and switch-on redshift. Results of computer modeling reveal an intergalactic medium with a density with respect to closure less than 0.35, temperature of 16,000 K and fractional ion abundances of approximately 0.000001 for H I and 0.001 for He II, under the assumption of a Hubble constant of 50 km/sec per Mpc. In addition, spectral indices of -0.4 and -1.3 are obtained for the regions of the emission spectrum above and below 4 x 10 to the 18th Hz.

Sherman, R. D.↗
NASA NTRS1988

Analysis of crack propagation in roller bearings using the boundary integral equation method - A mixed-mode loading problem

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions, making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

Ghosn, L. J.↗
NASA NTRS1993

Parameterization of sub-grid scale convection

The following topics are discussed: an overview of the cumulus parameterization problem; interactions between explicit and implicit processes in mesoscale models; effects of model grid size on the cumulus parameterization problem; parameterizing convective effects on momentum fields in mesoscale models; differences between slantwise and vertical cumulus parameterization; experiments with different closure hypotheses; and coupling cumulus parameterizations to boundary layer, stable cloud, and radiation schemes.

Frank, William↗
NASA NTRSDate not supplied

Coupling Online Conductivity with Offline Ion Chromatography Measurements using the Particle-into-Liquid Sampler

A particle-into-liquid sampler has been combined with a flow through conductivity cell to provide a continuous, non-destructive, online measurement in support of offline ion chromatography analysis. The conductivity measurement provides a rapid assessment of the total ion concentration augmenting the slower batch data from the offline analysis and is developed primarily to assist airborne measurements, where fast time response is essential. A model of the conductivity was developed for measured ions and excellent closure is derived for laboratory-generated aerosols. The PILS-conductivity measurement was extensively tested during the NASA Cloud, Aerosol and Monsoon Processes: Philippines Experiment (CAMP2Ex) across nineteen research flights and the conductivity data was found to augment the temporal capability of the PILS instrument for assessing ionic species, allowing sub-minute variability to be resolved. Through the sampling of a diverse range of ambient aerosol, including biomass burning, fresh and aged urban pollution, the conductivity measurement offered additional useful information to untangle the composition of complex aerosol mixtures, specifically for assessing acidic aerosol conditions.

Ewan Colin Crosbie↗
DOE OSTI2020

Verification and Validation of High Explosive Reactive Burn Models Implemented in LANL's EAP and LAP Code Base (FY2020 L2 Milestone MRT 7129) [Slides]

Advanced reactive burn models enable phenomena to be modeled that aren’t possible with programmed burn. Programmed burn propagates a burn front at a prescribed speed from a prescribed initiation point. Reactive burn uses the hydro variables (pressure, temperature, etc.) to determine when, how rapidly, and to what extent a zone detonates. This enables a completely new set of problems to be simulated. Our goal in this milestone has been to improve the usability of reactive burn modes. Thus, the closure criteria are: 1. Perform code and solution verification of the reactive burn models using appropriate solutions and simplified HE setups. 2. Perform validation assessments using small-scale experiments. 3. Evaluate sensitivities and quantify uncertainties due to model form variations and mesh dependencies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS2018

Preliminary Evaluation of the Simplified Higher Order Closure (SHOC) in the NASA Goddard Earth Observing System (GEOS)

Here we evaluate alternative planetary boundary layer (PBL) schemes in the NASA Goddard Earth Observing System (GEOS) model. The model currently calculates turbulent diffusivities with a local Richardson number-based scheme for stable layers, and a K-profile scheme with surface- and cloud-driven plumes based on Lock et al. (2000) for unstable layers. In this study, these are replaced with the TKE-based Simplified Higher Order Closure (SHOC). SHOC combines a 1.5-order turbulence closure with a double-Gaussian PDF serving as a subgrid condensation scheme. The new diffusivity calculations are paired with two new shallow mass flux schemes. First, the University of Washington shallow cumulus scheme of Park and Bretherton (2009), a buoyancy-sorting entraining-detraining bulk plume with overshooting entrainment and a TKE/CIN-based closure. Second, an Eddy Diffusivity Mass Flux (EDMF) scheme based on the multiple plume model of Suselj et al. (2013) with stochastic lateral entrainment. Global climatologies of low level temperature, humidity, and wind are evaluated against MERRA-2 and ERA-Interim reanalyses, and significant reductions in bias are seen with both mass flux schemes. The diurnal cycle over the continental United States is evaluated using MERRA-2 and an in situ network of surface Lidar profilers.

Arnold, Nathan↗
DOE OSTI2023

Data-Driven RANS Turbulence Closures for Forced Convection Flow in Reactor Downcomer Geometry

Recent progress in data-driven turbulence modeling has shown its potential to enhance or replace traditional equation-based Reynolds-averaged Navier-Stokes (RANS) turbulence models. Here, this work utilizes invariant neural network (NN) architectures to model Reynolds stresses and turbulent heat fluxes in forced convection flows (when the models can be decoupled). As the considered flow is statistically one dimensional, the invariant NN architecture for the Reynolds stress model reduces to the linear eddy viscosity model. To develop the data-driven models, direct numerical and RANS simulations in vertical planar channel geometry mimicking a part of the reactor downcomer are performed. Different conditions and fluids relevant to advanced reactors (sodium, lead, unitary-Prandtl-number fluid, and molten salt) constitute the training database. The models enabled accurate predictions of velocity and temperature, and compared to the baseline k–τ turbulence model with the simple gradient diffusion hypothesis, do not require tuning of the turbulent Prandtl number. The data-driven framework is implemented in the open-source graphics processing unit–accelerated spectral element solver nekRS and has shown the potential for future developments and consideration of more complex mixed convection flows.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

VERSE-LC TCCR Model Calibration using Batch 1A, Batch 2, and Batch 3 Column Data

The objective of this work is to utilize on-line gamma monitoring data collected during Tank Closure Cesium Removal (TCCR) Column Batch 3 operations to calibrate the parameters (essentially, the correction/dilution factors (CF or DF)) used in the VERSE-LC models for the prediction of cesium breakthrough while processing Savannah River Site (SRS) Tank 10H dissolved saltcake High-Level Waste through a packed cylindrical bed of Crystalline Silicotitanate (CST). Calibrated VERSE-LC models could be applied to aid in evaluating future TCCR operations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Stratus: An interactive steady state mixed layer model for personal computers

A steady-state, horizontally homogeneous, cloud-topped marine boundary layer model based primarily on the work of Lilly (1968) and Schubert et al., (1979) is presented. The conservative thermodynamic variables are equivalent potential temperature, theta(sub e), and total water mixing ratio, q + l. Some of the differences between this and Lilly's (1968) model are: radiation is allowed to penetrate into the boundary layer; cloud top values of longwave radiation, equivalent potential temperature, and water vapor mixing ratio are linear functions of height derived from climatological data at California coastal stations; and the closure assumption assumes a weighted average of Lilly's (1968) maximum and minimum entrainment theories. This model was programmed in FORTRAN and will run interactively on an IBM-compatible personal computer. The program allows the user to specify the geographical location, the wind speed, the sea-surface temperature, the large scale horizontal divergence, and the initial guess for cloud top height. Output includes the steady state values of cloud top and cloud base height, mixed layer equivalent potential temperature and total water mixing ratio, and the associated convective and radiative fluxes.

Guinn, Thomas A.↗
NASA NTRS1994

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗
DOE OSTI2024

Cost effectiveness of preemptive school closures to mitigate pandemic influenza outbreaks of differing severity in the United States

Background: Nonpharmaceutical interventions (NPIs) may be considered as part of national pandemic preparedness as a first line defense against influenza pandemics. Preemptive school closures (PSCs) are an NPI reserved for severe pandemics and are highly effective in slowing influenza spread but have unintended consequences. Methods: We used results of simulated PSC impacts for a 1957-like pandemic (i.e., an influenza pandemic with a high case fatality rate) to estimate population health impacts and quantify PSC costs at the national level using three geographical scales, four closure durations, and three dismissal decision criteria (i.e., the number of cases detected to trigger closures). At the Chicago regional level, we also used results from simulated 1957-like, 1968-like, and 2009-like pandemics. Our net estimated economic impacts resulted from educational productivity costs plus loss of income associated with providing childcare during closures after netting out productivity gains from averted influenza illness based on the number of cases and deaths for each mitigation strategy. Results: For the 1957-like, national-level model, estimated net PSC costs and averted cases ranged from $\$7.5$ billion (2016 USD) averting 14.5 million cases for two-week, community-level closures to $\$97$ billion averting 47 million cases for 12-week, county-level closures. We found that 2-week school-by-school PSCs had the lowest cost per discounted life-year gained compared to county-wide or school district–wide closures for both the national and Chicago regional-level analyses of all pandemics. The feasibility of spatiotemporally precise triggering is questionable for most locales. Theoretically, this would be an attractive early option to allow more time to assess transmissibility and severity of a novel influenza virus. However, we also found that county-wide PSCs of longer durations (8 to 12 weeks) could avert the most cases (31–47 million) and deaths (105,000–156,000); however, the net cost would be considerably greater ($\$88$-$\$103$ billion net of averted illness costs) for the national-level, 1957-like analysis. Conclusions: We found that the net costs per death averted ($\$180,000$-$\$4.2$ million) for the national-level, 1957-like scenarios were generally less than the range of values recommended for regulatory impact analyses ($\$4.6$ to 15.0 million). This suggests that the economic benefits of national-level PSC strategies could exceed the costs of these interventions during future pandemics with highly transmissible strains with high case fatality rates. In contrast, the PSC outcomes for regional models of the 1968-like and 2009-like pandemics were less likely to be cost effective; more targeted and shorter duration closures would be recommended for these pandemics.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗
NASA NTRS2002

Performance of advanced trauma life support procedures in microgravity

BACKGROUND: Medical operations on the International Space Station will emphasize the stabilization and transport of critically injured personnel and so will need to be capable of advanced trauma life support (ATLS). METHODS: We evaluated the ATLS invasive procedures in the microgravity environment of parabolic flight using a porcine animal model. Included in the procedures evaluated were artificial ventilation, intravenous infusion, laceration closure, tracheostomy, Foley catheter drainage, chest tube insertion, peritoneal lavage, and the use of telemedicine methods for procedural direction. RESULTS: Artificial ventilation was performed and appeared to be unaltered from the 1-G environment. Intravenous infusion, laceration closure, percutaneous dilational tracheostomy, and Foley catheter drainage were achieved without difficulty. Chest tube insertion and drainage were performed with no more difficulty than in the 1-G environment due to the ability to restrain patient, operator and supplies. A Heimlich valve and Sorenson drainage system were both used to provide for chest tube drainage collection with minimal equipment, without the risk of atmospheric contamination, and with the capability to auto-transfuse blood drained from a hemothorax. The use of telemedicine in chest tube insertion was demonstrated to be useful and feasible. Peritoneal lavage using a percutaneous technique, although requiring less training to perform, was found to be dangerous in weightlessness due to the additional pressure of the bowel on the anterior abdominal wall creating a high risk of bowel perforation. CONCLUSIONS: The performance of ATLS procedures in microgravity appears to be feasible with the exception of diagnostic peritoneal lavage. Minor modifications to equipment and techniques are required in microgravity to effect surgical drainage in the presence of altered fluid dynamics, to prevent atmospheric contamination, and to provide for the restraint requirements. A parabolic simulation system was developed for equipment and procedure verification, physiological research, and possible crew medical officer training in the future.

Parabolic Flight↗