Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear elimination”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fermilab PIP II Machine Protection System Digitized Data Noise Elimination Scheme and Its FPGA Implementation

In Fermilab's PIP-II machine protection system, beam loss signals from various detectors are digitized at 125 MS/s. Noise from both high-frequency sources and low-frequency 60 Hz AC power equipment can con-taminate the data. To suppress noise across these ranges especially 60 Hz and its harmonics, which overlap with beam loss signal frequencies advanced digital processing beyond standard filtering is re-quired. Several real-time functional blocks were simu-lated and tested on an FPGA: (1) a dual time-constant discharging integrator filter, (2) a de-ripple baseline extraction and storage block, and (3) a fast-recovery discharging integrator. The nonlinear IIR integrator filter removes high-frequency noise and feeds into the baseline extractor. Upon detecting abrupt beam loss, it switches to a longer time constant to prevent baseline distortion. The de-ripple block calculates a valid base-line by averaging over multiple 60 Hz periods, storing results in a 4096-word FPGA RAM. This baseline is subtracted from raw data before integration by the fast-recovery block, which resets quickly after use. All blocks achieved expected performance.

Wu, Jinyuan [Fermilab] (ORCID:0000000344329521)↗
DOE OSTI2021

Spatially Quasi-Periodic Water Waves of Infinite Depth

Abstract We formulate the two-dimensional gravity-capillary water wave equations in a spatially quasi-periodic setting and present a numerical study of solutions of the initial value problem. We propose a Fourier pseudo-spectral discretization of the equations of motion in which one-dimensional quasi-periodic functions are represented by two-dimensional periodic functions on a torus. We adopt a conformal mapping formulation and employ a quasi-periodic version of the Hilbert transform to determine the normal velocity of the free surface. Two methods of time-stepping the initial value problem are proposed, an explicit Runge–Kutta (ERK) method and an exponential time-differencing (ETD) scheme. The ETD approach makes use of the small-scale decomposition to eliminate stiffness due to surface tension. We perform a convergence study to compare the accuracy and efficiency of the methods on a traveling wave test problem. We also present an example of a periodic wave profile containing vertical tangent lines that is set in motion with a quasi-periodic velocity potential. As time evolves, each wave peak evolves differently, and only some of them overturn. Beyond water waves, we argue that spatial quasi-periodicity is a natural setting to study the dynamics of linear and nonlinear waves, offering a third option to the usual modeling assumption that solutions either evolve on a periodic domain or decay at infinity.

Wilkening, Jon (ORCID:0000000327827596)↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Posterior comparison of model dynamics in several hybrid turbulence model forms

Hybrid turbulence models that can accurately reproduce unsteady three-dimensional flow physics across the entire range of grid scales and turbulence dynamics from Reynolds-averaged Navier–Stokes (RANS), through large-eddy simulation (LES), down to direct numerical simulations (DNS) are of increasing interest to the turbulence modeling community. However, despite decades of research and development, the basic tasks of eliminating poor-performing hybrid RANS-LES models and accelerating adoption of superior models through well-designed validation and verification have yet to occur. As a step in this direction, in this work we evaluate thirteen different hybrid RANS-LES models via systematic grid refinement of decaying homogeneous isotropic turbulence. We further derive a novel mathematical framework for assessing the energy partitioning dynamics of each Hybrid RANS-LES model, wherein model-to-model variations in energy partitioning can be interpreted as different feedback mechanisms operating on a low-dimensional nonlinear dynamical system. We found that model forms similar to the flow simulation methodology—also often termed very-large eddy simulation—are dynamically inconsistent with DNS at all resolutions. Additionally, we found a strong dynamical similarity in the feedback mechanisms of all models related to detached eddy simulation and partially averaged Navier–Stokes that is inherent to their general model forms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Image charge interaction correction in charged-defect calculations

Charged-defect calculation using a periodic supercell is a significant class of problems in solid state physics. However, the finite supercell size induces an undesirable long-range image charge Coulomb interaction. Although a variety of methods have been proposed to eliminate such image Coulomb interaction, most of the previous schemes are based on a rough approximation of the defect charge screening. In this work, we present a rigorous derivation of the image charge interaction with a defect screening model where the use of a bulk macroscopic dielectric constant can be avoided. Here, we have verified this approach in comparison with a widely used approach for 12 different defects. Our correction scheme offers a much faster convergence concerning the supercell size for cases with considerable image charge interactions. In those cases, we also found that the nonlinear dielectric screening might play an important role. Our proposed defect screening model will also shed new light on understanding the defect screening properties and can be applied to other defect systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Semiclassical theory and the Koopman-van Hove equation *

Abstract The phase space Koopman-van Hove (KvH) equation can be derived from the asymptotic semiclassical analysis of partial differential equations. Semiclassical theory yields the Hamilton–Jacobi equation for the complex phase factor and the transport equation for the amplitude. These two equations can be combined to form a nonlinear semiclassical version of the KvH equation in configuration space. There is a natural injection of configuration space solutions into phase space and a natural projection of phase space solutions onto configuration space. Hence, every solution of the configuration space KvH equation satisfies both the semiclassical phase space KvH equation and the Hamilton–Jacobi constraint. For configuration space solutions, this constraint resolves the paradox that there are two different conserved densities in phase space. For integrable systems, the KvH spectrum is the Cartesian product of a classical and a semiclassical spectrum. If the classical spectrum is eliminated, then, with the correct choice of Jeffreys–Wentzel–Kramers–Brillouin (JWKB) matching conditions, the semiclassical spectrum satisfies the Einstein–Brillouin–Keller quantization conditions which include the correction due to the Maslov index. However, semiclassical analysis uses different choices for boundary conditions, continuity requirements, and the domain of definition. For example, use of the complex JWKB method allows for the treatment of tunneling through the complexification of phase space. Finally, although KvH wavefunctions include the possibility of interference effects, interference is not observable when all observables are approximated as local operators on phase space. Observing interference effects requires consideration of nonlocal operations, e.g. through higher orders in the asymptotic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

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

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Parasitic conduction loss of lithium niobate on insulator platform

Thin-film lithium niobate on insulator (LNOI) has recently emerged as an exciting platform for nonlinear integrated photonic devices, drawing parallels with the impact of silicon on insulator (SOI) in the electronic device realm. While the influence of interfacial states of SOI on electronic device performance, particularly their contribution to parasitic radio frequency (RF) losses, is well-documented, this aspect remains relatively unexplored in LNOI, despite their shared silicon/oxide interface structure. Here, we investigate the RF loss arising from mobile charges at the silicon/oxide interface of the LNOI platform. We show that this parasitic conduction can be exacerbated, especially during plasma etching, and leads to substantial RF losses. Through strategic incorporation of a trap-rich poly-silicon layer at this interface, we effectively eliminate this parasitic conduction, resulting in the suppression of RF losses. This work establishes crucial guidelines for harnessing LNOI in applications such as high-frequency electro-mechanical devices, high-speed electro-optic modulators, and high-efficiency microwave-to-optics converters.

Physics↗
DOE OSTI2022

A finite micro-rotation material point method for micropolar solid and fluid dynamics with three-dimensional evolving contacts and free surfaces

This paper introduces an explicit material point method designed specifically for simulating the micropolar continuum dynamics in the finite deformation and finite microrotation regime. The material point method enables us to simulate large deformation problems while circumventing the potential mesh distortion without remeshing. To eliminate rotational motion damping and loss of angular momentum during the projection, we introduce the mapping for microinertia and angular momentum between particles and grids through the affine particle-in-cell approach. The microrotation and the curvature at each particle are updated through zero-order forward integration of the microgyration and its spatial gradient. We show that the microinertia and the angular momentum are conserved during the projections between particles and grids in our formulation. We verify the formulation and implementation by comparing with the analytical dispersion relation of micropolar waves under the small strain and small microrotation, as well as the analytical soliton solution for solids undergoing large deformation and large microrotation. Additionally we also demonstrate the capacity of the proposed computational framework to handle a wide spectrum of simulations that exhibit size effects in the geometrical nonlinear regime through three representative numerical examples, i.e., a cantilever beam torsion problem, a fragment-impact penetration problem, and a micropolar fluid discharging problem.

42 ENGINEERING↗
DOE OSTI2021

Design and Analysis of High Heat Flux Plasma-Facing Components for NSTX Upgrade

The National Spherical Tokamak Experiment (NSTX) upgrade (NSTX-U) requirements lead to enhanced heat loads on plasma-facing components (PFCs) especially in the divertor regions, where normal heat flux density can reach up to 8 MW/m 2 . For these high heat flux (HHF) regions, the PFCs were redesigned, to use castellations which relieved the thermal stresses associated with high incident heat flux. Another design feature of HHF PFCs is the absence of front face mounting holes which create localized areas of high temperature and thermal stress concentrations. Optimized shaping of the front face of the HHF PFCs eliminates regions where the front face is perpendicular to the heat flux direction during normal helicity operation and spreads the heat load. A unique mechanism was designed to connect PFC tiles to the NSTX-U center stack casing using locking features accessible from the low heat flux regions. Isotropic graphite was selected as the HHF PFC tile material. Initial tile design parameters were accessed using analytical expressions for pulsed heat flux loading. A working prototype of the locking mechanism was created during the initial stages of the design to prove the concept performance. This article presents an overview of the divertor HHF PFC design and the results of the thermal and structural analyses performed using ANSYS software. The results of the analyses cover normal operating conditions and disruptions which impose electromagnetic (EM) loads from eddy and halo currents. The 3-D, transient, nonlinear analyses took into account the temperature-dependent properties of the materials, friction interfaces between the parts, and variable electric properties of the parts and interfaces. The results confirm that the HHF PFC tiles remain within the allowable limits for the loads defined by the NSTX-U Recovery Project. In addition, tolerance stack up analyses were performed to ensure tile performance in the worst possible assembly configuration. Finally, the design process was completed successfully, and the NSTX HHF PFC tiles are currently in the production phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Titanium Niobium Oxide: From Discovery to Application in Fast-Charging Lithium-Ion Batteries

Lithium-ion batteries are essential for portable technology and are now poised to disrupt a century of combustion-based transportation. The electrification revolution could eliminate our reliance on fossil fuels and enable a clean energy future; advanced batteries would facilitate this transition. However, owing to the demanding performance, cost, and safety requirements, it is challenging to translate new materials from laboratory prototypes to industrial-scale products. This Perspective describes that journey for a new lithium-ion battery anode material, TiNb 2 O 7 (TNO). TNO is intended as an alternative to graphite or Li 4 Ti 5 O 12 with better rate and safety characteristics than the former and higher energy density than the latter. The high capacity of TNO stems from the multielectron redox of Nb 5+ to Nb 3+ , its operating voltage window well above the Li + /Li reduction potential prevents lithium dendrite formation, and its open crystal structure leads to high-power performance. Nevertheless, the creation of a practical TNO anode was nonlinear and nontrivial. Its history is built on 30 years of fundamental science that preceded its application as a battery anode, and its battery development included a nearly 30-year gap. The insights and lessons contained in this Perspective, many of them acquired firsthand, serve two purposes: (i) to unite the disparate studies of TiNb 2 O 7 into a coherent modern understanding relevant to its application as a battery material and (ii) to highlight briefly some of the challenges faced when scaling up a new material that affect TiNb 2 O 7 as well as new electrode candidates more generally.

25 ENERGY STORAGE↗
DOE OSTI2024

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗
DOE OSTI2021

Investigating Thaw and Plant Productivity Constraints on Old Soil Carbon Respiration From Permafrost

Isotopic radiocarbon (Δ 14 C) signatures of ecosystem respiration (Reco) can identify old soil carbon (C) loss and serve as an early indicator of permafrost destabilization in a warming climate. Warming also stimulates plant productivity causing plant respiration to dominate Reco Δ 14 C signatures and potentially obscuring old soil C loss. Here, we investigate how a wide spatio-temporal gradient of permafrost thaw and plant productivity affects Reco Δ 14 C patterns and isotopic partitioning. Spatial gradients came from a warming experiment with doubling thaw depth and variable biomass, and a vegetation removal manipulation to eliminate plant contributions. We sampled in August and September to capture transitions from high to low plant productivity, decreased surface soil temperature, and relatively small seasonal thaw extensions. We found that surface processes dominate spatial variation in old soil C loss and a process-based partitioning approach was crucial for constraining old soil C loss. Resampling the same plots in different times of the year revealed that old soil C losses tripled with cooling surface temperature, and the largest old soil C losses were detected when the organic-to-mineral soil horizons thawed (~50–60 cm). We suggest that the measured increase in old soil respiration over the season and when the organic-to-mineral horizon thawed, may be explained by mobilization of nitrogen that stimulates microbial decomposition at depth. Our results suggest that soil C in the organic to mineral horizon may be an important source of soil C loss as the entire Arctic region warms and could lead to nonlinearities in projected permafrost climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Coherent anti-Stokes Raman scattering with squeezed light: CARS for quantum-enhanced spectroscopy and imaging

We theoretically investigate quantum-enhanced coherent anti-Stokes Raman scattering (CARS) using squeezed light to amplify vibrational transition rates at low photon flux. Quantum sensing approaches are needed for nondestructive nanometrology such as in bioimaging where reduced photodamage is desired while retaining resolution and sensitivity. We analyze both single-mode squeezing applied to the pump field and two-mode squeezing between the pump and Stokes fields. We also show that the ordering of displacement and squeezing operations—whether displacement precedes squeezing or squeezing precedes displacement—has an impact on the resulting CARS transition amplitudes due to a difference in the photon number and the quantum-enhancement coefficients, with the latter offering a stronger enhancement in the case of two modes squeezing of the pump and Stokes under experimentally accessible conditions. Furthermore, our calculations capture these quantum enhancements through the intrinsic photon-number correlations of squeezed light, eliminating the need for interferometric detection or higher pump powers that are otherwise required to reach comparable sensitivities in classical CARS. Finally, we outline a quantum plasmonic extension of our model in which local field enhancements caused by surface plasmon excitation in metallic nanoparticles can be incorporated via mode-selective field amplification factors, offering a pathway toward combining squeezed-light quantum optics with surface-enhanced nanoscale spectroscopy and imaging.

Atomic & molecular structure↗
DOE OSTI2024

Electrically interfaced Brillouin-active waveguide for microwave photonic measurements

New strategies for converting signals between optical and microwave domains could play a pivotal role in advancing both classical and quantum technologies. Traditional approaches to optical-to-microwave transduction typically perturb or destroy the information encoded on intensity of the light field, eliminating the possibility for further processing or distribution of these signals. In this paper, we introduce an optical-to-microwave conversion method that allows for both detection and spectral analysis of microwave photonic signals without degradation of their information content. This functionality is demonstrated using an optomechanical waveguide integrated with a piezoelectric transducer. Efficient electromechanical and optomechanical coupling within this system permits bidirectional optical-to-microwave conversion with a quantum efficiency of up to -54.16 dB. Leveraging the preservation of the optical field envelope in intramodal Brillouin scattering, we demonstrate a multi-channel microwave photonic filter by transmitting an optical signal through a series of electro-optomechanical waveguide segments, each with distinct resonance frequencies. Such electro-optomechanical systems could offer flexible strategies for remote sensing, channelization, and spectrum analysis in microwave photonics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗