Search NASA⌕ Search

SEARCH · Search NASA

Results for “force method”

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 1,009 records · Page 56

COSMIC-RAY CURRENT-DRIVEN TURBULENCE AND MEAN-FIELD DYNAMO EFFECT

We show that an {alpha} effect is driven by the cosmic-ray (CR) Bell instability exciting left-right asymmetric turbulence. Alfven waves of a preferred polarization have maximally helical motion, because the transverse motion of each mode is parallel to its curl. We show how large-scale Alfven modes, when rendered unstable by CR streaming, can create new net flux over any finite region, in the direction of the original large-scale field. We perform direct numerical simulations (DNSs) of a magnetohydrodynamic fluid with a forced CR current and use the test-field method to determine the {alpha} effect and the turbulent magnetic diffusivity. As follows from DNS, the dynamics of the instability has the following stages: (1) in the early stage, the small-scale Bell instability that results in the production of small-scale turbulence is excited; (2) in the intermediate stage, there is formation of larger-scale magnetic structures; (3) finally, quasi-stationary large-scale turbulence is formed at a growth rate that is comparable to that expected from the dynamo instability, but its amplitude over much longer timescales remains unclear. The results of DNS are in good agreement with the theoretical estimates. It is suggested that this dynamo is what gives weakly magnetized relativistic shocks such as those from gamma-ray bursts (GRBs) a macroscopic correlation length. It may also be important for large-scale magnetic field amplification associated with CR production and diffusive shock acceleration in supernova remnants (SNRs) and blast waves from GRBs. Magnetic field amplification by Bell turbulence in SNRs is found to be significant, but it is limited owing to the finite time available to the super-Alfvenicly expanding remnant. The effectiveness of the mechanisms is shown to be dependent on the shock velocity. Limits on magnetic field growth in longer-lived systems, such as the Galaxy and unconfined intergalactic CRs, are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of active region coronal magnetic field in determining coronal mass ejection propagation direction

We study the role of the coronal magnetic field configuration of an active region (AR) in determining the propagation direction of a coronal mass ejection (CME). The CME occurred in the AR 11944 (S09W01) near the disk center on 2014 January 7 and was associated with an X1.2 flare. A new CME reconstruction procedure based on a polarimetric technique is adopted, which shows that the CME changed its propagation direction by around 28° in latitude within 2.5 R{sub ⊙} and 43° in longitude within 6.5 R{sub ⊙} with respect to the CME source region. This significant non-radial motion is consistent with the finding of Möstl et al. We use nonlinear force-free field and potential field source surface extrapolation methods to determine the configurations of the coronal magnetic field. We also calculate the magnetic energy density distributions at different heights based on the extrapolations. Our results show that the AR coronal magnetic field has a strong influence on the CME propagation direction. This is consistent with the “channeling” by the AR coronal magnetic field itself, rather than deflection by nearby structures. These results indicate that the AR coronal magnetic field configuration has to be taken into account in order to determine CME propagation direction correctly.

79 ASTRONOMY AND ASTROPHYSICS↗

A study of electroactive polyvinyl chloride (PVC) gel actuators through the use of the electric modulus formalism and cyclic linear voltage sweeps

The conductivity and dielectric properties are integral to the function of polyvinyl chloride (PVC) gel actuators. The frequency-dependent properties of PVC gel actuators are investigated here in terms of their impedance, permittivity, and for the first time, electric modulus. The data shows that PVC gels’ conductive properties are just as, if not more, important as their dielectric properties in electromechanical transduction applications. The electrode polarization (EP) and its impact on the impedance and dielectric spectra of PVC gels, as well as the developed asymmetric space charge at the anode, are discussed. The electric modulus and tan δ spectra are used for the fitting of Cole-Cole (CC) and Debye relaxation models for gels of varying plasticizer content. The electrostatic adhesive force for PVC gels of varying plasticizer is also measured, indicating large electrostatic adhesion (>2 N cm -2 ). A cyclic linear voltage sweep is used to clarify the dynamics of space charge within the gels. The peak current (associated with space charge development) is seen to be concurrent with the onset of mechanical deformation, showing the asymmetric charge as the origin of electromechanical transduction. Additionally, the maximum charge transferred (as measured by the integration of current over time) before space charge development is found to correlate with the electrostatic adhesive force measured for the gels, pointing to a new method of characterizing PVC gels for actuation.

Instruments & Instrumentation↗

Halo clustering and g{sub NL}-type primordial non-gaussianity

A wide range of multifield inflationary models generate non-Gaussian initial conditions in which the initial adiabatic fluctuation is of the form (ζ{sub G}+g{sub NL}ζ{sub G}{sup 3}). We study halo clustering in these models using two different analytic methods: the peak-background split framework, and brute force calculation in a barrier crossing model, obtaining agreement between the two. We find a simple, theoretically motivated expression for halo bias which agrees with N-body simulations and can be used to constrain g{sub NL} from observations. We discuss practical caveats to constraining g{sub NL} using only observable properties of a tracer population, and argue that constraints obtained from populations whose observed bias is ∼<2.5 are generally not robust to uncertainties in modeling the halo occupation distribution of the population.

79 ASTRONOMY AND ASTROPHYSICS↗

Semi-microscopic model of pairing in nuclei

A semi-microscopic model for nucleon pairing in nuclei is presented starting from the ab initio BCS gap equation with the Argonne v{sub 18} force and the self-consistent energy density functional method basis characterized by the bare nucleon mass. The BCS theory is formulated in terms of the model space S{sub 0} with the effective pairing interaction calculated from the first principles in the subsidiary space S{sup '}. This effective interaction is supplemented with a small phenomenological addendum containing one phenomenological parameter, universal for all medium and heavy atomic nuclei. We consider the latter as a phenomenological way to take into account approximately both the many-body corrections to the BCS theory and the effective mass effects. For protons, the Coulomb interaction is introduced directly. Calculations made for several isotopic and isotonic chains of semi-magic nuclei are compared with empirical gap values derived from odd-even mass differences. The average disagreement is of the order of 0.1--0.2 MeV. The role of the self-consistent basis is analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal in-line force on a cylinder suspended at a fixed position were available from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to computational fluid dynamics (CFD) simulations based on the solution of Reynolds-averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM+. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM+ based on wave modeling, boundary conditions, turbulence modeling, and spatial and temporal discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order potential flow and viscous drag from Morison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Substantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗

Microscopy Techniques for Investigating the Control of Organic Constituents on Biomineralization

In this paper we address recent advances in the application of microscopy techniques to characterize crystallization processes as they relate to biomineralization and bio-inspired materials synthesis. In particular, we focus on studies aimed at revealing the role that organic macromolecules and functionalized surfaces have in modulating the mechanisms of nucleation and growth. In the nucleation studies, we explore the use of methods such asin situ transmission electron microscopy (TEM),atomic force microscopy (AFM), and cryogenic electron microscopy (cryoEM)to delineate formation pathways, phase stabilization, and the competing effects of free energy and kinetic barriers. In the growth studies, an emphasis is placed on understanding the interactions of macromolecular constituents with growing crystals and how to characterize the internal structures of the resulting composite crystals using techniques such as electron tomography, atom probe tomography (APT), and vibrational spectromicroscopy. Examples are drawn from both biological and bio-inspired synthetic systems.

Hendley, Coit T.↗

Laser-Induced Ultrafast Demagnetization and Spin Switching (Final Report)

The information technology revolution demands bigger and faster magnetic storage devices. The funded research focuses on a new frontier - Laser-induced ultrafast all-optical spin switching (AOS), where a single laser pulse can switch spins from one direction to another permanently, without a magnetic field, on a time scale of 1-10 picoseconds, 2-3 orders of magnitude faster than the existing technology. The funded research has successfully resolved a long-standing issue in ultrafast demagnetization whether the lattice contributes to the demagnetization. We developed a high level first-principles theory and employed a massively parallel computer to pin down the details of the physical mechanisms that are inaccessible to other methods. This unique feature allows us to compute the force on each atom in both bulk and monolayer systems. We find that it is indeed true there is a nonzero force on each atom, but the idea of vibration contribution is not crucial to our basic understanding of physics. This paves the way to future experimental, theoretical and technological developments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Microchannel-based Membrane-less Extraction of Li from Unconventional Lithium Sources & the Separation of REE

This final report provides an overview of the Project's entire duration, covering July 1, 2021 to December 31, 2023. It primarily focuses on the achievements, technological developments, and unique challenges the team faced while working on separating and extracting Lithium from produced waters. The project's primary aim was to create an integrated, high-throughput, membrane-less, and modular microfluidic platform that could extract Lithium from unconventional sources. We have successfully met all goals and milestones envisioned in the SOPO document. The most critical primary milestones, including the Go-No-Go milestone (refer to the Gantt chart in the Appendices), were successfully accomplished. We demonstrated phase separation (>90%) and extraction (>85%) performance in the MPSE using synthetic, and representative produced water composition feed at 50 ml/min total flow through MPSE 36. We have also performed a parametric study of the MPSE operations, beyond the scope of SOPO, exploring operating conditions of current and broader interest. The extended investigation of operational parameters is concurrent with our efforts to seek further development of the MPSE technology beyond the scope of the Project. Along these lines of development, we have made efforts to be responsive to DOE calls for technological developments of other types of resources (beyond PW) for the recovery of Critical Materials and higher TRL development (beyond TRL 4). During the work on this Project, we developed and implemented three innovative technical approaches that emerged from our efforts to successfully meet the Project milestones. The innovative & original technical approaches developed and implemented in this Project are now the contributions to process engineering that could be clearly credited to the Project. First, Convergent Design Approach is a comprehensive feedforward & feedback loop of four design phases: i) design for functionality, ii) design for manufacturing, iii) design for sustainability, and iv) design for market. Next was Process Intensification. A major aim of this Project was to create an innovative phase separation & extraction microscale-based technology for Li separation – thus the words microchannel-based in the Project title. A microscale-based technology is intrinsically in the center of the Process Intensification domain as defined by its unique principles. Therefore, Process Intensification was implicitly envisioned in the Project’s SOPO. Lastly, Time Scale Analysis is a novel tool for discovering the needs and directions of Process Intensification implementations in any process technology. This Project is fully credited for developing and implementing the three novel technical approaches mentioned above. These are general contributions to process engineering that emerged from this Project. Beyond the original SOPO scope, the OSU-U.Pitt research group utilized a Convergent Design methodology, integrating first-principles mathematical modeling with experimental validation on the Minimum Development Vehicle. By creating these Digital Twins, the team rapidly assessed manufacturing iterations to support TEA analysis. This framework further enabled the development of advanced Surface Modification Techniques, where hydrophobic and oleophobic coating strategies were optimized via Digital Twin tools and validated through rigorous 100-hour longevity testing. TEA Analysis: The closing efforts of this Project were focused on the TEA analysis. TEA analysis had two primary functions: i) enabling critical assessments of design variations withing 10 the Concurrent Design Approach, thus enabling evolution of the MPSE design to reach faster- better-cheaper alternatives; and ii) to create a bridge between the accomplishments of this Project and future projects of higher TRL, beyond TRL 6 level. It is important to note that the TEA model created in the Project stirred the technological solutions for the recovery of critical materials toward a vision of a very profitable modular plant that has unique zero-waste water discharge signature. More importantly, thanks to our experimental performance data and conservative assumptions, the TEA model predicts minimal technological and investment risks. Low cost of a modular unit of a nominal capacity of [1000 tons of Li 2 CO 3 /year] positions the MPSE based technology within the reach of community investors, thus offering a paradigm shift in the development of critical technologies. The project successfully navigated two primary challenges: solvent selection and manufacturing adaptation. Restricted by the SOPO to existing literature for lithium recovery, the team identified a critical need for a "material excellence program" to develop next-generation solvents, eventually concluding with a preliminary investigation into promising Ionic Liquids (ILs). Simultaneously, COVID-19 supply chain disruptions forced a pivot from traditional manufacturing to advanced additive methods at ATAMI-OSU. By transitioning from stainless steel to 3D-printed polymer substrates, the team achieved a transformative three-order-of- magnitude reduction in manufacturing costs and compressed prototyping timelines from several months to just two days. The MPSE technology offers significant energy, environmental, and economic advantages by overcoming the traditional bottlenecks of phase-separation hardware and contactor size. Unlike conventional mixer-settlers or membrane-based systems, MPSE operates without moving parts or fouling-prone membranes, achieving robust performance even with challenging, viscous, or particulate-heavy feeds. Key performance metrics include an energy intensity reduction of 5–50x (3–40 kJ/m 3 ) compared to incumbent technologies and a dramatic reduction of processing time to under 60 seconds, which drastically reduces the physical plant footprint. These technical efficiencies translate into superior economic outcomes; for a 100 t/year Li 2 CO 3 facility, implementing MPSE is projected to nearly halve contactor CAPEX (from $\$$6.08M to $\$$3.01M) and significantly increase the project's Net Present Value (NPV), derisking new investment and enabling distributed critical-mineral processing configurations. The commercialization of MPSE technology is being spearheaded by Vigsur Dynamics Inc., which has adopted a structured, parallel approach to technical and business development since its formation in January 2026. Following extensive customer discovery and engagement with the Oregon State University accelerator, Vigsur Dynamics is working to establish a business model that transitions from pilot demonstrations to modular hardware sales, ultimately aiming for a "build-own-operate" service strategy. Current technical milestones—including 100 hours of continuous operation, superior energy efficiency, and successful 6-unit modular scale-up— provide a foundation for this transition. Backed by ongoing IP licensing and a growing network of industrial and venture advisors, the company is actively de-risking the platform to replace conventional mixer-settler systems in the critical minerals market.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Soft actuators and soft actuating devices

A soft buckling linear actuator is described, including: a plurality of substantially parallel bucklable, elastic structural components each having its longest dimension along a first axis; and a plurality of secondary structural components each disposed between and bridging two adjacent bucklable, elastic structural components; wherein every two adjacent bucklable, elastic structural components and the secondary structural components in-between define a layer comprising a plurality of cells each capable of being connected with a fluid inflation or deflation source; the secondary structural components from two adjacent layers are not aligned along a second axis perpendicular to the first axis; and the secondary structural components are configured not to buckle, the bucklable, elastic structural components are configured to buckle along the second axis to generate a linear force, upon the inflation or deflation of the cells. Methods of actuation using the same are also described.

42 ENGINEERING↗

Soft actuators and soft actuating devices

A soft buckling linear actuator is described, including: a plurality of substantially parallel bucklable, elastic structural components each having its longest dimension along a first axis; and a plurality of secondary structural components each disposed between and bridging two adjacent bucklable, elastic structural components; wherein every two adjacent bucklable, elastic structural components and the secondary structural components in-between define a layer comprising a plurality of cells each capable of being connected with a fluid inflation or deflation source; the secondary structural components from two adjacent layers are not aligned along a second axis perpendicular to the first axis; and the secondary structural components are configured not to buckle, the bucklable, elastic structural components are configured to buckle along the second axis to generate a linear force, upon the inflation or deflation of the cells. Methods of actuation using the same are also described.

Yang, Dian↗

In-situ gloveport glove leak tester

The invention provides a system for detecting leaks in glovebox gloves, the system having a first seal within a glovebox aperture; a second seal contacting exterior surfaces of a glovebox, wherein the exterior surfaces define a periphery of the aperture; a first fluid supply for circumferentially expanding the first seal against center facing surfaces of the aperture to a first pressure; a second fluid supply for inflating the glove to a second pressure; and a pressure gauge for detecting a pressure decrease of the second pressure. Also provided is an in situ method for detecting leaks in glove box gloves, the method having the steps of applying a radially directed force to inside surfaces of cuffs of the gloves to hermetically seal the gloves from the exterior of the glovebox, thereby establishing a seal; inflating the gloves to a pressure; and monitoring the pressure for any decrease in pressure.

42 ENGINEERING↗

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code: Preprint

A passive yaw implementation is developed, validated, and explored for the Wave Energy Converter Simulator (WEC-Sim), an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01 was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than or equal to 1, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗

Model Structure and Climate Data Uncertainty in Historical Simulations of the Terrestrial Carbon Cycle (1850–2014)

Abstract The divergence among Earth system models in the terrestrial carbon cycle has prompted interest in how to reduce uncertainty. Previous studies have identified model structural uncertainty arising from process parameterizations and parameter values. The current study highlights the importance of climate forcing in generating carbon cycle uncertainty. We use simulations in which three models (Community Land Model version 4 (CLM4), CLM4.5, CLM5) with substantially different carbon cycles are forced with two climate reconstructions (CRUNCEPv7, Global Soil Wetness Project 3 version 1 (GSWP3v1)) to examine the contributions of model structure and climate to uncertainty in the carbon cycle over the period 1850–2014. Climate uncertainty for global annual net biome production exceeds one third of total uncertainty (defined as the sum of climate and model structure uncertainty) in the first half of the twentieth century, but declines after the 1950s. Global annual gross primary productivity, net primary productivity, heterotrophic respiration, and vegetation and soil carbon stocks have substantial climate uncertainty (relative to total uncertainty) throughout the simulation period. Climate forcing contributes more than one half of total uncertainty for these carbon cycle fluxes and stocks throughout boreal North America and Eurasia, some midlatitude regions, and in eastern Amazonia and western equatorial Africa during the decade 2000–2009. Comparison with observationally based data sets of the carbon cycle using model benchmarking methods provides insight into strengths and deficiencies among models and climate forcings, but we caution against overreliance on benchmarking to discriminate among models. The conceptualization of uncertainty arising from this study implies embracing multiple feasible model simulations rather than focusing on which model or simulation is best.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrical properties of films of zinc oxide nanoparticles and its hybrid with reduced graphene oxide

Free-standing films of ZnO nanoparticles (NPs) and reduced graphene oxide (rGO)-ZnO NPs hybrid are prepared at a liquid/liquid interface. The films are characterized by UV-visible spectroscopy, X-ray diffraction, scanning electron microscopy and atomic force microscopy. ZnO film consists of spherical aggregated NPs while the hybrid film contains folded sheets of rGO with embedded ZnO NPs. Electrical properties of the films and its photoresponse in presence of UV radiation are investigated using current sensing atomic force microscopy (CSAFM) at nanoscale and bulk measurements using two probe methods. Enhancement in photocurrent is observed in both cases and the current imaging reveals an inhomogeneous contribution by different ZnO grains in the film.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Quantifying sensitivity in numerical weather prediction-modeled offshore wind speeds through an ensemble modeling approach

A decade of research has shown that numerical weather prediction (NWP)-modeled wind speeds can be highly sensitive to the inputs and setups within the NWP model. For wind resource characterization applications, this sensitivity is often addressed by constructing a range of setups and selecting the one that best validates against observations. However, this approach is not possible in areas that lack high-quality hub height observations, especially offshore wind areas. In such cases, techniques to quantify and disseminate confidence in NWP-modeled wind speeds in the absence of observations are needed. We address this need in the present study and propose best practices for quantifying the spread in NWP-modeled wind speeds. We implement an ensemble approach in which we consider 24 different setups to the Weather Research and Forecasting (WRF) model. We construct the ensemble by considering variations in WRF version, WRF namelist, atmospheric forcing, and sea surface temperature (SST) forcing. Our analysis finds that the standard deviation produces more consistent estimates compared to the interquartile range and tends to be the more conservative estimator for ensemble variability. We further find that model spread increases closer to the surface and on shorter time scales. In conclusion, we explore methods to attribute total ensemble variability to the different ensemble components (e.g., atmospheric forcing and SST product) and find that contributions by components also vary depending on time scale. We anticipate that the methods and results presented in this paper will provide a reasonable foundation for further research into ensemble-based wind resource data sets.

17 WIND ENERGY↗

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗