SEARCH · Search NASA

Results for “THREE-DIMENSIONAL LATTICES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Effect of Crystal Orientation on Fatigue Failure of Single Crystal Nickel Base Turbine Blade Superalloys

High cycle fatigue (HCF) induced failures in aircraft gas turbine and rocket engine turbopump blades is a pervasive problem. Single crystal nickel turbine blades are being utilized in rocket engine turbopumps and jet engines throughout industry because of their superior creep, stress rupture, melt resistance, and thermomechanical fatigue capabilities over polycrystalline alloys. Currently the most widely used single crystal turbine blade superalloys are PWA 1480/1493, PWA 1484, RENE' N-5 and CMSX-4. These alloys play an important role in commercial, military and space propulsion systems. Single crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the part geometry a significant factor in the overall analysis. The failure modes of single crystal turbine blades are complicated to predict due to the material orthotropy and variations in crystal orientations. Fatigue life estimation of single crystal turbine blades represents an important aspect of durability assessment. It is therefore of practical interest to develop effective fatigue failure criteria for single crystal nickel alloys and to investigate the effects of variation of primary and secondary crystal orientation on fatigue life. A fatigue failure criterion based on the maximum shear stress amplitude /Delta(sub tau)(sub max))] on the 24 octahedral and 6 cube slip systems, is presented for single crystal nickel superalloys (FCC crystal). This criterion reduces the scatter in uniaxial LCF test data considerably for PWA 1493 at 1200 F in air. Additionally, single crystal turbine blades used in the alternate advanced high-pressure fuel turbopump (AHPFTP/AT) are modeled using a large-scale three-dimensional finite element model. This finite element model is capable of accounting for material orthotrophy and variation in primary and secondary crystal orientation. Effects of variation in crystal orientation on blade stress response are studied based on 297 finite element model runs. Fatigue lives at critical points in the blade are computed using finite element stress results and the failure criterion developed. Stress analysis results in the blade attachment region are also presented. Results presented demonstrates that control of secondary and primary crystallographic orientation has the potential to significantly increase a component S resistance to fatigue crack growth with- out adding additional weight or cost. [DOI: 10.1115/1.1413767]

Arakere, N. K.↗

NASA NTRS2004

Morphological changes in woody stem of Prunus jamasakura under simulated microgravity

When the four-week-old woody stem of Prunus jamasakura was grown under simulated microgravity condition on a three-dimensional clinostat, it bent at growth, and width of its secondary xylem decreased due to the reduction of fiber cell numbers and a smaller microfibril angle in the secondary cell wall, as reported in our previous paper. Gravity induces the development of the secondary xylem that supports the stem upward against the action of gravity. In this study, morphological changes of the tissues and cells were microscopically observed. Disorder was found in the concentric structure of tissues that organize the stem. The radial arrangement of the cells was also disturbed in the secondary xylem, and in the secondary phloem secondary cell walls of the bast fiber cells were undeveloped. These findings suggest that differentiation and development of the secondary xylem and the bast fiber cells are strongly controlled by terrestrial gravity. These tissue and cells functions to support the stem under the action of gravity. Furthermore, clinorotation induced disorder in the straight joint of vessel elements and the lattice-like structure of radial parenchyma cells, which is responsible for water transportation and storage, respectively. Gravity is an essential factor for keeping the division and differentiation normal in woody stem.

Prunus/cytology/physiology/ultrastructure↗

DOE OSTI2025

Radiation-resilient InAs quantum dot lasers

Displacement damage from particle radiation increases the threshold current of semiconductor lasers, impeding the deployment of photonic systems in harsh radiation environments. Replacing conventional quantum well (QW) gain regions with quantum dots (QDs) has been shown to suppress the radiation-induced increase in threshold current owing to three-dimensional carrier confinement, which effectively segregates charge carriers away from radiation-induced defects. However, radiation-induced change in threshold current is also affected by extrinsic design elements of a laser, such as photon lifetime, and therefore provides only a partial assessment of the radiation hardness. In this work, we advance the understanding of the radiation resiliency of InAs QD lasers by extracting the change in the nonradiative recombination rate after radiation exposure, which provides a more fundamental assessment of the intrinsic radiation hardness of InAs QD gain regions. We find that carriers that populate InAs QD states have a nonradiative recombination rate that is over an order-of-magnitude less sensitive to neutron radiation relative to carriers in InGaAsP QWs. InAs QD lasers must be designed such that carriers primarily occupy QD states, preventing carrier spillover into surrounding barrier layers to obtain the maximum benefit. Furthermore, we show that enhanced radiation hardness is observed in InAs QD lasers grown on both highly lattice-mismatched silicon substrates and native GaAs substrates, further broadening the appeal of InAs QD lasers as reliable light sources for photonic systems operating in harsh radiation environments.

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Alkali-Metal Interlocking of 2D V 4 O 10 Sheets Defines Discretized Interlayer Shear Relationships

Low-dimensional materials manifest structural anisotropy, quantum confinement, and tightly bound excitonic states, which make them attractive building blocks that can be assembled within three-dimensional laterally stitched heterostructures, stacked van der Waals solids, and complex moiré superlattices. Ion intercalation in the galleries between layered materials provides a means of modifying interlayer separation and coupling, but it is also known to drive the shearing of the layers. In this article, we explore the distinct ligand coordination environments afforded by vanadyl oxygens of singular [V 4 O 10 ] sheets and examine how the size, polarizability, and stoichiometry of Group I cations sandwiched between such layers determine the interlocking of the sheets in stacked structures. Based on the topochemical insertion of alkali-metal ions into the layered λ-V 2 O 5 , we identify seven types of guest ion coordination sites discretized into four distinct regimes of interlayer shear in units of half octahedral widths. The coordination preferences of intercalated cations govern how they interlock 2D [V 4 O 10 ] sheets and engender specific shear conformations. We present evidence that static and dynamic disorder in guest ion arrangement modulate the magnetic structure of the intercalated compounds based on electrostatic polarization, localization of charge and spin density, and lattice distortion. The results illustrate the use of topochemical ion insertion to modulate stacking relationships and magnetic transition characteristics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Tensor renormalization group approach to critical phenomena via symmetry-twisted partition functions

The locality of field theories strongly constrains the possible behaviors of symmetry-twisted partition functions, and thus they serve as order parameters to detect low-energy realizations of global symmetries, such as spontaneous symmetry breaking (SSB). We demonstrate that the tensor renormalization group (TRG) offers an efficient framework to compute the symmetry-twisted partition functions, which enables us to detect the symmetry-breaking transition and also to study associated critical phenomena. As concrete examples of SSB, we investigate the two-dimensional (2D) classical Ising model and the three-dimensional (3D) classical 𝑂⁡(2) nonlinear sigma model, and we identify their critical points solely from the twisted partition function. By employing the finite-size scaling argument, we find the critical temperature 𝑇 𝑐 = 2.2017⁢(2) with the critical exponent 𝜈 = 0.663⁢(33) for the 3D 𝑂⁡(2) model. In addition, we also study the Berezinskii–Kosterlitz–Thouless (BKT) criticality of the 2D classical 𝑂⁡(2) model by extracting the helicity modulus from the twisted partition functions, and we obtain the BKT transition temperature, 𝑇 BKT = 0.8928⁢(2).

lattice field theory↗

DOE OSTI2026

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗

DOE OSTI2025

Scalable mechanochemical synthesis of high-quality Prussian blue analogues for high-energy and durable potassium-ion batteries

Prussian blue analogues (PBAs) are recognized as promising cathode materials for potassium-ion batteries (PIBs), particularly the low-cost and high-energy K 2 Mn[Fe(CN) 6 ](KMnF). However, conventional solution-based synthesis inevitably introduces [Fe(CN) 6 ] 4− defects and lattice water while suffering low synthesis efficiency, unfavorable to the improvement of electrochemical performance and scalability. Here, in this work, we report a simple solvent-free mechanochemical strategy for the synthesis of a wide variety of K 2 M[Fe(CN) 6 ] (M = Mn, Mg, Ca, etc.) with negligible defects and water, and it is unprecedented to achieve kilogram-level products of high-quality KMnF within just 10 minutes. The as-prepared KMnF delivers a high energy density of 590 Wh kg −1 at 0.2 C and exhibits an astonishing stability over 10 000 cycles and rate ability up to 50 C in a potassium metal half-cell. Encouragingly, a high-areal-capacity pouch cell with 2.2 mAh cm −2 (16.5 mg cm −2 ) exhibits a capacity retention of 80.7% after 500 cycles. Furthermore, systematic in situ characterization reveals underlying mechanism insights into structure–performance relationships. Specifically, the fully coordinated Mn–N 6 octahedral configuration effectively suppresses Mn 3+ Jahn–Teller distortion, enabling reversible phase transitions under both high-voltage and long-term cycling conditions. In addition, minimal defects provide sufficient redox centers, while the continuous three-dimensional framework facilitates rapid K + diffusion kinetics. This work provides a new opportunity for the ultrafast, universal and scalable synthesis of high-quality PBAs, facilitating the practical application of PIBs while enabling precise structural and compositional design of novel PBAs.

Mechanochemical method↗

DOE OSTI2025

Revealing rotational symmetry breaking charge density wave order in the kagome superconductor (Rb, K)⁢V 3 ⁢Sb 5 by ultrafast pump-probe experiments

The recently discovered Kagome superconductor 𝐴⁢V 3 ⁢Sb 5 (where 𝐴 refers to K, Rb, Cs) has stimulated widespread research interest due to its interplay of nontrivial topology and unconventional correlated physics including charge-density waves (CDW) and superconductivity. The essential prerequisite to understanding the microscopic mechanisms of this complex electronic landscape is to unveil the configuration and symmetry of the charge-density wave order. As to now, little consensus has been made on what symmetry is broken. Herein, we clarify the microscopic structure and symmetry breaking of the CDW phase in RbV 3 ⁢Sb 5 and KV 3 ⁢Sb 5 by ultrafast time-resolved reflectivity. Our approach is based on extracting coherent phonon spectra induced by three-dimensional CDW and comparing them to calculated phonon frequencies via density-functional theory. The combination of these experimental results and calculations provides compelling evidence that the CDW structure of both compounds prevailing up to 𝑇 CDW is the 2 × 2 × 2 staggered inverse Star-of-David pattern with interlayer 𝜋 phase shift, in which the sixfold rotational symmetry is broken. Finally, these observations thus corroborate sixfold rotational symmetry breaking throughout the CDW phase of RbV 3 ⁢Sb 5 and KV 3⁢ Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Unveiling in-gap states and Majorana zero modes in superconductor-topological insulator bilayer model

Interfaces between topological insulators and superconductors (SCs) are promising platforms for realizing Majorana zero modes (MZMs) via the superconducting proximity effect. We introduce a bilayer model consisting of the surface states of a three-dimensional topological insulator (3DTI) coupled to an 𝑠-wave superconductor and systematically study the role of interlayer tunneling strength (𝑡⊥) motivated by the recent growth of the Fe⁡(Te,Se)/Bi 2 ⁢Te 3 heterostructure. We find that increasing 𝑡⊥ shifts the proximity-induced (PrI) gap minimum away from the Γ point, giving rise to momentum-selective interference patterns that manifest as spatial oscillations in the in-gap states. We introduce an antidot with a magnetic vortex in the model and investigate the resulting in-gap states, including MZMs and Caroli–de Gennes–Matricon (CdGM) modes. With increasing hybridization strength, the energy separation between MZMs and CdGM states increases, enhancing the isolation of MZMs. Importantly, in the strong hybridization limit, the leading CdGM separation remains large in spite of the decrease in the PrI gap. Spin- and spatially resolved wavefunction analysis reveals angular momentum asymmetries absent in conventional 𝑠-wave systems. A direct comparison with a standalone 𝑠-wave superconductor confirms the emergence of distinct 𝑝-wave-like features in the bilayer geometry. Our results provide experimentally relevant predictions for tuning the stability of MZMs and their differentiation from the CdGM modes in SC-3DTI heterostructures and offer a theoretical framework for probing unconventional superconductivity in engineered topological systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

On the Interplay of Electronic and Lattice Screening on Exciton Binding in Two-Dimensional Lead Halide Perovskites

We use path integral Monte Carlo to study the energetics of excitons in layered, hybrid organic–inorganic perovskites in order to elucidate the relative contributions of dielectric confinement and electron–phonon coupling. While the dielectric mismatch between polar perovskite layers and nonpolar ligand layers significantly increases the exciton binding energy relative to their three-dimensional bulk crystal counterparts, formation of exciton polarons attenuates this effect. The contribution from polaron formation is found to be a nonmonotonic function of the lead halide layer thickness, which is clarified by a general variational theory. Accounting for both of these effects provides a description of exciton binding energies in good agreement with experimental measurements. By studying isolated layers and stacked layered crystals of various thicknesses, with ligands of varying polarity, we provide a systematic understanding of the excitonic behavior of this class of materials and how to engineer their photophysics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Redox-Induced Microstructure and Phase Dynamics in Nickel: Insights from In Situ Synchrotron X-ray Diffraction

Using in situ synchrotron X-ray diffraction, we interrogate the microstructural and phase evolution of polycrystalline nickel (Ni) during redox cycling in O 2 , H 2 , and H 2 O environments. Oxidation in O 2 promotes strong (111) texturing in both the NiO overlayer and the underlying Ni substrate. However, this crystallographic alignment is lost following reduction in H 2 and subsequent reoxidation, demonstrating irreversible microstructural changes. H 2 exposure leads to proton dissolution into the Ni lattice, triggering a localized phase transition from face-centered cubic (FCC) to hexagonal close-packed (HCP) Ni in hydrogen-saturated regions. In H 2 O-containing atmospheres, dissociative H 2 O adsorption produces protons that permeate the NiO layer, forming γ-NiOOH within the NiO lattice and HCP Ni beneath the NiO overlayer as protons accumulate. Kinetic analysis via the Johnson-Mehl-Avrami–Kolmogorov model uncovers distinct growth mechanisms: preoxidized Ni surfaces follow one-dimensional (1D) kinetics for NiO, γ-NiOOH, and HCP growth, whereas pristine Ni exhibits three-dimensional (3D) kinetics due to island-like nucleation and growth of NiO. Furthermore, these results highlight the critical interplay between H 2 O dissociation, hydrogen permeation, and redox-driven phase transformations, with practical implications in engineering nickel-based catalysts and hydrogen storage systems through controlled microstructural and phase evolution.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Emergent quasi-one-dimensional antiferromagnetism in the distorted kagome magnet CePtPb

CePtPb hosts a distorted kagome lattice of Ce 3+ ions, providing a clean platform to investigate how reduced local symmetry and strong spin-orbit coupling reshape frustrated magnetism. Magnetization, specific heat, and magnetocaloric effect measurements, combined with a symmetry analysis of the single-ion anisotropy, demonstrate that the local 𝑚⁢2⁢𝑚 site symmetry selects a nearly Ising-like Kramers doublet with easy axes lying within the 𝑎⁢𝑏 plane. This results in three distinct in-plane Ising directions and an overall easy-plane anisotropy. The low-energy magnetic response is well captured by a three-sublattice Ising model, which quantitatively reproduces the saturation magnetization for arbitrary in-plane field orientations, including $[110]$ and $[1\bar{⁢1}⁢0]$, as well as the ratio of the field-induced critical fields. For 𝐵∥$[110]$, the phase diagram exhibits two quantum critical points at 𝐵 c⁢1 = 0.25T and 𝐵 c⁢2 = 0.55T, arising from the sequential polarization of the three Ising sublattices. In conclusion, these results reveal that the system develops quasi-one-dimensional spin chains along the 𝑐 axis, emerging from the nominally three-dimensional crystal structure composed of stacked kagome layers, and illustrate how reduced local symmetry can drive effective dimensional reduction in rare-earth Ising magnets.

Li, Fangli [Southern University of Science and Tec↗

DOE OSTI2025

Quantum Spin Ice in Three-Dimensional Rydberg Atom Arrays

Quantum spin liquids are exotic phases of matter whose low-energy physics is described as the deconfined phase of an emergent gauge theory. With recent theory proposals and an experiment showing preliminary signs of Z 2 topological order [G. Semeghini , ], Rydberg atom arrays have emerged as a promising platform to realize a quantum spin liquid. In this work, we propose a way to realize a U(1) quantum spin liquid in three spatial dimensions, described by the deconfined phase of U(1) gauge theory in a pyrochlore lattice Rydberg atom array. We study the ground state phase diagram of the proposed Rydberg system as a function of experimentally relevant parameters. Within our calculation, we find that by tuning the Rabi frequency, one can access both the confinement-deconfinement transition driven by a proliferation of “magnetic” monopoles and the Higgs transition driven by a proliferation of “electric” charges of the emergent gauge theory. We suggest experimental probes for distinguishing the deconfined phase from ordered phases. This work serves as a proposal to access a confinement-deconfinement transition in three spatial dimensions on a Rydberg-based quantum simulator. Published by the American Physical Society 2025

Shah, Jeet (ORCID:0000000158738129)↗

NASA NTRS1989

Numerical simulations - Some results for the 2- and 3-D Hubbard models and a 2-D electron phonon model

Numerical simulations on the half-filled three-dimensional Hubbard model clearly show the onset of Neel order. Simulations of the two-dimensional electron-phonon Holstein model show the competition between the formation of a Peierls-CDW state and a superconducting state. However, the behavior of the partly filled two-dimensional Hubbard model is more difficult to determine. At half-filling, the antiferromagnetic correlations grow as T is reduced. Doping away from half-filling suppresses these correlations, and it is found that there is a weak attractive pairing interaction in the d-wave channel. However, the strength of the pair field susceptibility is weak at the temperatures and lattice sizes that have been simulated, and the nature of the low-temperature state of the nearly half-filled Hubbard model remains open.

Scalapino, D. J.↗

NASA NTRS1993

Roughness-induced generation of crossflow vortices in three-dimensional boundary layers

The receptivity theory of Goldstein and Ruban is extended within the nonasymptotic (quasi-parallel) framework of Zavol'skii et al to predict the roughness-induced generation of stationary and nonstationary instability waves in three-dimensional, incompressible boundary layers. The influence of acoustic-wave orientation, as well as that of different types of roughness geometries, including isolated roughness elements, periodic arrays, and two-dimensional lattices of compact roughness shapes, as well as random, but spatially homogeneous roughness distributions, is examined. The parametric study for the Falkner-Skan-Cooke family of boundary layers supports our earlier conjecture that the initial amplitudes of roughness-induced stationary vortices are likely to be significantly larger than the amplitudes of similarly induced nonstationary vortices in the presence of acoustic disturbances in the free stream. Maximum unsteady receptivity occurs when the acoustic velocity fluctuation is aligned with the wavenumber vector of the unsteady vortex mode. On the other hand, roughness arrays that are oriented somewhere close to the group velocity direction are likely to produce higher instability amplitudes. Limitations of the nonasymptotic theory are discussed, and future work is suggested.

Choudhari, Meelan↗

NASA NTRS1997

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Tbeodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modem three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

NASA NTRS1997

Identification of Linear and Nonlinear Aerodynamic Impulse Responses Using Digital Filter Techniques

This paper discusses the mathematical existence and the numerically-correct identification of linear and nonlinear aerodynamic impulse response functions. Differences between continuous-time and discrete-time system theories, which permit the identification and efficient use of these functions, will be detailed. Important input/output definitions and the concept of linear and nonlinear systems with memory will also be discussed. It will be shown that indicial (step or steady) responses (such as Wagner's function), forced harmonic responses (such as Theodorsen's function or those from doublet lattice theory), and responses to random inputs (such as gusts) can all be obtained from an aerodynamic impulse response function. This paper establishes the aerodynamic impulse response function as the most fundamental, and, therefore, the most computationally efficient, aerodynamic function that can be extracted from any given discrete-time, aerodynamic system. The results presented in this paper help to unify the understanding of classical two-dimensional continuous-time theories with modern three-dimensional, discrete-time theories. First, the method is applied to the nonlinear viscous Burger's equation as an example. Next the method is applied to a three-dimensional aeroelastic model using the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code and then to a two-dimensional model using the CFL3D Navier-Stokes code. Comparisons of accuracy and computational cost savings are presented. Because of its mathematical generality, an important attribute of this methodology is that it is applicable to a wide range of nonlinear, discrete-time problems.

Silva, Walter A.↗

NASA NTRS2002

Force Evaluation in the Lattice Boltzmann Method Involving Curved Geometry

The present work investigates two approaches for force evaluation in the lattice Boltzmann equation: the momentum- exchange method and the stress-integration method on the surface of a body. The boundary condition for the particle distribution functions on curved geometries is handled with second order accuracy based on our recent works. The stress-integration method is computationally laborious for two-dimensional flows and in general difficult to implement for three-dimensional flows, while the momentum-exchange method is reliable, accurate, and easy to implement for both two-dimensional and three-dimensional flows. Several test cases are selected to evaluate the present methods, including: (i) two-dimensional pressure-driven channel flow; (ii) two-dimensional uniform flow past a column of cylinders; (iii) two-dimensional flow past a cylinder asymmetrically placed in a channel (with vortex shedding); (iv) three-dimensional pressure-driven flow in a circular pipe; and (v) three-dimensional flow past a sphere. The drag evaluated by using the momentum-exchange method agrees well with the exact or other published results.

Mei, Renwei↗