2025-06-18 | | Total: 99
This paper presents a social network analysis of the professional support networks of 100 LGBTQ+ and/or women PhD physicists, comparing the networks based on the career sectors of academia, industry, and government/nonprofit. The methods for constructing and analyzing the ego networks, which are novel in many ways, are explained in greater detail in an earlier publication (Hatcher et. al., 2025). We use statistical tests of independence to explore differences between sectors in terms of whole network metrics, network composition based on alter characteristics, and support types. We find that alters associated with groups (like affinity groups and personal and professional interest groups) are more likely to provide identity-based and community building support, participants in Academia have fewer personal friends in their networks while those in Industry have more, participants in Government report less instrumental support, and those in Academia report less material support. These results and others lead to suggestions for employers in these sectors on how to better support these physicists, including continuing to promote participation in affinity and interest groups, providing more material support and/or personal time in the academic sector, and more instrumental support in the form of professional development or training in the government sector.
We use a quantum-logic spectroscopy (QLS) protocol to control the quantum state of a CaH+ ion in a cryogenic environment, in which reduced thermal radiation extends rotational state lifetimes by an order of magnitude over those at room temperature. By repeatedly and adaptively probing the molecule, detecting the outcome of each probe via an atomic ion, and using a Bayesian update scheme to quantify confidence in the molecular state, we demonstrate state preparation and measurement (SPAM) in a single quantum state with infidelity less than $6\times10^{-3}$ and measure Rabi flopping between two states with greater than 99% contrast. The protocol does not require any molecule-specific lasers and the detection scheme is non-destructive. To our knowledge, this result represents the highest SPAM fidelity of a single molecular quantum state demonstrated to date.
The air flows in the proximal and distal portions of the human lungs are interconnected: the lower Reynolds number in the deeper generations causes a progressive flow regularization, while mass conservation requires flow rate oscillations to propagate through the airway bifurcations. To explain how these two competing effects shape the flow state in the deeper generations, we have performed the first high-fidelity numerical simulations of the air flow in a lung model including 23 successive bifurcations of a single planar airway. Turbulence modelling or assumptions on flow regimes are not required. The chosen flow rate is stationary (steady on average), and representative of the peak inspiratory flow reached by adult patients breathing through therapeutical inhalers. As expected, advection becomes progressively less important after each bifurcation, until a time-dependent Stokes regime governed solely by viscous diffusion is established in the smallest generations. However, fluctuations in this regime are relatively fast and large with respect to the mean flow, which is in contrast with the commonly agreed picture that only the breathing frequency is relevant at the scale of the alveoli. We demonstrate that the characteristic frequency and amplitude of these fluctuations are linked to the flow in the upper part of the bronchial tree, as they originate from the time-dependent flow splitting in the upper bifurcations. Even though these fluctuations are observed here in an idealized, rigid lung model, our findings suggest that the assumptions usually adopted in many of the current lung models might need to be revised.
The features of delayed desorption from CH4-doped Ar matrices irradiated with an electron beam of subthreshold energy were studied. Radiolysis products were detected by emission spectroscopy. The total desorption yield was monitored by recording the pressure in the experimental chamber. Based on the analysis of the concentration dependence of delayed desorption bursts and their structure, an assumption was made about the formation of CH4 clusters in Ar matrices. At a high dopant concentration of 10%, up to three consecutive bursts were recorded. Delayed desorption from a sample doped with 0.1% CH4 was registered for the first time. A correlation has been found between the burst of particles and the flash of luminescence of H atoms. This seemingly contradictory observation of the correlation of the H atom content with the particle explosion due to their recombination was explained by analyzing the energy transfer and capture processes and the features of the energy structure of H2 molecule. A linear dependence of the total particle yield on the CH4 concentration in the range of 1-10% was found. The dynamics of intensity changes in the sequence of main flashes and their delay time were discussed.
We investigate the triangular instability of a Batchelor vortex subjected to a stationary triangular strain field generated by three satellite vortices, in the presence of weak axial flow. The analysis combines theoretical predictions with numerical simulations. Theoretically, the instability arises from resonant coupling between two quasi-neutral Kelvin modes with azimuthal wavenumbers $m$ and $m+3$ with the background strain. Numerically, we solve the linearized Navier-Stokes equations around a quasi-steady base flow to identify the most unstable modes, and compare their growth rates and frequencies with theoretical predictions for a Reynolds number $Re = 10^4$ and a straining strength $\epsilon = 0.008$. In the absence of axial flow, only the mode pair $(m_A, m_B) = (-1,2)$ (and its symmetric counterpart) is unstable. However, we show that additional combinations such as $(0,3)$, $(1,4)$, and $(2,5)$, which are otherwise strongly damped by the critical layer in the absence of axial flow, also become unstable once axial flow exceeds a certain threshold, as the critical layer damping is significantly reduced. Furthermore, we show that the most unstable mode in the no-axial-flow case, originating from the second branch of $m = -1$ and the first branch of $m = 2$, becomes less unstable as axial flow increases. It is eventually overtaken by a mode from the first branches of both wavenumbers, which then remains the dominant unstable mode across a wide range of axial flow strengths and Reynolds numbers. A comprehensive instability diagram as a function of the axial flow parameter is presented.
Processing in memory (PIM) has received significant attention due to its high efficiency, low latency, and parallelism. In optical computation, coherent memory is a crucial infrastructure for PIM frameworks. This study presents an all-optical convolution experiment conducted within computational storage based on a cold atomic ensemble. By exploiting the light-atom phase transfer facilitated by the electromagnetically induced transparency, we demonstrated spiral phase contrast processing of photon images in memory, resulting in the edge enhancement of retrieved images recorded using time-correlated photon imaging. In particular, adopting state-of-the-art atomic techniques provides a coherent memory lifetime exceeding 320 us for PIM operations. Our results highlight the significant potential of cold atomic ensembles as computational storage for developing all-optical PIM systems.
Underwater explosions produce complex fluid phenomena relevant to diverse applications including maritime engineering, medical therapeutics, and inertial confinement fusion. These systems exhibit multiphase flows, chemical kinetics, and highly compressible dynamics that challenge traditional computational approaches. Current hydrodynamic solvers, while accurate, are computationally expensive and non-differentiable, limiting their use in design optimization and real-time applications. Here we show that deep neural networks can predict underwater explosion-induced fluid dynamics 4,025 times faster than traditional solvers while maintaining mean absolute percent errors below 0.005\% across all fluid state variables. Our approach maps from explosive material thermodynamic parameters to the temporal evolution of shock fronts and material interfaces, enabling rapid prediction of system behavior for a broad range of ideal explosive materials. Feature importance analysis reveals that exponential decay parameters of the explosive equation of state are the primary drivers of system dynamics, uncovering a previously unknown relationship between thermodynamic compatibility and energy transfer efficiency at material interfaces. Furthermore, we demonstrate an inverse design framework that leverages the differentiability of our neural surrogate to perform parameter discovery, recovering unknown explosive material properties to within 1\% accuracy through gradient-based optimization. This combination of rapid inference, physical insight, and inverse design capabilities provides a route to engineering controlled fluid behavior in underwater explosive systems through material design. We anticipate our approach could enable new applications in defense systems, underwater manufacturing, and medical procedures where precise control of shock waves and bubble dynamics is essential.
This study aims to develop and evaluate a digital twin (DT) framework to enhance adaptive proton therapy for prostate stereotactic body radiotherapy (SBRT), focusing on improving treatment precision for dominant intraprostatic lesions (DILs) while minimizing organ-at-risk (OAR) toxicity. We propose a decision-theoretic (DT) framework combining deep learning (DL)-based deformable image registration (DIR) with a prior treatment database to generate synthetic CTs (sCTs) for predicting interfractional anatomical changes. Using daily CBCT from five prostate SBRT patients with DILs, the framework precomputes multiple plans with high (DT-H) and low (DT-L) similarity sCTs. Plan optimization is performed in RayStation 2023B, assuming a constant RBE of 1.1 and robustly accounting for positional and range uncertainties. Plan quality is evaluated via a modified ProKnow score across two fractions, with reoptimization limited to 10 minutes. Daily CBCT evaluation showed clinical plans often violated OAR constraints (e.g., bladder V20.8Gy, rectum V23Gy), with DIL V100 < 90% in 2 patients, indicating SIFB failure. DT-H plans, using high-similarity sCTs, achieved better or comparable DIL/CTV coverage and lower OAR doses, with reoptimization completed within 10 min (e.g., DT-H-REopt-A score: 154.3-165.9). DT-L plans showed variable outcomes; lower similarity correlated with reduced DIL coverage (e.g., Patient 4: 84.7%). DT-H consistently outperformed clinical plans within time limits, while extended optimization brought DT-L and clinical plans closer to DT-H quality. This DT framework enables rapid, personalized adaptive proton therapy, improving DIL targeting and reducing toxicity. By addressing geometric uncertainties, it supports outcome gains in ultra-hypofractionated prostate RT and lays groundwork for future multimodal anatomical prediction.
Density Functional Theory (DFT) is the most widely used electronic structure method for predicting the properties of molecules and materials. Although DFT is, in principle, an exact reformulation of the Schrödinger equation, practical applications rely on approximations to the unknown exchange-correlation (XC) functional. Most existing XC functionals are constructed using a limited set of increasingly complex, hand-crafted features that improve accuracy at the expense of computational efficiency. Yet, no current approximation achieves the accuracy and generality for predictive modeling of laboratory experiments at chemical accuracy -- typically defined as errors below 1 kcal/mol. In this work, we present Skala, a modern deep learning-based XC functional that bypasses expensive hand-designed features by learning representations directly from data. Skala achieves chemical accuracy for atomization energies of small molecules while retaining the computational efficiency typical of semi-local DFT. This performance is enabled by training on an unprecedented volume of high-accuracy reference data generated using computationally intensive wavefunction-based methods. Notably, Skala systematically improves with additional training data covering diverse chemistry. By incorporating a modest amount of additional high-accuracy data tailored to chemistry beyond atomization energies, Skala achieves accuracy competitive with the best-performing hybrid functionals across general main group chemistry, at the cost of semi-local DFT. As the training dataset continues to expand, Skala is poised to further enhance the predictive power of first-principles simulations.
Polarized solid targets produced via Dynamic Nuclear Polarization rely on Continuous-Wave Nuclear Magnetism Resonance measurements to accurately determine the degree of polarization of bulk samples polarized to nearly 100%. Since the late 1970's phase sensitive detection methods have been utilized to observe the magnetization of a sample as a small change in inductance under RF excitation near the Larmor frequency of the nuclear species of interest, using a device known as a Q-meter. Liverpool Q-meters, produced in the UK in the 80's and 90's, have been the workhorse devices for these targets for decades, however their age and scarcity has meant new systems are needed. We describe a Q-meter system designed and built at Jefferson Lab in the Liverpool style to have comparable electronic performance with several improvements to update and adapt the devices for modern use.
Many biological materials must couple mechanical strength with the ability to rapidly self-assemble at a specific location. In particular, biomolecular condensates readily self-assemble via phase separation, but may also need to resist external forces to fulfil their function. Spatial localisation of condensate formation can be controlled by active cores that preferentially drive the production of condensate material at a particular point, while resistance to external forces can be facilitated by viscoelastic material properties. To investigate the interplay of these two processes, we develop a continuum model of viscoelastic growth around an active core. We find that viscoelastic stresses restrict condensate growth, but also impart resistance to deformation. We investigate the effect of different incorporation schemes on growth dynamics, and test the influence of mechanical properties on condensate strength. Finally, we compare the predictions of our model to experimental data from centrosomes in C. elegans embryos, identifying a parameter regime in which rapid growth can be combined with appropriate mechanical strength, and studying how strain-dependent material incorporation may lead to isotropic growth of scaffold material. Our results provide general design principles for other materials that must reconcile rapid, localised self-assembly with mechanical strength, such as focal adhesions.
We present a quantum-assisted generative algorithm for synthetic tracks of Lagrangian tracer particles in a turbulent shear flow. The parallelism and sampling properties of quantum algorithms are used to build and optimize a parametric quantum circuit, which generates a quantum state that corresponds to the joint probability density function of the classical turbulent velocity components, p(u_1^{\prime}, u_2^{\prime}, u_3^{\prime}). Velocity samples are drawn by one-shot measurements on the quantum circuit. The hybrid quantum-classical algorithm is validated with two classical methods, a standard stochastic Lagrangian model and a classical sampling scheme in the form of a Markov-chain Monte Carlo approach. We consider a homogeneous turbulent shear flow with a constant shear rate S as a proof of concept for which the velocity fluctuations are Gaussian. The generation of the joint probability density function is also tested on a real quantum device, the 20-qubit IQM Resonance quantum computing platform for cases of up to 10 qubits. Our study paves the way to applications of Lagrangian small-scale parameterizations of turbulent transport in complex turbulent flows by quantum computers.
Integrated photonic neural networks (PNNs) have demonstrated significant potential to complement the digital electronic counterparts [1-3]. Nevertheless, robust and repeatable performance of scalable integrated PNNs is directly tied to the quality of their training. Error backpropagation (BP), which relies on nonlinear activation gradient computation, is the mainstream algorithm to train digital neural networks due to its scalability, versatility, and implementation efficiency [4]. Consequently, its adoption is highly desirable for the training of scalable PNNs. Despite such benefits and due to the lack of scalable on-chip activation gradient [5], PNNs have mostly been trained using a digital computer to run BP, which is inadequate in addressing device variations, or through gradient-free algorithms that do not fully benefit from the versatility of BP training. Here, we report the demonstration of an integrated photonic deep neural network with end-to-end on-chip gradient-descent BP training. All linear and nonlinear computations are performed on a single photonic chip, leading to scalable and robust training despite the considerable--but typical--fabrication-induced device variations. Two nonlinear data classification tasks are demonstrated in which the chip performance matches that of the ideal digital model, both in accuracy and robustness. Integrating the advantages of BP training with PNNs allows for generalization to various PNN architectures, paving the way for scalable and reliable next-generation photonic computing systems.
Bacterial flagellar swarming enables dense microbial populations to migrate collectively across surfaces, often resulting in emergent, coordinated behaviors. However, probing the underlying energetics of swarming at the single cluster level remains a challenge. Here, we combine optical tweezers and multiparticle tracking within a stochastic thermodynamic framework to characterize the active motility of confined Proteus mirabilis clusters. Using the Photon Momentum Method to directly measure trapping forces, we show that swarming clusters generate persistent, dissipative flows indicative of non equilibrium stationary motility within confined solenoidal mesostructures. These flagellar rotational dynamics break detailed balance in mesoscopic force space and exceed the limits of passive friction, as evidenced by force velocity correlations and vortex like circulations. By coarse graining cluster trajectories into an active Brownian phase space, we quantify the work performed by bacterial swarms at cooperative coupling to thermal fluctuations, resulting in dissipative Ohmic like currents overcoming conservative trapping. Our findings establish a generalizable approach to quantify collective motility and energetic dissipation in active bacterial clusters, offering new insights into the physical principles governing microbial cooperativity.
We present and compare three approaches for accurately retrieving depth-resolved temperature distributions within materials from their thermal-radiation spectra, based on: (1) a nonlinear equation solver implemented in commercial software, (2) a custom-built nonlinear equation solver, and (3) a deep neural network (DNN) model. These methods are first validated using synthetic datasets comprising randomly generated temperature profiles and corresponding noisy thermal-radiation spectra for three different structures: a fused-silica substrate, an indium antimonide substrate, and a thin-film gallium nitride layer on a sapphire substrate. We then assess the performance of each approach using experimental spectra collected from a fused-silica window heated on a temperature-controlled stage. Our results demonstrate that the DNN-based method consistently outperforms conventional numerical techniques on both synthetic and experimental data, providing a robust solution for accurate depth-resolved temperature profiling.
Electromagnetically Induced Transparency (EIT) cooling is a well-established method for preparing trapped ion systems in their motional ground state. However, isolating a three-level system, as required for EIT cooling, is often challenging or impractical. In this work, we extend the EIT cooling framework to multilevel systems where the number of ground states exceeds the number of excited states, ensuring the presence of at least one dark state. We develop a formalism to accurately determine the cooling rate in the weak sideband coupling regime and provide an approximate estimate for cooling rates beyond this regime, without the need for explicit simulation of the motional degree of freedom. We clarify the connection between the cooling rate and the absorption spectrum, offering a pathway for efficient near-ground-state cooling of ions with complex electronic structures.
Holographic displays have significant potential in virtual reality and augmented reality owing to their ability to provide all the depth cues. Deep learning-based methods play an important role in computer-generated holograms (CGH). During the diffraction process, each pixel exerts an influence on the reconstructed image. However, previous works face challenges in capturing sufficient information to accurately model this process, primarily due to the inadequacy of their effective receptive field (ERF). Here, we designed complex-valued deformable convolution for integration into network, enabling dynamic adjustment of the convolution kernel's shape to increase flexibility of ERF for better feature extraction. This approach allows us to utilize a single model while achieving state-of-the-art performance in both simulated and optical experiment reconstructions, surpassing existing open-source models. Specifically, our method has a peak signal-to-noise ratio that is 2.04 dB, 5.31 dB, and 9.71 dB higher than that of CCNN-CGH, HoloNet, and Holo-encoder, respectively, when the resolution is 1920$\times$1072. The number of parameters of our model is only about one-eighth of that of CCNN-CGH.
Defect-based quantum emitters in hexagonal boron nitride (hBN) are promising building blocks for scalable quantum photonics due to their stable single-photon emission at room temperature. However, enhancing their emission intensity and controlling the decay dynamics remain significant challenges. This study demonstrates a low-cost, scalable fabrication approach to integrate plasmonic nanocavities with defect-based quantum emitters in hBN nanoflakes. Using the thermal dewetting process, we realize two distinct configurations: stochastic Ag nanoparticles (AgNPs) on hBN flakes and hybrid plasmonic nanocavities formed by AgNPs on top of hBN flakes supported on gold/silicon dioxide (Au/SiO2) substrates. While AgNPs on bare hBN yield up to a two-fold photoluminescence (PL) enhancement with reduced emitter lifetimes, the hybrid nanocavity architecture provides a dramatic, up to 100-fold PL enhancement and improved uniformity across multiple. emitters, all without requiring deterministic positioning. Finite-difference time-domain (FDTD) simulations and time-resolved PL measurements confirm size-dependent control over decay dynamics and cavity-emitter interactions. Our versatile solution overcomes key quantum photonic device development challenges, including material integration, emission intensity optimization, and spectral multiplexity. Future work will explore potential applications in integrated photonic circuits hosting on-chip quantum systems and hBN-based label-free single-molecule detection through such quantum nanoantennas.
Although the Smoothed Particle Hydrodynamics (SPH) method has been demonstrated as a promising numerical solver for multiphase flow problems due to its Lagrangian nature, its application to complex channel flow may encounter additional issues such as the open boundary condition and numerical instability in the extreme flow state. The present work aims to establish a general SPH algorithm for accurate and stable simulation of complex multiphase flows in configurations with open boundaries. The general scheme of weakly compressible SPH is adopted with special treatments implemented to ease the numerical oscillation in the density discontinuity scenario, and the turbulent model is implemented for interpretations of extreme flow conditions in high Reynolds numbers. Then, the conventional open boundary condition is fine-tuned by two new algorithms to guarantee the numerical stability in the inflow and the outflow regions. Firstly, a density relaxation is proposed to alleviate the pressure instability in the inflow region, which improves the smoothness of the particle pre-processing procedure. Secondly, the particle shifting technique with adaptive damper is implemented to adjust the magnitude of correction in the outflow region, which helps to suppress the velocity oscillation near the outlet. Validations of the proposed algorithms are carried out through four classic numerical examples, presenting appealing agreements with the analytical solutions and thus demonstrating the versatility in various flow conditions. Then, the robustness of the method is established through the turbulent multiphase channel flow cases and the horizontal slug channel flow with large density ratios. These results shed light on the value of the proposed algorithm as a general solver for complex multiphase flow problems with open boundaries.
Tracking the multifarious ultrafast electronic and structural changes occurring in a molecule during a photochemical transformation is a challenging endeavor that benefits from recent experimental and computational progress in time-resolved techniques. Measurements of valence electronic states, which provide a global picture of the bonding structure of the molecule, and core electronic states, which provide insight into the local environment, traditionally require different approaches and are often studied separately. Here, we demonstrate that X-ray pulses from a seeded free-electron laser (FEL) enable the measurement of high-resolution, time-resolved X-ray photoelectron spectra (XPS) that capture weak satellite states resulting from shake-down processes in a valence-excited molecule. This approach effectively combines the advantages of both valence- and core-state investigations. We applied this method to investigate photoexcited CS$_2$ molecules, where the role of internal conversion (IC) and intersystem crossing (ISC) in determining the pre-dissociation dynamics is controversial. We present XPS spectra from photoexcited CS$_2$, obtained at the FERMI FEL. High-resolution measurements, compared to the corresponding spectra obtained from accurate multireference quantum chemical calculations, reveal that shake-down satellite channels are highly sensitive to both valence electronic and geometric changes. Previous studies of the pre-dissociation dynamics have led to uncertain assignments around the branching between singlet and triplet excited states. We derived a propensity rule that demonstrates the spin-selectivity of the shake-downs. This selectivity allows us to unequivocally assign contributions from the bright and dark singlet excited states with populations tracked along the pre-dissociation dynamic pathway.
Accurate thermochemical data with sub-chemical accuracy (i.e., within $\pm$1 kcal mol$^{-1}$ from sufficiently accurate experimental or theoretical reference data) is essential for the development and improvement of computational chemistry methods. Challenging thermochemical properties such as heats of formation and total atomization energies (TAEs) are of particular interest because they rigorously test the ability of computational chemistry methods to accurately describe complex chemical transformations involving multiple bond rearrangements. Yet, existing thermochemical datasets that confidently reach this level of accuracy are limited in either size or scope. Datasets with highly accurate reference values include a small number of data points, and larger datasets provide less accurate data or only cover a narrow portion of the chemical space. The existing datasets are therefore insufficient for developing data-driven methods with predictive accuracy over a large chemical space. The Microsoft Research Accurate Chemistry Collection (MSR-ACC) will address this challenge. Here, it offers the MSR-ACC/TAE25 dataset of 76,879 total atomization energies obtained at the CCSD(T)/CBS level via the W1-F12 thermochemical protocol. The dataset is constructed to exhaustively cover chemical space for all elements up to argon by enumerating and sampling chemical graphs, thus avoiding bias towards any particular subspace of the chemical space (such as drug-like, organic, or experimentally observed molecules). With this first dataset in MSR-ACC, we enable data-driven approaches for developing predictive computational chemistry methods with unprecedented accuracy and scope.
Accurately modeling wind turbine wakes is essential for optimizing wind farm performance but remains a persistent challenge. While the dynamic wake meandering (DWM) model captures unsteady wake behavior, it suffers from near-wake inaccuracies due to empirical closures. We propose a Symbolic Regression-enhanced DWM (SRDWM) framework that achieves equation-level closure by embedding symbolic expressions for volumetric forcing and boundary terms explicitly into governing equations. These physically consistent expressions are discovered from LES data using symbolic regression guided by a hierarchical, domain-informed decomposition strategy. A revised wake-added turbulence formulation is further introduced to enhance turbulence intensity predictions. Extensive validation across varying inflows shows that SRDWM accurately reproduces both mean wake characteristics and turbulent dynamics, achieving full spatiotemporal resolution with over three orders of magnitude speedup compared to LES. The results highlight symbolic regression as a bridge between data and physics, enabling interpretable and generalizable modeling.
This study aims to offer a new perspective on complex network representation of real-world systems. Currently, the most well-known transformation algorithms in the literature treat each data point in a time series as a node and transform the time series into a network. In this study, we present a new approach converting spatio-temporal series into a complex network. We focus on studying this transformation by grounding it in the context of physics , with the aim of adapting it to real-world problems, which often manifest as complex systems across various domains. We introduce the Gravitational Graph (GG) algorithm, which is grounded in the concept of gravitational force from fundamental physics. We consider air pollution concentrations, which represent a global environmental health risk, as an example of a complex environmental system, and apply the GG algorithm to particulate matter of 10 microns (PM10) recorded by 21 air quality monitoring stations located in various regions of Istanbul, Turkiye. While the GG algorithm allows for the conversion of spatio-temporal series -- rather than time series -- into networks, it also enables the analysis of the statistical properties and characteristics of the converted networks, thereby uncovering hidden relationships and dependencies that may not be apparent in the original time series.
Reliable fundamental parameters are essential for accurate X-ray fluorescence analysis of sodium. This work presents updated experimental values with reliable uncertainties for sodium K-shell fundamental parameters such as: fluorescence yield, photoionization cross sections, and Auger yields. Using a physically calibrated setup and a thin NaCl layer on a silicon nitride membrane, a holistic determination approach was applied to reduce uncertainties. The new values represent a significant improvement on the ones widely used in databases. All data are available via Zenodo to support precise sodium quantification in scientific and industrial applications.
An immersed-boundary method for the incompressible Navier--Stokes equations is presented. It employs discrete forcing for a sharp discrimination of the solid-fluid interface, and achieves second-order accuracy, demonstrated in examples with highly complex three-dimensional geometries. The method is implicit, meaning that the point in the solid which is nearest to the interface is accounted for implicitly, which benefits stability and convergence properties; the correction is also implicit in time (without requiring a matrix inversion), although the temporal integration scheme is fully explicit. The method stands out for its simplicity and efficiency: when implemented alongside second-order finite differences, only the weight of the center point of the Laplacian stencil in the momentum equation is modified, and no corrections for the continuity equation and the pressure are required. The immersed-boundary method, its performance and its accuracy are first verified on simple problems, and then put to test on a simple laminar, two-dimensional flow and on two more complex examples: the turbulent flow in a channel with a sinusoidal wall, and the flow in a human nasal cavity, whose extreme anatomical complexity mandates an accurate treatment of the boundary.