Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial optimization”

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

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

At least 1,675 records · Page 93

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part II: linesweep ALE rezoning for nonconformal block-structured AMR meshes

The simulation of flows presenting contact discontinuities, vorticity, and large variations in spatial scales can be performed in a framework coupling Arbitrary Lagrangian Eulerian (ALE) algorithms and Adaptive Mesh Refinement (AMR). This coupling requires adaptation of ALE rezoning techniques to meshes containing nonconformal nodes arising from both the AMR topology and the junction of mesh blocks. Here, in this paper, we present an ALE rezoning strategy that is compatible with such meshes, and that can also act as a disentangling algorithm. Emphasis is put on an algorithm that respects intrinsic Lagrangian mesh properties in order to preserve accuracy around discontinuities. To that end, we adapt the weighted linesweep algorithm to nonconformal block-structured AMR meshes. Then, we present control parameters introduced in the method for it to be applicable in practical situations. Notably, the method is coupled to a specific metric optimization in order to palliate some shortcomings of the linesweep method. Finally, numerical test cases are presented that feature the capabilities of the ALE-AMR algorithm for flows that present discontinuities, vorticity, and a variety of scales. Notably, we show that our ALE-AMR algorithm gives results at least similar to Euler-AMR, but provides better accuracy in cases where discontinuities are involved, thanks to a method that respects the Lagrangian features of the mesh. Additionally, it enables Euler-AMR-like computations on domains with temporally varying domain boundaries.

Adaptive mesh refinement↗

Building partnerships for development of sustainable energy systems with atmospheric measurements

Atmospheric dynamics often play a critical role in the sustainability and reliability of diverse forms of energy production. This is especially true for the growing number of renewable energy deployments that harness aspects of the environment for power production. While the University of Memphis has a strong research background in energy systems, we have little experience working with the Earth and Environmental Systems Science Division (EESSD) and their associated User Facilities. Of particular interest to us is the Atmospheric Science Research and the Atmospheric Radiation Measurement (ARM) user facility to address surface-boundary layer interactions and physical phenomena. One of the major challenges for understanding and developing energy systems and management platforms is accurate modeling/forecasting of atmospheric conditions across disparate spatial and temporal scales. These conditions are often required to understand the lowest levels of the atmospheric boundary layer, but are also important to understand higher atmospheric conditions where aerosols affect cloud development. The objective of this work was to develop partnerships with national laboratories for collaboration on environmental science and its intersection with sustainable energy systems, as well as to leverage the ARM user facility data repositories to enhance our research capabilities in energy systems and their inter-dependence on environmental systems for future engagement with EESSD. Specifically, we accomplished these objectives by (1) developing collaborations with Oakridge National Laboratory ARM Data Science and Integration Group which resulted in student internships, (2) employed ARM data to develope modeling of the atmospheric boundary layer optical turbulence, and (3) optimally-sized large-scale renewable energy systems and their associated energy storage systems with ARM repository data.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Interannual Variability of Tropical Precipitation: How Well Do Climate Models Agree With Current Satellite Estimates?

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. Recent availability of the TRMM (Tropical Rainfall Measuring Mission) has further aided in narrowing uncertainties in rainfall over die tropics and subtropics. Although climate modeling efforts have long relied on space-based precipitation estimates for validation, we now are in a position to make more quantitative assessments of model performance, particularly in tropical regions. An integration of the CCM3 using observed SSTs as a lower boundary condition is used to examine how well this model responds to ENSO forcing in terms of anomalous precipitation. An integration of the NCEP spectral model used for the Reanalysis-H effort is also examined. This integration is run with specified SSTs, but with no data assimilation. Our analysis focuses on two aspects of inter-annual variability. First are the spatial anomalies that are indicative of dislocations in Hadley and Walker circulations. Second, we consider the ability of models to replicate observed increases in oceanic precipitation that are noted in satellite observations for large ENSO events. Finally, we consider a slab ocean version of the CCM3 model with prescribed ocean beat transports that mimic upwelling anomalies, but which still allows the surface energy balance to be predicted. This less restrictive experiment is used to understand why model experiments with specified SSTs seem to have noticeably less interannual variability in precipitation than do the satellite observations.

Robertson, Franklin R.↗

Wave optics simulations for new soft x-ray beamlines at NSLS-II

In this work, we performed fully- and partially-coherent synchrotron emission and propagation simulations with the "Synchrotron Radiation Workshop" computer code to analyze the performance of two soft X-ray beamlines under development at the National Synchrotron Light Source II: Soft X-ray Nanoprobe (SXN), and Angle-Resolved Photoemission Spectroscopy (ARPES) and Resonant Inelastic X-ray Scattering (RIXS) Imaging (ARI). The SXN beamline intends to provide high flux and high spatial resolution coherent soft X-ray imaging capabilities using both zone plate and lensless coherent imaging techniques. The ARI beamline aims to perform high flux ARPES and RIXS experiments with a focal spot size at the sample approaching 100 nm using highly-demagnifying mirrors in Kirkpatrick-Baez geometry. To accurately calculate the resolution and the degree of X-ray coherence provided by the two state-of-the-art beamlines, partial coherence effects are required to be taken into account in wave optics simulations for these two beamlines. In this talk, beamline performance parameters such as spot size, degree of coherence, flux, and energy resolution at the sample are presented. The effects of mirror surface slope errors on beamline performance were studied and some suggestions for further optimization are discussed.

36 MATERIALS SCIENCE↗

Spatially and Temporally Resolved Velocimetry for Hypersonic Flows

The development of new hypersonic flight vehicles is limited by the physical understanding that may be obtained from ground test facilities. This has motivated the present development of a temporally and spatially resolved velocimetry measurement for Sandia National Laboratories (SNL) Hypersonic Wind Tunnel (HWT) using Femtosecond Laser Electronic Excitation Tagging (FLEET). First, a multi-line FLEET technique has been created for the first time and tested in a supersonic jet, allowing simultaneous measurements of velocities along multiple profiles in a flow. Secondly, two different approaches have been demonstrated for generating dotted FLEET lines. One employs a slit mask pattern focused into points to yield a dotted line, allowing for two- or three-component velocity measurements free of contamination between components. The other dotted-line approach is based upon an optical wedge array and yields a grid of points rather than a dotted line. Two successful FLEET measurement campaigns have been conducted in SNL’s HWT. The first effort established optimal diagnostic configurations in the hypersonic environment based on earlier benchtop reproductions, including validation of the use of a 267 nm beam to boost the measurement signal-to-noise ratio (SNR) with minimal risk of perturbing the flow and greater simplicity than a comparable resonant technique at 202 nm. The same FLEET system subsequently was reconstituted to demonstrate the ability to make velocimetry measurements of hypersonic turbulence in a realistic flow field. Mean velocity profiles and turbulence intensity profiles of the shear layer in the wake of a hypersonic cone model were measured at several different downstream stations, proving the viability of FLEET as a hypersonic diagnostic.

33 ADVANCED PROPULSION SYSTEMS↗

Influence of Resonances on the Noise Performance of SQUID Susceptometers

Scanning Superconducting Quantum Interference Device (SQUID) Susceptometry simultaneously images the local magnetic fields and susceptibilities above a sample with sub-micron spatial resolution. Further development of this technique requires a thorough understanding of the current, voltage, and flux ( I V Φ ) characteristics of scanning SQUID susceptometers. These sensors often have striking anomalies in their current–voltage characteristics, which we believe to be due to electromagnetic resonances. The effect of these resonances on the performance of these SQUIDs is unknown. To explore the origin and impact of the resonances, we develop a model that qualitatively reproduces the experimentally-determined I V Φ characteristics of our scanning SQUID susceptometers. We use this model to calculate the noise characteristics of SQUIDs of different designs. We find that the calculated ultimate flux noise is better in susceptometers with damping resistors that diminish the resonances than in susceptometers without damping resistors. Such calculations will enable the optimization of the signal-to-noise characteristics of scanning SQUID susceptometers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Quenching

For about a half century, the stabilization of a steady planar deflagration on a heat-sink-type flat-flame burner has been of extraordinary service for the theoretical modeling and diagnostic probing of combusting gaseous mixtures. However, most engineering devices and most unwanted fire involve the burning of initially unmixed reactants. The most vigorous burning of initially separated gaseous fuel and oxidizer is the diffusion flame. In this useful idealization (limiting case), the reactants are converted to product at a mathematically thin interface, so no interpenetration of fuel and oxidizer occurs. This limit is of practical importance because it often characterizes the condition of optimal performance (and sometimes environmentally objectionable operation) of a combustor. A steady planar diffusion flame is most closely approached in the laboratory in the counterflow apparatus. The utility of this simple-strain-rate flow for the modeling and probing of diffusion flames was noted by Pandya and Weinberg 35 years ago, though only in the last decade or so has its use become internationally common place. However, typically, as the strain rate a is reduced below about 20 cm(exp -1), and the diffusion-flame limit (reaction rate much faster than the flow rate) is approached, the burning is observed to become unstable in earth gravity. The advantageous steady planar flow is not available in the diffusion-flame limit in earth gravity. This is unfortunate because the typical spatial scale in a counterflow is (k/a)(sup 1/2), where k denotes a characteristic diffusion coefficient; thus, the length scale becomes large, and the reacting flow is particularly amenable to diagnostic probing, as the diffusion-flame limit is approached. The disruption of planar symmetry is owing the fact that, as the strain rate a decreases, the residence time (l/a) of the throughput in the counterflow burner increases. Observationally, when the residence time exceeds about 50 msec, the inevitably present convective (Rayleigh-Benard) instabilities, associated with hot-under-cold (flame-under-fresh-reactant) stratification of fluid in a gravitational field, have time to grow to finite amplitude during transit of the burner.

Fendell, Francis↗

High-throughput micro-scale bandgap mapping for perovskite-inspired materials with complex composition space

Abstract To realize the full promise of high-throughput experimental workflows, the rate of sample synthesis must be matched by that of characterization. Of growing interest are contactless optical techniques that can rapidly measure material homogeneity and properties. Here, we present a hyperspectral imaging method to measure local optical bandgap distributions within samples, utilizing spatially-resolved reflectance spectra coupled with automated data analysis. We collect approximately one million optical bandgap data across the compositional space of Cs 3 (Bi x Sb 1-x ) 2 (Br y I 1-y ) 9 perovskite-inspired materials. Our results show non-monotonic bandgap variations (i.e., bandgap bowing) along six composition gradient sequences, in addition to identifying samples with multiple bandgaps in statistics. High-throughput transient absorption spectroscopy reveals that within these compositions, the depletion of the ground state carriers to excited states occurred at discrete energy levels with independent carrier dynamics, consistent with the bandgap observation and indicative of phase separation. This work demonstrates the potential for rapid optical measurements to assess material quality and homogeneity in a high-throughput experimental setting, supporting screening and recipe optimization of optoelectronic material candidates with desired carrier dynamics and optical properties.

Science & Technology - Other Topics↗

Characterization of latent image of electron beam resist via critical-dimension resonant soft X-ray scattering

Metrology plays a crucial role in semiconductor manufacturing by providing accurate and precise measurement and characterization of critical parameters. With the development of high-resolution extreme ultraviolet lithography (EUVL) processes, critical dimensions are shrinking to sub-10 nm. Resist materials encounter the challenge of providing heightened sensitivity and a handle on exacerbating stochastic variations. A comprehensive understanding of the chemical profile of the latent image is pivotal for mitigating stochastic effects and optimizing pattern quality. However, the subtle differences in chemistry between the exposed and unexposed regions of the resists make it extremely challenging to characterize the latent images with sub-nanometer precision. Here, we develop the metrology with critical-dimension resonant soft X-ray scattering (CD-RSoXS) to probe the chemical profiles of latent images stored in resist after exposure. The combination of absorption spectroscopy and enhanced scattering contrast makes it possible to characterize the subtle structural and chemical variations in the latent image. Moreover, the results of the measurements are compared with the simulations with a finite element method–based Maxwell solver to extract a detailed profile of the latent and developed images. We demonstrate that the CD-RSoXS technique can provide valuable insights into the high spatial resolution and local chemical sensitivity simultaneously, which is crucial to understanding the resolution limits and stochastic effects in EUVL processes.

36 MATERIALS SCIENCE↗

Design of a composite filter realizable on practical spatial light modulators

Hybrid optical correlator systems use two spatial light modulators (SLM's), one at the input plane and the other at the filter plane. Currently available SLM's such as the deformable mirror device (DMD) and liquid crystal television (LCTV) SLM's exhibit arbitrarily constrained operating characteristics. The pattern recognition filters designed with the assumption that the SLM's have ideal operating characteristic may not behave as expected when implemented on the DMD or LCTV SLM's. Therefore it is necessary to incorporate the SLM constraints in the design of the filters. In this report, an iterative method is developed for the design of an unconstrained minimum average correlation energy (MACE) filter. Then using this algorithm a new approach for the design of a SLM constrained distortion invariant filter in the presence of input SLM is developed. Two different optimization algorithms are used to maximize the objective function during filter synthesis, one based on the simplex method and the other based on the Hooke and Jeeves method. Also, the simulated annealing based filter design algorithm proposed by Khan and Rajan is refined and improved. The performance of the filter is evaluated in terms of its recognition/discrimination capabilities using computer simulations and the results are compared with a simulated annealing optimization based MACE filter. The filters are designed for different LCTV SLM's operating characteristics and the correlation responses are compared. The distortion tolerance and the false class image discrimination qualities of the filter are comparable to those of the simulated annealing based filter but the new filter design takes about 1/6 of the computer time taken by the simulated annealing filter design.

Rajan, P. K.↗

Recent Developments in Transition-Edge Strip Detectors for Solar X-Rays

LMSAL and NIST are developing position-sensitive x-ray strip detectors based on Transition Edge Sensor (TES) microcalorimeters optimized for solar physics. By combining high spectral (E/ delta E approximately equals 1600) and temporal (single photon delta t approximately equals 10 micro s) resolutions with imaging capabilities, these devices will be able to study high-temperature (>l0 MK) x-ray lines as never before. Diagnostics from these lines should provide significant new insight into the physics of both microflares and the early stages of flares. Previously, the large size of traditional TESs, along with the heat loads associated with wiring large arrays, presented obstacles to using these cryogenic detectors for solar missions. Implementing strip detector technology at small scales, however, addresses both issues: here, a line of substantially smaller effective pixels requires only two TESs, decreasing both the total array size and the wiring requirements for the same spatial resolution. Early results show energy resolutions of delta E(sub fwhm) approximately equals 30 eV and spatial resolutions of approximately 10-15 micron, suggesting the strip-detector concept is viable.

Rausch, Adam J.↗

Critical Soil Moisture Thresholds of Plant Water Stress in Terrestrial Ecosystems

Plant water stress occurs at the point when soil moisture (SM) limits transpiration, defining a critical SM threshold (θ crit ). Knowledge of the spatial distribution of θ crit is crucial for future projections of climate and water resources. Here, we use global eddy covariance observations to quantify θcrit and evaporative fraction (EF) regimes. Three canonical variables describe how EF is controlled by SM: the maximum EF (EF max ), θ crit , and slope (S) between EF and SM. We find systematic differences of these three variables across biomes. Variation in θ crit , S, and EFmax is mostly explained by soil texture, vapor pressure deficit, and precipitation, respectively, as well as vegetation structure. Dryland ecosystems tend to operate at low θ crit and show adaptation to water deficits. The negative relationship between θ crit and S indicates that dryland ecosystems minimize θcrit through mechanisms of sustained SM extraction and transport by xylem. Our results further suggest an optimal adaptation of local EF-SM response that maximizes growing-season evapotranspiration and photosynthesis.

Zheng Fu↗

Cooperation in Transmission Expansion Planning: Enhancing Grid Reliability and Efficiency Under a Changing Climate

Electricity grids are challenged to maintain reliability during more intense and frequent extreme weather events due to climate change. This challenge is exacerbated by multi-sector electrification and power sector decarbonization through increased reliance on variable renewable energy, which necessitates the expansion of transmission infrastructure. However, transmission expansion planning is often complicated by intertwined planning authorities and jurisdictions, and allocation of large capital investment needs. These factors cause authorities to manage transmission investments individually (i.e., only/mostly intraregional planning), which can lead to suboptimal transmission networks. This study investigates the potential benefits of cooperative transmission expansion planning (i.e., both intraregional and interregional planning that optimizes transmission investments across the entire physical system). Using sectoral and economic optimization, and machine learning models, it analyzes the impact of different levels of cooperation among transmission planning regions within U.S. Western Interconnection in 2019 and 2059 via an iterative investment process. Furthermore, it examines the effects of future climate change on transmission cooperation by simulating historical heat waves from 2019 under conditions of 2059. The results indicate that cooperative transmission planning leads to lower wholesale electricity prices, decreased energy outages, and reduced greenhouse gas emissions. However, the advantages of collaboration diminish during widespread heat waves, despite remaining beneficial especially for regions like California Independent System Operator with substantial solar installations. The study underscores the importance of transmission cooperation in reducing costs and enhancing reliability, emphasizing the need for strategic investments in storage to address challenges posed by future extreme weather events with varying spatial scales.

Capacity Expansion Model↗

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗

Conformal duality of the nonlinear Schrödinger equation: Theory and applications to parameter estimation

The nonlinear Schrödinger equation (NLSE) in one spatial dimension has stationary solutions similar to those of the linear Schrödinger equation (LSE) as well as more exotic solutions such as solitary waves and quantum droplets. Here, we present a newly discovered conformal duality which unifies the stationary and time-dependent traveling-wave solutions of the one-dimensional cubic-quintic NLSE, the cubic NLSE and LSE. Any two systems that are classified by the same single number called the cross ratio are related by this symmetry. Notably, the conformal duality can also be adapted in Newtonian mechanics and serves as a powerful tool for investigating physical systems that otherwise cannot be directly accessed in experiments. Further, we show that the conformal symmetry is a valuable resource to substantially improve NLSE parameter estimation from noisy empirical data by introducing an optimization afterburner. The new method therefore has far reaching practical applications for nonlinear physical systems. Published by the American Physical Society 2025

Reinhardt, David B. (ORCID:0009000409812838)↗

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗

Dynamically Reconfigurable Data Readout of Pixel Detectors for Automatic Synchronization with Data Acquisition Systems

Reconfigurable detectors with dynamically selectable sensing and readout modes are highly desirable for implementing edge computing as well as enabling advanced imaging techniques such as foveation. The concept of a camera system capable of simultaneous passive imaging and dynamic ranging in different regions of the detector is presented. Such an adaptive-autonomous detector with both spatial and temporal control requires programmable window of exposure (time frames), ability to switch between readout modes such as full-frame imaging and zero-suppressed data, modification of the number of pixel data bits and independent programmability for distinct detector regions. In this work, a method is presented for seamlessly changing time frames and readout modes without data corruption while still ensuring that the data acquisition system (DAQ) does not need to stop and resynchronize at each change of setting, thus avoiding significant dead time. Data throughput is maximized by using a minimum unique data format, rather than lengthy frame headers, to differentiate between consecutive frames. A data control and transmitter (DCT) synchronizes data transfer from the pixel to the periphery, reconfigures the data to transmit it serially off-chip, while providing optimized decision support based on a DAQ definable mode. Measurements on a test structure demonstrate that the DCT can operate at 1 GHz in a 65 nm LP CMOS process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗