SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Iterative solution techniques in boundary element analysis

Iterative techniques for the solution of the algebraic equations associated with the direct boundary element analysis (BEA) method are discussed. Continuum structural response analysis problems are considered, employing single- and multizone boundary element models with and without zone condensation. The impact on convergence rate and computer resource requirements associated with the sparse and blocked matrices, resulting in multizone BEA, is studied. Both conjugate gradient and generalized minimum residual preconditioned iterative solvers are applied for these problems and the performance of these algorithms is reported. Included is a quantification of the impact of the preconditioning utilized to render the boundary element matrices solvable by the respective iterative methods in a time competitive with direct methods.

Kane, J. H.↗

DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗

DOE OSTI2021

Investigation of potential aerosol transmission and infectivity of SARS-CoV-2 through central ventilation systems

Here we evaluate the concentrations and probabilities of infection for both building interior and exterior exposure sources using a well-mixed model in a connected multizone building. As a central hub of many community and economic activities, buildings provide social connectivity, but the COVID-19 pandemic has reduced social connectivity due to concerns of viral spread within buildings. Although single zone models of infectious spread are well studied, the impact of aerosolized spread of SARS-CoV-2 via air handling systems in multizone buildings remains unexplored. Here we evaluate the influence of filtration, air exchange rates, and the fraction of outdoor air on the probability of infection using the well-known well-mixed modeling approach for a multizone. We find filtration lowers the concentration and probability of infection in both source and connected rooms provided at least some air is recirculated, but that probability is not zero. Filtration has no influence without recirculation or unless the outdoor air contains virus. We find that increasing the air exchange rate removes virus from the source room faster but also increases the rate of exposure to connected rooms. Therefore, slower air exchange rates reduce infectivity in connected rooms at shorter durations, but higher air exchange rates reduce infectivity at longer durations. We further find that when outdoor air is virus free, increasing the fraction of outdoor air is helpful, but, when outdoor air is infective, pathogen exposure inside can persist for hours after a short-term release.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1995

Application of a Navier-Stokes Solver to the Analysis of Multielement Airfoils and Wings Using Multizonal Grid Techniques

A computational study was performed to determine the predictive capability of a Reynolds averaged Navier-Stokes code (CFL3D) for two-dimensional and three-dimensional multielement high-lift systems. Three configurations were analyzed: a three-element airfoil, a wing with a full span flap and a wing with a partial span flap. In order to accurately model these complex geometries, two different multizonal structured grid techniques were employed. For the airfoil and full span wing configurations, a chimera or overset grid technique was used. The results of the airfoil analysis illustrated that although the absolute values of lift were somewhat in error, the code was able to predict reasonably well the variation with Reynolds number and flap position. The full span flap analysis demonstrated good agreement with experimental surface pressure data over the wing and flap. Multiblock patched grids were used to model the partial span flap wing. A modification to an existing patched- grid algorithm was required to analyze the configuration as modeled. Comparisons with experimental data were very good, indicating the applicability of the patched-grid technique to analyses of these complex geometries.

Jones, Kenneth M.↗

DOE OSTI2026

Presupernova O–C Shell Mergers Could Produce More 44 Ti than the Explosion

The formation of 44 Ti in massive stars is thought to occur during explosive nucleosynthesis; however, recent studies have shown it can be produced during O–C shell mergers prior to core collapse. We investigate how mixing according to 3D macrophysics derived from hydrodynamic simulations impacts presupernova O–C shell merger nucleosynthesis and if it can dominate explosive supernova production of 44 Ti and other radioactive isotopes. We compare a range of observations and models of explosive 44 Ti yields to preexplosive multizone mixing–burning nucleosynthesis simulations of an O–C shell merger in a 15 M ⊙ Z = 0.02 stellar model with mixing conditions corresponding to different 3D hydrodynamic mixing scenarios. Radioactive species produced in the O shell have a multidex spread in preexplosive yield predictions across different 3D mixing scenarios of 1.54 and 2.14 dex on average, depending on mass cut. 44 Ti has the largest spread of 4.78 and 4.81 dex depending on mass cut. Further, we show that the preexplosive production of 44 Ti can be larger than the explosive production of models and can match observations. Our results also show that 3D mixing physics enhances 44 Ti in 1D models without modifying 56 Ni yields. We conclude that quantitative predictions of 44 Ti and other radioactive species more broadly require an understanding of the 3D hydrodynamic mixing conditions present during the O–C shell merger.

Issa, Joshua [University of Victoria, BC (Canada);↗

DOE OSTI2023

Evaluation of thermostat location for multizone commercial building performance

In multi-zone buildings, it is often found that a single shared thermostat controls more than one conditioned zones. Although these shared zones are supposed to have similar thermal needs (e.g., cooling and heating load), in reality, they are not mainly due to different orientations, sizes of windows, occupancy, space types, etc. This can cause unnecessary energy waste or thermal discomfort for the occupants. How to quantify this impact in multizone buildings remains a research gap. Therefore, this study aims to evaluate the impact of different sensor (i.e., thermostat) locations for multizone commercial buildings through a comprehensive modeling study. Here, two different scenarios for the sensor locations were selected to evaluate the impact in terms of energy and thermal comfort. The scenario (1) is that one to five sensors distributed among the five zones, but the sensor readings from selected zones will be used for no-sensor zones, which is no-mean sensor scenario. The scenario (2) is that one to five sensors distributed among the five zones, but the average temperature from the shared zones will be used for each of the shared zones, which is a mean sensor scenario. The uncertainty analysis was performed for different sensor location scenarios.(a)The major findings from an energy perspective, for scenario (1), the differences of cooling energy go as high as 17% more or 12% less, compared with the baseline. For heating energy consumption, the discrepancies go as high as 51% more or 52% less, compared with baseline. For site energy consumption, the discrepancies go as high as 3.2% more or 3.2% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.2% more or as low as 1.0% less, compared with baseline. For scenario (2), the discrepancies of cooling energy go as high as 3% more, or 0.5% less, compared with the baseline. For heating energy consumption, the discrepancies, are go as high as 10.1% more or as low as 3.0% less, compared with baseline. For site energy consumption, the discrepancies go as high as 1.3% more or 0.3% less, compared with baseline. For fan energy consumption, the discrepancies go as high as 3.1% more or as low as 1.0% less, compared with baseline.(b) In terms of the indoor thermal comfort, for the no-mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 1,200 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 740 h, compared with the baseline. For the mean-sensor scenarios, the discrepancies of unmet hours for cooling mode can be as high as 750 h, compared with the baseline. The discrepancies of unmet hours for heating mode can be as high as 50 h, compared with the baseline.

42 ENGINEERING↗

NASA NTRS1990

Investigation of surface water behavior during glaze ice accretion

A series of experimental investigations that focused on isolating the primary factors that control the behavior of unfrozen surface water during glaze ice accretion were conducted. Detailed microvideo observations were made of glaze ice accretions on 2.54 cm diam cylinders in a closed-loop refrigerated wind tunnel. Distinct zones of surface water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film, a rough zone where surface tension effects caused coalescence of surface water into stationary beads, and a zone where surface water ran back as rivulets. The location of the transition from the smooth to the rough zone was found to migrate towards the stagnation point with time. Comparative tests were conducted to study the effect of the substrate thermal and roughness properties on ice accretion. The importance of surface water behavior was evaluated by the addition of a surface tension reducing agent to the icing tunnel water supply, which significantly altered the accreted glaze ice shape. Measurements were made to determine the contact angle behavior of water droplets on ice. A simple multizone modification to current glaze ice accretion models was proposed to include the observed surface roughness behavior.

Hansman, R. John, Jr.↗

NASA NTRS2003

Defect, Kinetics and Heat Transfer of CDTE Bridgman Growth without Wall Contact

A detached growth mechanism has been proposed, which is similar to that proposed by Duffar et al. and used to study the current detached growth system. From numerical results, we can conclude that detached growth will more likely appear if the growth and wetting angles are large and meniscus is flat. Detached thickness is dependent on growth angle, wetting angle, and gap width and shape of the fins. The model can also explain why the detached growth will not happen for metals in which the growth angle is almost zero. Since the growth angle of CdZnTe cannot be changed, to promote detached growth, the number density of the fins should be low and the wetting angle should be high. Also, a much smaller gap width of the fins should be used in the ground experiment and the detached gap width is much smaller. The shape of the fins has minor influence on detached growth. An integrated numerical model for detached solidification has been developed combining a global heat transfer sub-model and a wall contact sub-model. The global heat transfer sub-model accounts for heat and mass transfer in the multiphase system, convection in the melt, macro-segregation, and interface dynamics. The location and dynamics of the solidification interface are accurately tracked by a multizone adaptive grid generation scheme. The wall contact sub-model accounts for the meniscus dynamics at the three-phase boundary. Simulations have been performed for crystal growth in a conventional ampoule and a designed ampoule to understand the benefits of detached solidification and its impacts on crystalline structural quality, e.g., stoichiometry, macro-segregation, and stress. From simulation results, both the Grashof and Marangoni numbers will have significant effects on the shape of growth front, Zn concentration distribution, and radial segregation. The integrated model can be used in designing apparatus and determining the optimal geometry for detached solidification in space and on the ground.

Larson, D. J., Jr.↗

NASA NTRS1988

Investigation of surface water behavior during glaze ice accretion

Microvideo observations of glaze ice accretions on 1-in-diameter cylinders in a closed-loop refrigerated wind tunnel were obtained to study factors controlling the behavior of unfrozen surface water during glaze ice accretion. Three zones of surface water behavior were noted, each with a characteristic roughness. The effect of substrate thermal and roughness properties on ice accretions was also studied. The contact angle and hysteresis were found to increase sharply at temperatures just below 0 C, explaining the high resistance to motion of water beads observed on accreting glaze ice surfaces. Based on the results, a simple multizone modification to the current glaze ice accretion model is proposed.

Hansman, R. John, Jr.↗

NASA NTRS1990

A user guide for the EMTAC-MZ CFD code

The computer code (EMTAC-MZ) was applied to investigate the flow field over a variety of very complex three-dimensional (3-D) configurations across the Mach number range (subsonic, transonic, supersonic, and hypersonic flow). In the code, a finite volume, multizone implementation of high accuracy, total variation diminishing (TVD) formulation (based on Roe's scheme) is used to solve the unsteady Euler equations. In the supersonic regions of the flow, an infinitely large time step and a space-marching scheme is employed. A finite time step and a relaxation or 3-D approximate factorization method is used in subsonic flow regions. The multizone technique allows very complicated configurations to be modeled without geometry modifications, and can easily handle combined internal and external flow problems. An elliptic grid generation package is built into the EMTAC-MZ code. To generate the computational grid, only the surface geometry data are required. Results obtained for a variety of configurations, such as fighter-like configurations (F-14, AVSTOL), flow through inlet, multi-bodies (shuttle with external tank and SRBs), are reported and shown to be in good agreement with available experimental data.

Szema, Kuo-Yen↗

NASA NTRS1993

Effect of leading-edge porosity on blade-vortex interaction noise

The effect of the porous leading-edge of an airfoil on the blade-vortex interaction noise, which dominates far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near-field are calculated directly from the computation. The porosity effect on the surface is modeled. Results show leading-edge transpiration can suppress pressure fluctuations at the leading-edge during BVI, and consequently reduce the amplitude of propagating noise by 30 percent at maximum in the near-field. The effect of porosity factor on the noise level is also investigated.

Lee, Soogab↗

NASA NTRS1994

Reduction of blade-vortex interaction noise through porous leading edge

The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Lee, Soogab↗

DOE OSTI2026

ComStock Measure Documentation: Fan Static Pressure Reset for Multizone Variable Air Volume Systems

This report assesses the potential for nationwide adoption of a duct static pressure reset in MZ VAV systems in appropriate applications. Building on the 3-year End-Use Load Profiles project to calibrate and validate the U.S. Department of Energy’s ResStock™ and ComStock™ models, this work produces national datasets that enable cities, states, utilities, and other stakeholders to answer a broad range of questions regarding their commercial building stock. ComStock is a highly granular, bottom-up model that uses various data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual sub-hourly energy consumption of the commercial building stock across the United States. The “baseline” model intends to represent the U.S. commercial building stock as it existed in 2018. The methodology of the baseline model is discussed in the ComStock Reference Documentation. The goal of this work is to develop energy efficiency and demand flexibility measures that cover market-ready technologies and study their mass-adoption impact on the baseline building stock. “Measures” refers to various “what-if” scenarios that can be applied to buildings. The results for the baseline and measure scenario simulations are published in public datasets that provide insights into building stock characteristics, operational behaviors, utility bill impacts, and annual and sub-hourly energy usage by fuel type and end use. This report describes the modeling methodology for a single ComStock measure scenario—Fan Static Pressure Reset for Multizone Variable Air Volume (VAV) Systems—and briefly introduces key results. The full public dataset can be accessed on the ComStock data lake or via the Data Viewer at comstock.nrel.gov. The public dataset enables users to create custom aggregations of results for their use case (e.g., filter to a specific county or building type).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1990

Multizone accretional evolution of planetesimal swarms

The general features of a new numerical simulation of planetesimal accretion which models multiple heliocentric distance zones, together with a detailed model for the planetesimal size and orbital distribution in each zone, are described. A restricted version of this model which allows only a single heliocentric distance zone has been used to test the validity of the code by comparing with results from earlier authors when the same physical phenomena are included. Generally, very good agreement is found.

Spaute, D.↗

NASA NTRS1991

Boundary-element shape sensitivity analysis for thermal problems with nonlinear boundary conditions

Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems are presented that have closed-form solutions.

Kane, James H.↗

NASA NTRS1992

Boundary formulations for shape sensitivity of temperature dependent conductivity problems

Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.

Kane, James H.↗

NASA NTRS1999

User's Guide for ENSAERO_FE Parallel Finite Element Solver

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaces ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multizonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRANKOSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft model. It is used to create the stiffness matrices for each substructure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the substructure boundary nodes.

Eldred, Lloyd B.↗

DOE OSTI2021

Characterization of the indoor near-field aerosol transmission in a model commercial office building

To evaluate the exposure potential of infectious aerosols containing SARS-CoV-2 in an office building setting, synthetic test aerosols were used to experimental study airborne particle transmission in a multizone small office test building at the Oak Ridge National Laboratory. Nine measurement points in a single zone using active aerosol impactors report that the coefficient of variation of the time-averaged concentration is <10% in two campaigns and < 15% in one campaign, so a nearly well-mixed condition was noted. To understand the effect of HVAC system operation on the dynamic concentration of aerosols in office spaces, an aerosol transport model that includes factors such as outside air (OA) ratio, filtration, return air fraction, transport loss in air ducts, and particle deposition has been developed. The results of model fitting demonstrate strong agreement with experimental data. Overall, our investigation finds the return air fraction effects outweigh other mechanisms for the aerosol recirculation in this study, and the impact of air change rate (ACR) is more important than the small particle deposition for aerosol removal. Because ACR dominates the aerosol transport, the full model can be simplified to just one factor, the ACR, while maintaining an acceptable representation of the experimental data.

60 APPLIED LIFE SCIENCES↗