Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak boundary conditions”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Acoustic Source Localization in Aircraft Interiors Using Microphone Array Technologies

Using three microphone array configurations at two aircraft body stations on a Boeing 777-300ER flight test, the acoustic radiation characteristics of the sidewall and outboard floor system are investigated by experimental measurement. Analysis of the experimental data is performed using sound intensity calculations for closely spaced microphones, PATCH Inverse Boundary Element Nearfield Acoustic Holography, and Spherical Nearfield Acoustic Holography. Each method is compared assessing strengths and weaknesses, evaluating source identification capability for both broadband and narrowband sources, evaluating sources during transient and steady-state conditions, and quantifying field reconstruction continuity using multiple array positions.

Sklanka, Bernard J.↗
NASA NTRS2019

A New Method to Retrieve the Diurnal Variability of Planetary Boundary Layer Height from Lidar Under Different Thermodynamic Stability Conditions

The planetary boundary layer height (PBLH) is an important parameter for understanding the accumulation of pollutants and the dynamics of the lower atmosphere. Lidar has been used for tracking the evolution of PBLH by using aerosol backscatter as a tracer, assuming aerosol is generally well-mixed in the PBL; however, the validity of this assumption actually varies with atmospheric stability. This is demonstrated here for stable boundary layers (SBL), neutral boundary layers (NBL), and convective boundary layers (CBL) using an 8-year dataset of micropulse lidar (MPL) and radiosonde (RS) measurements at the ARM Southern Great Plains, and MPL at the GSFC site. Due to weak thermal convection and complex aerosol stratification, traditional gradient and wavelet methods can have difficulty capturing the diurnal PBLH variations in the morning and forenoon, as well as under stable conditions generally. A new method is developed that combines lidar-measured aerosol backscatter with a stability dependent model of PBLH temporal variation (DTDS). The latter helps “recalibrate” the PBLH in the presence of a residual aerosol layer that does not change in harmony with PBL diurnal variation. The hybrid method offers significantly improved PBLH detection, with better correlation and smaller biases, under most thermodynamic conditions, especially for SBL and CBL. Relying on the physical process of PBL diurnal development, different schemes are developed for growing, maintenance, and decaying periods. Comprehensive evaluation of this new method shows much better tracking of diurnal PBLH variation and significantly smaller biases under various pollution levels.

Tianning Su↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
NASA NTRS2022

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2024

Quasifragmentation functions in the massive Schwinger model

We introduce the concept of the quark quasifragmentation function (qFF) using an equal-time and spatially boosted form of the Collins-Soper fragmentation function where the out-meson fragment is replaced by the current asymptotic condition. We derive the qFF for a fermion in two-dimensional quantum electrodynamics (QED2) using the Kogut-Susskind Hamiltonian after a mapping onto spin qubits in a spatial lattice with open boundary conditions. This form is suitable for quantum computations. We compute the qFF by exact diagonalization of the spin Hamiltonian. The results are compared to the qFF following from the Drell-Levy-Yan result for QED2, both at strong and weak coupling, and to two-dimensional quantum chromodynamics in the lowest Fock approximation.

Fragmentation functions↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
NASA NTRS1984

Numerical simulations of the Ross Sea tides

Tidal currents below the floating Ross Ice shelf are reconstructed by using a numerical tidal model. They are predominantly diurnal, achieve maximum strength in regions near where the ice shelf runs aground, and are significantly enhanced by topographic Rossby wave propagation along the ice front. A comparison with observations of the vertical motion of the ice shelf surface indicates that the model reproduces the diurnal tidal characteristics within 20 percent. Similar agreement for the relatively weak semidiurnal tides was not obtained, and this calls attention to possible errors of the open boundary forcing obtained from global-ocean tidal simulations and to possible errors in mapping zones of ice shelf grounding. Air-sea contact below the ice shelf is eliminated by the thick ice cover. The dominant sub-ice-shelf circulation may thus be tidally induced. A preliminary assessment of sub-ice-shelf conditions based on the numerical tidal simulations suggests that (1) strong barotropic circulation is driven along the ice front and (2) tidal fronts may form in the sub-ice-shelf cavity where the water column is thin and where the buoyancy input is weak.

Macayeal, D. R.↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

Magnetotail Boundary Crossings at Lunar Distances: ARTEMIS Observations

We compare results from a preliminary analysis of two years of ARTEMIS magnetopause boundary crossings at lunar distances with available empirical models. We remove the effects of variable solar wind flow directions and aberration angles to study the magnetotail cross-section as a function of solar wind conditions. The average magnetopause distance from the central axis is 26 RE, but this distance ranges from 10 RE for high solar wind dynamic pressures and strong northward IMF orientations to 39 RE for low solar wind dynamic pressures and weak southward IMF orientations compared to the nominal solar wind conditions. The time-independent Howe and Binsack (1972) model describes the average location of the crossings very well. For high solar wind dynamic pressures, the Lu et al., (2011) model performs best, while for low pressures the Petrinec and Russell (1996) model gives the closest prediction. As predicted by theory and seen in past studies, the magnetotail cross-section is suggestive of prolate during intervals of strong IMF By, but oblate during intervals of strong IMF Bz. Any asymmetric variation of the tail boundary with respect to the sign of IMF By was not observed. The decreasing size of boundary with the increasing dynamic pressure was found when dynamic pressures are smaller than 2nPa. Although the scatter is larger, the tail size for pressures larger than 2nPa suggests a constant radius. The tail boundary size decreases as the strength of IMF Bz increases regardless of its polarity. However, it was also observed that an even stronger southward IMF Bz can cause larger magnetopause size in the presence of large dynamic pressures.

Magnetotail↗
NASA NTRS1991

Analysis of the linear stability of compressible boundary layers using the PSE

The application of linearized parabolic stability equations (PSE) to compressible flow is considered. The effect of mean-flow nonparallelism is found to be weak on 2D waves and strong on 3D waves. Results for a single choice of free-stream parameters that corresponds to the atmospheric conditions at 15,000 m above sea level are presented.

Bertolotti, F. P.↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Emerging low-cloud feedback and adjustment in global satellite observations

From mid-2003 to mid-2024, a global decrease in low-cloud amount enhanced the absorption of solar radiation by 0.22±0.07 W m −2 per decade (±1σ range), accelerating the energy imbalance trend during that period (0.44 W m −2 per decade). Through controlling factor analysis, here we show that the low-cloud trend is due to a combination of cloud feedback and adjustments to greenhouse gases and aerosols (respectively 0.09±0.02, 0.05±0.03, and 0.03±0.03 W m −2 per decade), which jointly account for 74 % of the trend. The contribution of natural climate variability is weak but uncertain (0.01±0.08 W m −2 per decade), owing to a poorly constrained trend in boundary-layer inversion strength. Importantly, the observed low-cloud radiative trend lies well within the range of values simulated by contemporary global climate models under conditions close to present day. Any systematic model error in the representation of present-day global energy imbalance trends is thus likely to originate in processes unrelated to low clouds.

Geosciences↗
NASA NTRS1992

Computational aspects of helicopter trim analysis and damping levels from Floquet theory

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Gaonkar, Gopal H.↗
NASA NTRS1993

Helicopter trim analysis by shooting and finite element methods with optimally damped Newton iterations

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used, and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Achar, N. S.↗
DOE OSTI2026

Influence of the as-built microstructure on the recrystallization of an additively manufactured Inconel939 Ni-based superalloy

This study investigates the influence of the as-built microstructure on the recrystallization (RX) behavior and mechanical properties of the Ni-based superalloy Inconel 939 produced by laser powder bed fusion (PBF-LB/M). Two distinct as-built microstructures were obtained by varying the hatch distance (h d ): a columnar, strongly textured condition (h d =50, termed h d 50) and an equiaxed, weakly textured condition (h d =70, termed h d 70)). Both were subjected to nine solution treatments combining three temperatures (1100, 1150, and 1200 °C) and three holding times (1, 4, and 8 h). Comprehensive microstructural characterization was conducted to assess grain morphology, texture, grain boundary character, dislocation density, and precipitate distribution. Recrystallization was found to be significantly slower than in cast counterparts, requiring higher temperatures and longer times for completion. The initial microstructure plays a decisive role: full RX was achieved only in hd70 specimens after treatment at 1200 °C for 8 h, whereas hd50 samples exhibited delayed and incomplete RX under identical conditions. This behavior is attributed to the finer grain size and higher fraction of high-angle grain boundaries in hd70, which promote recrystallization. Mechanical testing revealed that hd70 samples subjected to a 1200 °C/8 h treatment followed by standard double ageing show higher yield and tensile strengths across the investigated temperature range than both printed and cast Inconel939 processed under conventional conditions, albeit with slightly reduced ductility. The enhanced mechanical performance is attributed to the larger grain size, which limits grain boundary sliding. These results demonstrate the critical importance of controlling the as-built microstructure and tailoring post-processing strategies to optimize high-temperature performance of PBF-LB/M Inconel939.

Inconel939↗
NASA NTRS2012

Characterization of Turbulent Latent and Sensible Heat Flux Exchange Between the Atmosphere and Ocean in MERRA

Turbulent fluxes of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth's energy and water balance. Characterizing both the spatiotemporal variability and the fidelity of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. This study examines the veracity of the recently completed NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) product with respect to its representation of the surface turbulent heat fluxes. A validation of MERRA turbulent heat fluxes and near-surface bulk variables at local, high-resolution space and time scales is achieved by making comparisons to a large suite of direct observations. Both in situ and satellite-observed gridded surface heat flux estimates are employed to investigate the spatial and temporal variability of the surface fluxes with respect to their annual mean climatologies, their seasonal covariability of near-surface bulk parameters, and their representation of extremes. The impact of data assimilation on the near-surface parameters is assessed through evaluation of incremental analysis update tendencies produced by the assimilation procedure. It is found that MERRA turbulent surface heat fluxes are relatively accurate for typical conditions but have systematically weak vertical gradients in moisture and temperature and have a weaker covariability between the near-surface gradients and wind speed than found in observations. This results in an underestimate of the surface latent and sensible heat fluxes over the western boundary current and storm track regions. The assimilation of observations mostly acts to bring MERRA closer to observational products by increasing moisture and temperature near the surface and decreasing the near-surface wind speeds. The major patterns of spatial and temporal variability of the turbulent heat fluxes produced by MERRA compare favorably to observationally based estimates. However, MERRA is distinct in terms of amplitude. These results suggest that MERRA is likely to be a valuable resource for a number of research applications though, as with all turbulent flux estimates, systematic issues should be taken into account

Robert, J. Brent↗
NASA NTRS2012

Characterization of Turbulent Latent and Sensible Heat Flux Exchange Between the Atmosphere and Ocean in MERRA

Turbulent fluxes of heat and moisture across the atmosphere-ocean interface are fundamental components of the Earth s energy and water balance. Characterizing both the spatiotemporal variability and the fidelity of these exchanges of heat and moisture is critical to understanding the global water and energy cycle variations, quantifying atmosphere-ocean feedbacks, and improving model predictability. This study examines the veracity of the recently completed NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) product with respect to its representation of the surface turbulent heat fluxes. A validation of MERRA turbulent heat fluxes and near-surface bulk variables at local, high-resolution space and time scales is achieved by making comparisons to a large suite of direct observations. Both in situ and satellite-observed gridded surface heat flux estimates are employed to investigate the spatial and temporal variability of the surface fluxes with respect to their annual mean climatologies, their seasonal covariability of near-surface bulk parameters, and their representation of extremes. The impact of data assimilation on the near-surface parameters is assessed through evaluation of incremental analysis update tendencies produced by the assimilation procedure. It is found that MERRA turbulent surface heat fluxes are relatively accurate for typical conditions but have systematically weak vertical gradients in moisture and temperature and have a weaker covariability between the near-surface gradients and wind speed than found in observations. This results in an underestimate of the surface latent and sensible heat fluxes over the western boundary current and storm track regions. The assimilation of observations mostly acts to bring MERRA closer to observational products by increasing moisture and temperature near the surface and decreasing the near-surface wind speeds. The major patterns of spatial and temporal variability of the turbulent heat fluxes produced by MERRA compare favorably to observationally based estimates. However, MERRA is distinct in terms of amplitude. These results suggest that MERRA is likely to be a valuable resource for a number of research applications though, as with all turbulent flux estimates, systematic issues should be taken into account.

Roberts, J. Brent↗