Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “velocity model evaluation”

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 883 records · Page 49

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗

Development and Evaluation of Velocity Predictive Optimal Energy Management Strategies in Intelligent and Connected Hybrid Electric Vehicles

In this study, a thorough and definitive evaluation of Predictive Optimal Energy Management Strategy (POEMS) applications in connected vehicles using 10 to 20 s predicted velocity is conducted for a Hybrid Electric Vehicle (HEV). The presented methodology includes synchronous datasets gathered in Fort Collins, Colorado using a test vehicle equipped with sensors to measure ego vehicle position and motion and that of surrounding objects as well as receive Infrastructure to Vehicle (I2V) information. These datasets are utilized to compare the effect of different signal categories on prediction fidelity for different prediction horizons within a POEMS framework. Multiple artificial intelligence (AI) and machine learning (ML) algorithms use the collected data to output future vehicle velocity prediction models. The effects of different combinations of signals and different models on prediction fidelity in various prediction windows are explored. All of these combinations are ultimately addressed where the rubber meets the road: fuel economy (FE) enabled from POEMS. FE optimization is performed using Model Predictive Control (MPC) with a Dynamic Programming (DP) optimizer. FE improvements from MPC control at various prediction time horizons are compared to that of full-cycle DP. All FE results are determined using high-fidelity simulations of an Autonomie 2010 Toyota Prius model. The full-cycle DP POEMS provides the theoretical upper limit on fuel economy (FE) improvement achievable with POEMS but is not currently practical for real-world implementation. Perfect prediction MPC (PP-MPC) represents the upper limit of FE improvement practically achievable with POEMS. Real-Prediction MPC (RP-MPC) can provide nearly equivalent FE improvement when used with high-fidelity predictions. Constant-Velocity MPC (CV-MPC) uses a constant speed prediction and serves as a “null” POEMS. Results showed that RP-MPC, enabled by high-fidelity ego future speed prediction, led to significant FE improvement over baseline nearly matching that of PP-MPC.

33 ADVANCED PROPULSION SYSTEMS↗

Hot gas ingestion characteristics and flow visualization of a vectored thrust STOVL concept

A 9.2 percent scale short takeoff and vertical landing (STOVL) hot gas ingestion model was designed and built by McDonnell Douglas Corporation (MCAIR) and tested in the NASA Lewis Research Center 9- by 15-Foot Low Speed Wind Tunnel (LSWT). Hot gas ingestion, the entrainment of heated engine exhaust into the inlet flow field, is a key development issue for advanced short takeoff and vertical landing aircraft. The Phase 1 test program, conducted by NASA Lewis and McDonnell Douglas Corporation, evaluated the hot ingestion phenomena and control techniques and Phase 2 test program which was conducted by NASA Lewis are both reported. The Phase 2 program was conducted at exhaust nozzles temperatures up to 1460 R and utilized a sheet laser system for flow visualization of the model flow field in and out of ground effects. Hot gas ingestion levels were measured for the several forward nozzle splay configurations and with flow control/lift improvement devices which reduced the hot gas ingestion. The model support system had four degrees of freedom, heated high pressure air for nozzle flow, and a suction system exhaust for inlet flow. The headwind (freestream) velocity for Phase 1 was varied from 8 to 90 kn, with primary data taken in the 8 to 23 kn headwind velocity range. Phase 2 headwind velocity varied from 10 to 23 kn. Results of both Phase 1 and 2 are presented. A description of the model, facility, a new model support system, and a sheet laser illumination system are also provided. Results are presented over a range of main landing gear height (model height) above the ground plane at a 10 kn headwind velocity. The results contain the compressor face pressure and temperature distortions, total pressure recovery, compressor face temperature rise, and the environmental effects of the hot gas. The environmental effects include the ground plane temperature and pressure distributions, model airframe heating, and the location of the ground flow separation. Results from the sheet laser flow visualization test are also shown.

Johns, Albert L.↗

Electron impact polarization expected in solar EUV lines from flaring chromospheres/transition regions

We have evaluated lower bounds on the degree of impact Extreme Ultraviolet/Ultraviolet (EUV/UV) line polarization expected during solar flares. This polarization arises from collisional excitation by energetic electrons with non-Maxwellian velocity distributions. Linear polarization was observed in the S I 1437 A line by the Ultraviolet Spectrometer and Polarimeter/Solar Maximum Mission (UVSP/SMM) during a flare on 15 July 1980. An early interpretation suggested that impact excitation by electrons propagating through the steep temperature gradient of the flaring transition region/high chromosphere produced this polarization. Our calculations show that the observed polarization in this UV line cannot be due to this effect. We find instead that, in some flare models, the energetic electrons can produce an impact polarization of a few percent in EUV neutral helium lines (i.e., lambda lambda 522, 537, and 584 A).

Fineschi, S.↗

Improved Finite Element Modeling of the Turbofan Engine Inlet Radiation Problem

Improvements have been made in the finite element model of the acoustic radiated field from a turbofan engine inlet in the presence of a mean flow. The problem of acoustic radiation from a turbofan engine inlet is difficult to model numerically because of the large domain and high frequencies involved. A numerical model with conventional finite elements in the near field and wave envelope elements in the far field has been constructed. By employing an irrotational mean flow assumption, both the mean flow and the acoustic perturbation problem have been posed in an axisymmetric formulation in terms of the velocity potential; thereby minimizing computer storage and time requirements. The finite element mesh has been altered in search of an improved solution. The mean flow problem has been reformulated with new boundary conditions to make it theoretically rigorous. The sound source at the fan face has been modeled as a combination of positive and negative propagating duct eigenfunctions. Therefore, a finite element duct eigenvalue problem has been solved on the fan face and the resulting modal matrix has been used to implement a source boundary condition on the fan face in the acoustic radiation problem. In the post processing of the solution, the acoustic pressure has been evaluated at Gauss points inside the elements and the nodal pressure values have been interpolated from them. This has significantly improved the results. The effect of the geometric position of the transition circle between conventional finite elements and wave envelope elements has been studied and it has been found that the transition can be made nearer to the inlet than previously assumed.

Roy, Indranil Danda↗

Developments in Human Centered Cueing Algorithms for Control of Flight Simulator Motion Systems

The authors conducted further research with cueing algorithms for control of flight simulator motion systems. A variation of the so-called optimal algorithm was formulated using simulated aircraft angular velocity input as a basis. Models of the human vestibular sensation system, i.e. the semicircular canals and otoliths, are incorporated within the algorithm. Comparisons of angular velocity cueing responses showed a significant improvement over a formulation using angular acceleration input. Results also compared favorably with the coordinated adaptive washout algorithm, yielding similar results for angular velocity cues while eliminating false cues and reducing the tilt rate for longitudinal cues. These results were confirmed in piloted tests on the current motion system at NASA-Langley, the Visual Motion Simulator (VMS). Proposed future developments by the authors in cueing algorithms are revealed. The new motion system, the Cockpit Motion Facility (CMF), where the final evaluation of the cueing algorithms will be conducted, is also described.

Houck, Jacob A.↗

Satellite Soil Moisture Data Assimilation Impacts on Modeling Weather Variables and Ozone in the Southeastern US – Part 2: Sensitivity to Dry-Deposition Parameterizations

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

Min Huang↗

Development of biomechanical models for human factors evaluations

Previewing human capabilities in a computer-aided engineering mode has assisted greatly in planning well-designed systems without the cost and time involved in mockups and engineering models. To date, the computer models have focused on such variables as field of view, accessibility and fit, and reach envelopes. Program outputs have matured from simple static pictures to animations viewable from any eyepoint. However, while kinematics models are available, there are few biomechanical models available for estimating strength and motion patterns. Those, such as Crew Chief, that are available are based on strength measurements taken in specific positions. Johnson Space Center is pursuing a biomechanical model which will use strength data collected on single joints at two or three velocities to attempt to predict compound motions of several joint simultaneously and the resulting force at the end effector. Two lines of research are coming together to produce this result. One is an attempt to use optimal control theory to predict joint motion in complex motions, and another is the development of graphical representation of human capabilities. The progress to date in this research is described.

Woolford, Barbara↗

The study of PDF turbulence models in combustion

In combustion computations, it is known that the predictions of chemical reaction rates are poor if conventional turbulence models are used. The probability density function (pdf) method seems to be the only alternative that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus is the only viable approach for more accurate turbulent combustion calculations. The fact that the pdf equation has a very large dimensionality renders finite difference schemes extremely demanding on computer memories and thus impractical. A logical alternative is the Monte Carlo scheme. Since CFD has a certain maturity as well as acceptance, it seems that the use of a combined CFD and Monte Carlo scheme is more beneficial. Therefore, a scheme is chosen that uses a conventional CFD flow solver in calculating the flow field properties such as velocity, pressure, etc., while the chemical reaction part is solved using a Monte Carlo scheme. The discharge of a heated turbulent plane jet into quiescent air was studied. Experimental data for this problem shows that when the temperature difference between the jet and the surrounding air is small, buoyancy effect can be neglected and the temperature can be treated as a passive scalar. The fact that jet flows have a self-similar solution lends convenience in the modeling study. Futhermore, the existence of experimental data for turbulent shear stress and temperature variance make the case ideal for the testing of pdf models wherein these values can be directly evaluated.

Hsu, Andrew T.↗

Evaluation of a New Cloud Droplet Activation Parameterization wtih in Situ Data from CRYSTAL-FACE and CSTRIPE

The accuracy of the 2003 prognostic, physically based aerosol activation parameterization of A. Nenes and J. H. Seinfeld (NS) with modification introduced by C. Fountoukis and A. Nenes in 2005 (modified NS) is evaluated against extensive microphysical data sets collected on board the Center for Interdisciplinary Remotely Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft for cumuliform and stratiform clouds of marine and continental origin. The cumuliform cloud data were collected during NASA's Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE, Key West, Florida, July 2002), while the stratiform cloud data were gathered during Coastal Stratocumulus Imposed Perturbation Experiment (CSTRIPE, Monterey, California, July 2003). In situ data sets of aerosol size distribution, chemical composition, and updraft velocities are used as input for the NS parameterization, and the evaluation is carried out by comparing predicted cloud droplet number concentrations (CDNC) with observations. This is the first known study in which a prognostic cloud droplet activation parameterization has been evaluated against a wide range of observations. On average, predicted droplet concentration in adiabatic regions is within -20% of observations at the base of cumuliform clouds and -30% of observations at different altitudes throughout the stratiform clouds, all within experimental uncertainty. Furthermore, CDNC is well parameterized using either a single mean updraft velocity w or by weighting droplet nucleation rates with a Gaussian probability density function of w. This study suggests that for nonprecipitating warm clouds of variable microphysics, aerosol composition, and size distribution the modified NS parameterization can accurately predict cloud droplet activation and can be successfully implemented for describing the aerosol activation process in global climate models.

Meskhidze, Nicholas↗

Development of optical diagnostics for performance evaluation of arcjet thrusters

Laser and optical emission-based measurements have been developed and implemented for use on low-power hydrogen arcjet thrusters and xenon-propelled electric thrusters. In the case of low power hydrogen arcjets, these laser induce fluorescence measurements constitute the first complete set of data that characterize the velocity and temperature field of such a device. The research performed under the auspices of this NASA grant includes laser-based measurements of atomic hydrogen velocity and translational temperature, ultraviolet absorption measurements of ground state atomic hydrogen, Raman scattering measurements of the electronic ground state of molecular hydrogen, and optical emission based measurements of electronically excited atomic hydrogen, electron number density, and electron temperature. In addition, we have developed a collisional-radiative model of atomic hydrogen for use in conjunction with magnetohydrodynamic models to predict the plasma radiative spectrum, and near-electrode plasma models to better understand current transfer from the electrodes to the plasma. In the final year of the grant, a new program aimed at developing diagnostics for xenon plasma thrusters was initiated, and results on the use of diode lasers for interrogating Hall accelerator plasmas has been presented at recent conferences.

Cappelli, Mark A.↗

Wind-Driven Montgolfiere Balloons for Mars

Solar Montgolfiere balloons, or solar-heated hot air balloons have been evaluated by use on Mars for about 5 years. In the past, JPL has developed thermal models that have been confirmed, as well as developed altitude control systems to allow the balloons to float over the landscape or carry ground sampling instrumentation. Pioneer Astronautics has developed and tested a landing system for Montgolfieres. JPL, together with GSSL. have successfully deployed small Montgolfieres (<15-m diameter) in the earth's stratosphere, where conditions are similar to a Mars deployment. Two larger Montgolfieres failed, however, and a series of larger scale Montgolfieres is now planned using stronger, more uniform polyethylene bilaminate, combined with stress-reducing ripstitch and reduced parachute deceleration velocities. This program, which is presently under way, is a joint effort between JPL, WFF, and GSSL, and is planned for completion in three years.

Mars↗

The Plasma Environment at Mercury

Mercury is the least explored terrestrial planet, and the one subjected to the highest flux of solar radiation in the heliosphere. Its highly dynamic, miniature magnetosphere contains ions from the exosphere and solar wind, and at times may allow solar wind ions to directly impact the planet's surface. Together these features create a plasma environment that shares many features with, but is nonetheless very different from, that of Earth. The first in situ measurements of plasma ions in the Mercury space environment were made only recently, by the Fast Imaging Plasma Spectrometer (FIPS) during the MESSENGER spacecraft's three flybys of the planet in 2008-2009 as the probe was en route to insertion into orbit about Mercury earlier this year. Here. we present analysis of flyby and early orbital mission data with novel techniques that address the particular challenges inherent in these measurements. First. spacecraft structures and sensor orientation limit the FIPS field of view and allow only partial sampling of velocity distribution functions. We use a software model of FIPS sampling in velocity space to explore these effects and recover bulk parameters under certain assumptions. Second, the low densities found in the Mercury magnetosphere result in a relatively low signal-to-noise ratio for many ions. To address this issue, we apply a kernel density spread function to guide removal of background counts according to a background-signature probability map. We then assign individual counts to particular ion species with a time-of-flight forward model, taking into account energy losses in the carbon foil and other physical behavior of ions within the instrument. Using these methods, we have derived bulk plasma properties and heavy ion composition and evaluated them in the context of the Mercury magnetosphere.

Raines, James M.↗

Model assessment of synthetic jets for turbulent combustion experiments

Abstract Understanding turbulent premixed flames is essential to predict and optimize advanced combustion strategies, but critical capability gaps exist for collecting and validating measurements such as turbulent flame speed. Here, we evaluate synthetic jets as a new, promising turbulence generation device for constant-volume combustion chambers, quantitatively assessing turbulence intensity and spatial uniformity in a hypothetical 4,189-cm 3 vessel for various premixture conditions. Graphic Abstract

33 ADVANCED PROPULSION SYSTEMS↗

Impact of molten salt inflow on the temperature distribution in thermal energy storage tanks at startup for central receiver concentrating solar power plants

Concentrating Solar Power (CSP) systems with molten salt thermal energy storage (TES) tanks are one of the most promising, renewable-based energy conversion technologies for larger-scale power generation. The TES tank is one of the most critical components in CSP plants due to its high-temperature operation (up to 565 °C), daily thermal cycling, and intermittent solar radiation conditions. The plant startup is one of the most challenging operation conditions that could lead to damaging thermal gradients due to low salt inventory levels. In this study an analytical model for the sparger ring was developed and integrated with a detailed computational fluid dynamics model of a commercial-scaled molten salt tank. The integrated model allows an accurate representation of the tank operation to evaluate the effect of molten salt inflow on the mixing process. The tank filling process during plant startup was analyzed considering sparger rings with variations in design features, including inlet orifice configurations, number of orifices, direction of the inlets, and orifice diameter. The results demonstrated that higher temperature gradients are obtained in the tank floor during the plant startup. A sparger ring configuration with a predetermined orifice inlet inclination (30°, 45° and 60°) leads to significant temperature differences in the floor, between 57 °C and 62 °C, but better homogeneity in the temperature of the salt inventory. Lower salt inflow velocities result in a more homogeneous floor temperature, with maximum temperature differences under 38 °C. The sparger ring configuration with 52 orifices of 1-in. diameter and vertical flow showed better homogeneity in the temperature differences as a function of the salt level and lower temperature gradients in the tank floor. Because large temperature gradients in the tank's floor have been identified as one of the main factors contributing to tank failures, assessing various sparger ring design features is fundamental to determining proper inflow conditions that lead to low-temperature gradients and reducing failure susceptibility.

14 SOLAR ENERGY↗

A simple model of a strong shock driven by a spherical or cylindrical piston

A simple model of piston-driven spherical and cylindrical shocks is suggested. The model is based on a consistent use of two factors: (a) an almost uniform pressure across the shocked layer and, (b) continuous geometrical stretching of the surface elements of the expanding piston. It turns out that for a uniform pre-shock medium the gas between the piston and the shock behaves essentially as an incompressible fluid. An algebraic equation for the shock vs piston position is obtained. Detailed evaluation of the accuracy of the proposed solution shows that its accuracy is a few percent for the adiabatic index γ = 5/3 (as in ideal plasma). A closed-form solution describing enhancement of a weak ambient magnetic field by the shock is presented. The proposed model of piston-driven shocks goes beyond the classical self-similar solutions in that it: (1) naturally covers an early, non-asymptotic dynamics and its transition to asymptotic regime; (2) allows for smooth radial density variation of an ambient gas of the form of bumps, dimples or ramps between two constant values; and (3) allows for smooth temporal variation of piston velocity of the form of bumps, dimples, or ramps. Furthermore, this simple and versatile model provides some new insights into a classical hydrodynamical problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗