Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Velocity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Daily and Seasonal Variations of Shallow Seismic Velocities in Southern California From Joint Analysis of H/V Ratios and Autocorrelations of Seismic Waveforms

Establishing a baseline of ongoing secular velocity variations at the subsurface can improve the accuracy of detecting and interpreting short-term velocity changes, and advance the understanding of observed seismic motions and the behavior of subsurface materials. Toward these goals, we develop and apply a deconvolved autocorrelation (DA) method to estimate regional daily and seasonal changes of seismic velocities in southern California. The DA method combines advantages of traditional autocorrelation and Horizontal-to-Vertical Spectral Ratio, and is used to analyze over 10 years of data recorded by 50 stations. The results indicate widespread daily and seasonal changes of up to 10% and 4%, respectively, in the top tens of meters of the crust. The thickness of the surface layer, distance from the coast, and topographic variations are important factors controlling the amplitudes of the resolved velocity variations. The results suggest that changes of soil moisture and thermoelastic strain are likely dominant factors affecting the daily and seasonal variations, respectively. The developed DA method can improve the accuracy and robustness of estimated changes of subsurface materials at other locations.

58 GEOSCIENCES↗
DOE OSTI2022

Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy

Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.

3D velocity models↗
DOE OSTI2023

Transmission electron microscopy of the rapid solidification microstructure evolution and solidification interface velocity determination in hypereutectic Al-20at.%Cu after laser melting

The evolution of rapid solidification (RS) microstructure and solidification interface velocity has been studied experimentally by in-situ transmission electron microscopy and postmortem characterization for hypereutectic Al-20at.%Cu (37 wt.%Cu) after laser melting. Four morphologically distinct regions formed via growth modes changing from θ-Al 2 Cu phase dendrites to eutectic cell growth, possibly α-Al dendrite growth (α-cell), banded growth, and α-Al plane front growth. Tendency for faceting and low capacity for solute trapping limited the θ-phase dendrite growth to solidification interface velocity v < 0.07 m/s. Consequently, formation of pro-eutectic micro-constituent was suppressed, and RS microstructure formation was dominated by eutectic, α-cell, and banded morphology grains for the Al-20at.%Cu alloy. Eutectic growth operated for interface velocity of 0.1 m/s ≤ v < 0.3 m/s, with a transition from regular lamellar, 2-λ and 1-λ mode to a dense irregular morphology dominated by α-phase at v = 0.3 m/s. Interface temperature calculations indicated feasibility of α-cell growth mode for 0.3 m/s ≤ v < 0.7 m/s. Banded growth occurred for 0.7 m/s ≤ v < 1.3 m/s. Plane front α-phase growth was evident for interface velocities v ≥ 1.3 m/s. Previous work on rapid solidification microstructure development in hypereutectic Al-Cu alloys reported a regime of α-cell growth subsequently to eutectic and prior to transition to banded growth for Al-19at.%Cu (36 wt.%Cu), while for Al-22at.%Cu (40 wt.%Cu) eutectic growth transitioned directly to α-plane front growth without emergence of a banded regime. Finally, based on the current study the disappearance of the banded growth regime occurs for a composition larger than 20at.%Cu.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data reduction considerations for the burning velocity of spherical constant volume flames of R32 (CH 2 F 2 ) with air

Here, the present work explores data reduction techniques for the measurement of the laminar burning velocities of R32(CH 2 F 2 )-air mixtures using a constant volume combustion device, in which the pressure-time history is the only measured parameter. To allow clear assessment of the accuracy of the data reduction methods, the pressure-time histories used for analysis are synthetically generated via a detailed numerical simulation employing full kinetics and with and without an optically-thin radiation model. Various data reduction models are employed, including a two-zone model and two multi-zone models, and these are compared with the results from the burning velocity obtained from the output of the numerical simulation. The data reduction schemes are shown to be accurate if the same radiation model is employed in the data reduction as was used in the flame simulation to generate the pressure trace used for post-processing. If the incorrect radiation model is employed, however, the errors can be quite large. The effects of stretch, radiation, and different data post-processing methodologies are explored and the errors quantified. Stretch is shown to be important for the early stages and the selected data range that is used for extrapolation has a significant effect on the extrapolated burning velocity. However, with an appropriate choice of data considered for extrapolation, the prediction of the unstretched burning velocity can be quite accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pervasive low-velocity layer atop the 410-km discontinuity beneath the northwest Pacific subduction zone: Implications for rheology and geodynamics

Regional triplication waveforms of five intermediate-depth events are modeled to simultaneously obtain the compressional (P) and shear (SH) wave velocity structure beneath northwestern Pacific subduction zone. Both the P- and SH-wave velocity models for three different sub-regions show a low-velocity layer (LVL) with a thickness of ~55-80 km lying above the 410-km discontinuity with a ~900 km lateral extent from the Japan Sea to the northeastern Asian continental margin. With the dihedral angle approaching to zero around 400 km, a minute amount of melt atop the 410-km discontinuity caused by the hydrous slab might completely wet olivine grain boundaries and result in a low seismic velocity layer in this specific subduction zone. This mechanism suggests that the 410-LVL is a low viscosity zone that would partially decouple the upper mantle from the transition zone. We infer that the widespread 410-LVL provides evidence for a water-bearing mantle transition zone beneath the western Pacific subduction zone.

58 GEOSCIENCES↗
DOE OSTI2023

Shear band velocity and activation volume during shear deformation by acoustic emission in a Zr-based bulk metallic glass

Here, an analysis is presented for the shear band velocity of plastic deformation and activation volume of the shear transformation zone (STZ) in a Zr 60 Cu 30 Al 10 bulk metallic glass (BMG). Compared to the values calculated by the stick-slip model, the values of shear band velocity calculated after optimization by the acoustic emission technique are much larger than those reported previously. Plastic flow in BMGs is found to be described by a power-law correlation between the maximum velocity of the shear band and the strain rates. Combining the free volume theory and the STZ theory, the relationship between the volume change in the STZ and the shear band velocity is derived as $Ω=A$(27–ln$v/l$). The value of A, which is a material dependent constant, is 5.43 × 10 –3 for the current BMG at room temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparative study between experimental measurements and model predictions of bubble rise velocity in molten LiCl-KCl

This research investigated the shape and rise velocities of gas bubbles (helium, nitrogen, argon, and krypton) in molten LiCl-KCl at 500 °C using a specialized optical visualization apparatus. The apparatus features a custom quartz rectangular prism cell housed in a furnace, enabling the application of the shadowgraph measurement technique. Bubbles with diameters between 0.45 and 3.2 mm were produced via a capillary, and their equivalent diameters and rise velocities were analyzed using high-resolution imaging and centroid tracking algorithms. Important non-dimensional parameters were calculated to characterize behavior. Bubbles with 0.03 < Eo < 1.32, 17 < Re < 718, and Mo order of 10 -11 were generated in the molten LiCl-KCl, revealing that bubbles with diameters less than 1.46 mm remain spherical, whereas those with larger diameters are ellipsoidal. The transition and oscillation onset for bubble shape occurred at larger diameters than those previously observed in air–water systems. Comparative analysis of the gas-molten salt data against existing air–water bubble rise velocity correlations highlighted inaccuracies in predictions, particularly across the different bubble regimes. An existing correlation was calibrated to the molten salt data to improve the bubble velocity predictions. In conclusion, this study contributes essential experimental data for the design and safety evaluation of molten salt reactors and offers insights for the optimization of sparge systems for fission product removal.

Bubble rise velocity↗
DOE OSTI2024

Ion velocity effect governs damage annealing process in defective KTaO 3

Effects of electronic to nuclear energy losses (S e /S n ) ratio on damage evolution in defective KTaO 3 have been investigated by irradiating pre-damaged single crystal KTaO 3 with intermediate energy O ions (6 MeV, 8 MeV and 12 MeV) at 300 K. By exploring these processes in pre-damaged KTaO 3 containing a fractional disorder level of 0.35, the results demonstrate the occurrence of a precursory stage of damage production before the onset of damage annealing process in defective KTaO 3 that decreases with O ion energy. The observed ionization-induced annealing process by ion channeling analysis has been further mirrored by high resolution transmission electron microscopy analysis. In addition, the reduction of disorder level is accompanied by the broadening of the disorder profiles to greater depth with increasing ion fluence, and enhanced migration is observed with decreasing O ion energy. Since S e (~3.0 keV nm –1 ) is nearly constant for all 3 ion energies across the pre-damaged depth, the difference in behavior is due to the so-called 'velocity effect': the lower ion velocity below the Bragg peak yields a confined spread of the electron cascade and hence an increased energy deposition density. Here, the inelastic thermal spike calculation has further confirmed the existence of a velocity effect, not previously reported in KTaO 3 or very scarcely reported in other materials for which the existence of ionization-induced annealing has been reported. In other words, understanding of ionization-induced annealing has been advanced by pointing out that ion velocity effect governs the healing of pre-existing defects, which may have significant implication for the creation of new functionalities in KTaO 3 through atomic-level control of microstructural modifications, but may not be limited to KTaO 3 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Wide-angle effects for peculiar velocities

ABSTRACT The line-of-sight peculiar velocities of galaxies contribute to their observed redshifts, breaking the translational invariance of galaxy clustering down to a rotational invariance around the observer. This becomes important when the line-of-sight direction varies significantly across a survey, leading to what are known as ‘wide-angle’ effects in redshift-space distortions. Wide-angle effects will also be present in measurements of the momentum field, i.e. the galaxy density-weighted velocity field, in upcoming peculiar velocity surveys. In this work, we study how wide-angle effects modify the predicted correlation function and power spectrum for momentum statistics, both in autocorrelation and in cross-correlation with the density field. Using both linear theory and the Zel'dovich approximation, we find that deviations from the plane-parallel limit are large and could become important in data analysis for low-redshift surveys. We point out that even multipoles in the cross-correlation between density and momentum are non-zero regardless of the choice of line of sight, and therefore contain new cosmological information that could be exploited. We discuss configuration space, Fourier space, and spherical analyses; providing exact expressions in each case rather than relying on an expansion in small angles. We hope these expressions will be of use in the analysis of upcoming surveys for redshift-space distortions and peculiar velocities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

JET EROSION OF PARTICLE BEDS: PROJECTING CRITICAL SUSPENSION VELOCITIES FROM EFFECTIVE CLEARING / CLEANING RADII

Here we explore the relationship between effective cleaning/clearing radii, ECR, and critical suspension velocities, Ucs. Although one set of physics controls both the radial extent of erosion (i.e., the ECR) and the velocity needed to erode from the center of jet impingement to any given location (i.e., Ucs, typically defined from nozzle center to the vessel center), the relationship between these two remains unexplored quantitatively. Here we advance the model of Kuhn, et al., (PNNL-22816, 2013) as described by Pease, et al. (FEDSM2017-69444, 2017) to evaluate the relationship between the effective clearing/cleaning radius and the nozzle velocity on flat surfaces including flat bottomed vessels. Two governing dimensionless groups are identified. We present both a closed form analytical but transcendental solution and a non-iterative approximation modeled after the Serghides approximation of Colebrook’s nonlinear equation for both ECR versus nozzle velocity and Ucs. Comparison of the model to data from flume testing on a flat surface finds reasonable agreement.

jet, erosion, critical suspension velocity, effect↗
DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Kimberlina 1.2 Velocity Models and Seismic Data

Kimberlina 1.2 Velocity model and synthetic seismic data, produced in collaboration of teams at the National Energy Technology Laboratory, Los Alamos National Laboratory, and Lawrence Livermore National Laboratory through the National Risk Assessment Partnership. Data is associated with the following publication: Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019. The Kimberlina 1.2 Velocity models were produced by Zan Wang, Robert Dilmore, William Harbert, and Lianjie Huang at NETL. The following citations are directly related to the creation of the velocity models: Wang, Z. Harbert, W., Dilmore, R., Huang, L. Modeling of time-lapse seismic monitoring using CO2 leakage simulations for a model CO2 storage site with realistic geology: Application in assessment of early leak-detection capabilities. International Journal of Greenhouse Gas Control. V. 76, September 2018, Pages 39-52. https://doi.org/10.1016/j.ijggc.2018.06.011 Wang, Z., Dilmore, R., Harbert, W. Inferring CO2 saturation from synthetic surface seismic and downhole monitoring data using machine learning for leakage detection at CO2 sequestration sites. International Journal of Greenhouse Gas Control, V. 100, September 2020. https://doi.org/10.1016/j.ijggc.2020.103115 The velocity models were built based on the Kimberlina 1.2 aquifer impact data which is associated with the following publications: Buscheck, T.A., Mansoor, K., Yang, X., Wainwright, H., and Carroll, S. (2019). Downhole pressure and chemical monitoring for CO2 and brine leak detection in aquifers above a CO2 storage reservoir. International Journal of Greenhouse Gas Control. 91. 102812. 10.1016/j.ijggc.2019.102812. Xianjin Yang, Thomas A. Buscheck, Kayyum Mansoor, Zan Wang, Kai Gao, Lianjie Huang, Delphine Appriou, Susan A. Carroll, Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, Volume 90, 2019, 102803, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2019.102803 The synthetic seismic data was produced by Youzuo Lin and team at LANL, and are associated with the following citations: Jordan, P. D., and J. L. Wagoner. Characterizing Construction of Existing Wells to a CO2 Storage Target: The Kimberlina Site, California. Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019.

CO2 Leakage↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A Decade of Radial-velocity Monitoring of Vega and New Limits on the Presence of Planets

We present an analysis of 1524 spectra of Vega spanning 10 yr, in which we search for periodic radial-velocity variations. A signal with a periodicity of 0.676 day and a semi-amplitude of ∼10 m s{sup −1} is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The activity signal appears to evolve on long timescales, which may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega’s photometric rotational modulation for the first time, with a total amplitude of only 10 ppm. A comparison of the spectroscopic and photometric amplitudes suggests that the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial-velocity injection recovery tests exclude the presence of transiting planets larger than 2 R {sub ⊕} and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial-velocity signal with a period of 2.43 days and a semi-amplitude of 6 m s{sup −1}. If caused by an orbiting companion, its minimum mass would be ∼20 M {sub ⊕}; because of Vega’s pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

60 Candidate High-velocity Stars Originating from the Sagittarius Dwarf Spheroidal Galaxy in Gaia EDR3

Using proper motions from Gaia Early Data Release 3 (Gaia EDR3) and radial velocities from several surveys, we identify 60 candidate high-velocity stars with a total velocity greater than 75% of the escape velocity that probably originated from the Sagittarius dwarf spheroidal galaxy (Sgr) by orbital analysis. Sgr’s gravity has little effect on the results and the Large Magellanic Cloud’s gravity has a nonnegligible effect on only a few stars. The closest approach of these stars to the Sgr occurred when the Sgr passed its pericenter (~38.2 Myr ago), which suggests they were tidally stripped from the Sgr. The positions of these stars in the Hertzsprung–Russell diagram and the chemical properties of 19 of them with available [Fe/H] are similar to the Sgr stream member stars. This is consistent with the assumption of their accretion origin. Two of the 60 are hypervelocity stars, which may also be produced by the Hills mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗