Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filtering”

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,081 records · Page 60

Data-based filtered dissipation rate modelling for multi-modal turbulent combustion: evaluating a priori model generalizability

Manifold-based models offer a computationally efficient alternative to directly transporting the thermochemical state in computational simulations of turbulent reacting flows, projecting the high-dimensional thermochemical state-space onto a low-dimensional manifold. Recent efforts have yielded a manifold-based model applicable to multi-modal combustion, enabling reconstruction of the thermochemical state from solutions to two-dimensional manifold equations in mixture fraction and generalized progress variable that are parameterised by three scalar dissipation rates. In coarse-grained simulations such as Large Eddy Simulation (LES), closure of the multi-modal manifold equations and subfilter variances/covariance requires closure of three filtered scalar dissipation rates. Here, the present work adopts a data-based approach, providing closure for the three filtered scalar dissipation rates via deep neural networks (DNNs). High-fidelity datasets corresponding to an autoigniting n-dodecane jet flame and a bluff body swirl-stabilized confined lifted spray flame of two aviation fuels (Jet-A and C1) with different ignition propensities are leveraged to generate training data that spans a diverse range of thermodynamic conditions and combustion modes, including low- and high-temperature ignition regimes in addition to premixed and nonpremixed behaviour. A final DNN model is trained to enforce inherent physical constraints by learning nonlinear functional transformations of the three filtered scalar dissipation rates. The generalizability of this constrained DNN model is demonstrated a priori via conditional statistics evaluated on the lifted spray flame with C1–a configuration that had not been included in the training data. Excellent DNN agreement with conditional DNS statistics is observed, and integrated gradients are computed to identify the most sensitive input variables. The similarity of the marginal PDFs of the most informative input variables and outputs across configurations are quantified via the Wasserstein metric, demonstrating that data-based models may successfully generalize to unseen parametric conditions so long as the most informative input variables share similar distributions across training and testing datasets.

Data-based modelling↗

Rapid constriction of the selectivity filter underlies C-type inactivation in the KcsA potassium channel

C-type inactivation is a time-dependent process observed in many K+ channels whereby prolonged activation by an external stimulus leads to a reduction in ionic conduction. While C-type inactivation is thought to be a result of a constriction of the selectivity filter, the local dynamics of the process remain elusive. Here, we use molecular dynamics (MD) simulations of the KcsA channel to elucidate the nature of kinetically delayed activation/inactivation gating coupling. Microsecond-scale MD simulations based on the truncated form of the KcsA channel (C-terminal domain deleted) provide a first glimpse of the onset of C-type inactivation. We observe over multiple trajectories that the selectivity filter consistently undergoes a spontaneous and rapid (within 1–2 µs) transition to a constricted conformation when the intracellular activation gate is fully open, but remains in the conductive conformation when the activation gate is closed or partially open. Multidimensional umbrella sampling potential of mean force calculations and nonequilibrium voltage-driven simulations further confirm these observations. Electrophysiological measurements show that the truncated form of the KcsA channel inactivates faster and greater than full-length KcsA, which is consistent with truncated KcsA opening to a greater degree because of the absence of the C-terminal domain restraint. Together, these results imply that the observed kinetics underlying activation/inactivation gating reflect a rapid conductive-to-constricted transition of the selectivity filter that is allosterically controlled by the slow opening of the intracellular gate.

Li, Jing↗

TOMOGRAPHY OF PLASMA FLOWS IN THE UPPER SOLAR CONVECTION ZONE USING TIME-DISTANCE INVERSION COMBINING RIDGE AND PHASE-SPEED FILTERING

The consistency of time-distance inversions for horizontal components of the plasma flow on supergranular scales in the upper solar convection zone is checked by comparing the results derived using two k-{omega} filtering procedures-ridge filtering and phase-speed filtering-commonly used in time-distance helioseismology. I show that both approaches result in similar flow estimates when finite-frequency sensitivity kernels are used. I further demonstrate that the performance of the inversion improves (in terms of a simultaneously better averaging kernel and a lower noise level) when the two approaches are combined together in one inversion. Using the combined inversion, I invert for horizontal flows in the upper 10 Mm of the solar convection zone. The flows connected with supergranulation seem to be coherent only for the top {approx}5 Mm; deeper down there is a hint of change of the convection scales toward structures larger than supergranules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Simons Observatory: validation of reconstructed power spectra from simulated filtered maps for the small aperture telescope survey

We present a transfer function-based method to estimate angular power spectra from filtered maps for cosmic microwave background (CMB) surveys. This is especially relevant for experiments targeting the faint primordial gravitational wave signatures in CMB polarisation at large scales, such as the Simons Observatory (SO) small aperture telescopes. While timestreams can be filtered to mitigate the contamination from low-frequency noise, usual methods that calculate the mode coupling at individual multipoles can be challenging for experiments covering large sky areas or reaching few-arcminute resolution. The method we present here, although approximate, is more practical and faster for larger data volumes. We validate it through the use of simulated observations approximating the first year of SO data, going from half-wave plate-modulated timestreams to maps, and using simulations to estimate the mixing of polarisation modes induced by an example of time-domain filtering. We show its performance through an example null test and with an end-to-end pipeline that performs inference on cosmological parameters, including the tensor-to-scalar ratio r. The performance demonstration uses simulated observations at multiple frequency bands. We find that the method can recover unbiased parameters for our simulated noise levels.

CMBR experiments↗

Parallelized Kalman-Filter-Based Reconstruction of Particle Tracks on Many-Core Architectures

Faced with physical and energy density limitations on clock speed, contemporary microprocessor designers have increasingly turned to on-chip parallelism for performance gains. Algorithms should accordingly be designed with ample amounts of fine-grained parallelism if they are to realize the full performance of the hardware. This requirement can be challenging for algorithms that are naturally expressed as a sequence of small-matrix operations, such as the Kalman filter methods widely in use in high-energy physics experiments. In the High-Luminosity Large Hadron Collider (HL-LHC), for example, one of the dominant computational problems is expected to be finding and fitting charged-particle tracks during event reconstruction; today, the most common track-finding methods are those based on the Kalman filter. Experience at the LHC, both in the trigger and offline, has shown that these methods are robust and provide high physics performance. Previously we reported the significant parallel speedups that resulted from our efforts to adapt Kalman-filter-based tracking to many-core architectures such as Intel Xeon Phi. Here we report on how effectively those techniques can be applied to more realistic detector configurations and event complexity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A matched-filter approach to radio variability and transients: searching for orphan afterglows in the VAST Pilot Survey

Radio transient searches using traditional variability metrics struggle to recover sources whose evolution time-scale is significantly longer than the survey cadence. Motivated by the recent observations of slowly evolving radio afterglows at gigahertz frequency, we present the results of a search for radio variables and transients using an alternative matched-filter approach. We designed our matched-filter to recover sources with radio light curves that have a high-significance fit to power-law and smoothly broken power-law functions; light curves following these functions are characteristic of synchrotron transients, including ‘orphan’ gamma-ray burst afterglows, which were the primary targets of our search. Applying this matched-filter approach to data from Variables and Slow Transients Pilot Survey conducted using the Australian SKA Pathfinder, we produced five candidates in our search. Subsequent Australia Telescope Compact Array observations and analysis revealed that: one is likely a synchrotron transient; one is likely a flaring active galactic nucleus, exhibiting a flat-to-steep spectral transition over 4 months; one is associated with a starburst galaxy, with the radio emission originating from either star formation or an underlying slowly evolving transient; and the remaining two are likely extrinsic variables caused by interstellar scintillation. The synchrotron transient, VAST J175036.1–181454, has a multifrequency light curve, peak spectral luminosity, and volumetric rate that is consistent with both an off-axis afterglow and an off-axis tidal disruption event; interpreted as an off-axis afterglow would imply an average inverse beaming factor $\langle f^{-1}_{\text{b}} \rangle = 860^{+1980}_{-710}$, or equivalently, an average jet opening angle of $\langle \theta _{\textrm {j}} \rangle = 3^{+4}_{-1}\,$ deg.

79 ASTRONOMY AND ASTROPHYSICS↗

Automatic Kalman-filter-based wavelet shrinkage denoising of 1D stellar spectra

We propose a non-parametric method to denoise 1D stellar spectra based on wavelet shrinkage followed by adaptive Kalman thresholding. Wavelet shrinkage denoising involves applying the discrete wavelet transform (DWT) to the input signal, ‘shrinking’ certain frequency components in the transform domain, and then applying inverse DWT to the reduced components. The performance of this procedure is influenced by the choice of base wavelet, the number of decomposition levels, and the thresholding function. Typically, these parameters are chosen by ‘trial and error’, which can be strongly dependent on the properties of the data being denoised. We here introduce an adaptive Kalman-filter-based thresholding method that eliminates the need for choosing the number of decomposition levels. We use the ‘Haar’ wavelet basis, which we found to provide excellent filtering for 1D stellar spectra, at a low computational cost. We introduce various levels of Poisson noise into synthetic PHOENIX spectra, and test the performance of several common denoising methods against our own. It proves superior in terms of noise suppression and peak shape preservation. Finally, we expect it may also be of use in automatically and accurately filtering low signal-to-noise galaxy and quasar spectra obtained from surveys such as SDSS, Gaia, LSST, PESSTO, VANDELS, LEGA-C, and DESI.

79 ASTRONOMY AND ASTROPHYSICS↗

Valley filtering in strain-induced α-$\mathcal{T}_3$ quantum dots

We test the valley filtering capabilities of a quantum dot inscribed by locally straining an α-$\mathcal{T}_3$ lattice. Specifically, we consider an out-of-plane Gaussian bump in the center of a four-terminal configuration and calculate the generated pseudo-magnetic field having opposite direction for electrons originating from different valleys, the resulting valley-polarized currents, and the conductance between the injector and collector situated opposite one another. Depending on the quantum dot’s width and width-to-height ratio, we detect different transport regimes with and without valley filtering for both the α-$\mathcal{T}_3$ and Dice lattice structures. In addition, we analyze the essence of the conductance resonances with high valley polarization in terms of related (pseudo) Landau levels, the spatial distribution of the local density of states and the local current densities. The observed local charge and current density patterns reflect the local inversion symmetry-breaking by the strain, besides the global inversion symmetry breaking due to the scaling parameter α. By this way we can also filter out different sublattices

36 MATERIALS SCIENCE↗

Soil bacterial communities are shaped by temporal and environmental filtering: evidence from a long‐term chronosequence

Summary Soil microbial communities are abundant, hyper‐diverse and mediate global biogeochemical cycles, but we do not yet understand the processes mediating their assembly. Current hypothetical frameworks suggest temporal (e.g. dispersal limitation) and environmental (e.g. soil pH ) filters shape microbial community composition; however, there is limited empirical evidence supporting this framework in the hyper‐diverse soil environment, particularly at large spatial (i.e. regional to continental) and temporal (i.e. 100 to 1000 years) scales. Here, we present evidence from a long‐term chronosequence (4000 years) that temporal and environmental filters do indeed shape soil bacterial community composition. Furthermore, nearly 20 years of environmental monitoring allowed us to control for potentially confounding environmental variation. Soil bacterial communities were phylogenetically distinct across the chronosequence. We determined that temporal and environmental factors accounted for significant portions of bacterial phylogenetic structure using distance‐based linear models. Environmental factors together accounted for the majority of phylogenetic structure, namely, soil temperature (19%), pH (17%) and litter carbon:nitrogen (C:N; 17%). However, of all individual factors, time since deglaciation accounted for the greatest proportion of bacterial phylogenetic structure (20%). Taken together, our results provide empirical evidence that temporal and environmental filters act together to structure soil bacterial communities across large spatial and long‐term temporal scales.

Freedman, Zachary↗

Alumina-based filters made via binder jet 3D printing of alumina powder, colloidal silica infiltration, and sintering

Alumina-based, porous filter media was made via a binder jet 3D printing process consisting of an alumina powder printing step with subsequent heating, colloidal silica infiltration, drying, and sintering to consolidate particles yet retain a net open porous microstructure. The composites made were alumina-silica or alumina-mullite, where the silica sintering aid was used to densify and join the alumina particles. The resulting composite structures had open porosities in the 25–31 vol% range as measured by Archimedes density. Pressure drops were measured across the filter media at constant flow rates to compare disc shapes and complex, 3D printed filters based on the N95 design requirements. Complex, 3D-printed alumina composites were produced with acceptable pressure drops for N95 implementation.

36 MATERIALS SCIENCE↗

Solid-state lighting life prediction using extended Kalman filter

Solid-state lighting (SSL) luminaires containing light emitting diodes (LEDs) have the potential of seeing excessive temperatures when being transported across country or being stored in non-climate controlled warehouses. They are also being used in outdoor applications in desert environments that see little or no humidity but will experience extremely high temperatures during the day. This makes it important to increase our understanding of what effects high temperature exposure for a prolonged period of time will have on the usability and survivability of these devices. The U.S. Department of Energy has made a long term commitment to advance the efficiency, understanding and development of solid-state lighting (SSL) and is making a strong push for the acceptance and use of SSL products to reduce overall energy consumption attributable to lighting. Traditional light sources “burn out” at end-of-life. For an incandescent bulb, the lamp life is defined by B50 life. However, the LEDs have no filament to “burn”. The LEDs continually degrade and the light output decreases eventually below useful levels causing failure. Presently, the TM-21 test standard is used to predict the L70 life of SSL Luminaires from LM-80 test data. The TM-21 model uses an Arrhenius Equation with an Activation Energy, Pre-decay factor and Decay Rates. Several failure mechanisms may be active in a luminaire at a single time causing lumen depreciation. The underlying TM-21 Arrhenius Model may not capture the failure physics in presence of multiple failure mechanisms. Correlation of lumen maintenance with underlying physics of degradation at system-level is needed. In this paper, a Kalman Filter and Extended Kalman Filters have been used to develop a 70% Lumen Maintenance Life Prediction Model for a LEDs used in SSL luminaires. This model can be used to calculate acceleration factors, evaluate failure-probability and identify ALT methodologies for reducing test time. Ten-thousand hour LM-80 test data for various LEDs have been used for model development. System state has been described in state space form using the measurement of the feature vector, velocity of feature vector change and the acceleration of the feature vector change. System state at each future time has been computed based on the state space at preceding time step, system dynamics matrix, control vector, control matrix, measurement matrix, measured vector, process noise and measurement noise. The future state of the lumen depreciation has been estimated based on a second order Kalman Filter model and a Bayesian Framework. The measured state variable has been related to the underlying damage using physics-based models. Life prediction of L70 life for the LEDs used in SSL luminaires from KF and EKF based models have been compared with the TM-21 model predictions and experimental data.

54 ENVIRONMENTAL SCIENCES↗

Prediction of L70 lumen maintenance and chromaticity for LEDs using extended Kalman filter models

Solid-state lighting (SSL) luminaires containing light emitting diodes (LEDs) have the potential of seeing excessive temperatures when being transported across country or being stored in non-climate controlled warehouses. They are also being used in outdoor applications in desert environments that see little or no humidity but will experience extremely high temperatures during the day. This makes it important to increase our understanding of what effects high temperature exposure for a prolonged period of time will have on the usability and survivability of these devices. Traditional light sources “burn out” at end-of-life. For an incandescent bulb, the lamp life is defined by B50 life. However, the LEDs have no filament to “burn”. The LEDs continually degrade and the light output decreases eventually below useful levels causing failure. Presently, the TM-21 test standard is used to predict the L70 life of LEDs from LM-80 test data. Several failure mechanisms may be active in a LED at a single time causing lumen depreciation. The underlying TM-21 Model may not capture the failure physics in presence of multiple failure mechanisms. Correlation of lumen maintenance with underlying physics of degradation at system-level is needed. In this paper, Kalman Filter (KF) and Extended Kalman Filters (EKF) have been used to develop a 70-percent Lumen Maintenance Life Prediction Model for LEDs used in SSL luminaires. Ten-thousand hour LM-80 test data for various LEDs have been used for model development. System state at each future time has been computed based on the state space at preceding time step, system dynamics matrix, control vector, control matrix, measurement matrix, measured vector, process noise and measurement noise. The future state of the lumen depreciation has been estimated based on a second order Kalman Filter model and a Bayesian Framework. The measured state variable has been related to the underlying damage using physics-based models. Life prediction of L70 life for the LEDs used in SSL luminaires from KF and EKF based models have been compared with the TM-21 model predictions and experimental data.

Lall, Pradeep↗

Flattening filter removal for improved image quality of megavoltage fluoroscopy

Purpose: Removal of the linear accelerator (linac) flattening filter enables a high rate of dose deposition with reduced treatment time. When used for megavoltage imaging, an unflat beam has reduced primary beam scatter resulting in sharper images. In fluoroscopic imaging mode, the unflat beam has higher photon count per image frame yielding higher contrast-to-noise ratio. The authors’ goal was to quantify the effects of an unflat beam on the image quality of megavoltage portal and fluoroscopic images.Methods: 6 MV projection images were acquired in fluoroscopic and portal modes using an electronic flat-panel imager. The effects of the flattening filter on the relative modulation transfer function (MTF) and contrast-to-noise ratio were quantified using the QC3 phantom. The impact of FF removal on the contrast-to-noise ratio of gold fiducial markers also was studied under various scatter conditions.Results: The unflat beam had improved contrast resolution, up to 40% increase in MTF contrast at the highest frequency measured (0.75 line pairs/mm). The contrast-to-noise ratio was increased as expected from the increased photon flux. The visualization of fiducial markers was markedly better using the unflat beam under all scatter conditions, enabling visualization of thin gold fiducial markers, the thinnest of which was not visible using the unflat beam.Conclusions: The removal of the flattening filter from a clinical linac leads to quantifiable improvements in the image quality of megavoltage projection images. These gains enable observers to more easily visualize thin fiducial markers and track their motion on fluoroscopic images.

61 RADIATION PROTECTION AND DOSIMETRY↗

Filtered-backprojection reconstruction for a cone-beam computed tomography scanner with independent source and detector rotations

Purpose: A new cone-beam CT scanner for image-guided radiotherapy (IGRT) can independently rotate the source and the detector along circular trajectories. Existing reconstruction algorithms are not suitable for this scanning geometry. The authors propose and evaluate a three-dimensional (3D) filtered-backprojection reconstruction for this situation. Methods: The source and the detector trajectories are tuned to image a field-of-view (FOV) that is offset with respect to the center-of-rotation. The new reconstruction formula is derived from the Feldkamp algorithm and results in a similar three-step algorithm: projection weighting, ramp filtering, and weighted backprojection. Simulations of a Shepp Logan digital phantom were used to evaluate the new algorithm with a 10 cm-offset FOV. A real cone-beam CT image with an 8.5 cm-offset FOV was also obtained from projections of an anthropomorphic head phantom. Results: The quality of the cone-beam CT images reconstructed using the new algorithm was similar to those using the Feldkamp algorithm which is used in conventional cone-beam CT. The real image of the head phantom exhibited comparable image quality to that of existing systems. Conclusions: The authors have proposed a 3D filtered-backprojection reconstruction for scanners with independent source and detector rotations that is practical and effective. This algorithm forms the basis for exploiting the scanner’s unique capabilities in IGRT protocols.

60 APPLIED LIFE SCIENCES↗

SU-F-T-511: Feasibility Study of Using Flattening-Filter-Free Photon Beams to Deliver Conventional Flat Beam

Purpose: Various dosimetric benefits such as increased dose rate, and reduced leakage and out of field dose has led to the growth of FFF technology in the clinic. In this study, we concentrate on investigating the feasibility of using flattening-filter-free (FFF) beams to deliver conventional flat beams (CFB), since completely getting rid of the flattening-filter module from the gantry head can not only simplify the gantry design but also decrease the workload on machine maintain and quality assurance. Methods: The sliding window based IMRT technique was utilized to generate the CFB from the FFF beam for various beam configurations on the Elekta Versa HD. The flat beam reproducibility and MU efficiency were compared for each beam configuration among FFF planning, delivery and CFB planning. Results: Compared to the CFB plan, the 3%3mm passing rates of the FFF beams from both measurement and plan are 100% and 95%(or better) for 15×15 cm{sup 2} or smaller field size and for any field size greater than 15×15 cm{sup 2}at 10 cm depth, respectively. The largest discrepancy is about 5% and typically appears at the field shoulder area. The MU increase average was 80% for FFF compared to CFB, however has a minimal impact on treatment delivery time. Conclusion: The ability to deliver conventional flat treatments is not absent when operating in FFF mode. With proper TPS manipulation and beam modulation, FFF mode can achieve reasonable flat profiles and comparable dose coverage as CFB does for various conventional treatment techniques, such as four field box, or long spine treatment techniques. The ability to deliver most clinical treatments from the same treatment unit, will allow for less quality assurance as well as maintenance, and completely eliminate the need for the flattening filter on modern linacs.

60 APPLIED LIFE SCIENCES↗

SU-G-TeP2-11: Initial Evaluation of a Novel Split-Filter Dual-Energy CT for Use in Radiation Oncology

Purpose: To perform an initial evaluation of a novel split-filter dual-energy CT (DECT) system with the goal of understanding the clinical utility and limitations of the system for radiation therapy. Methods: Several phantoms were imaged using the split-filter DECT technique on the Siemens Edge CT scanner using a range of clinically-relevant doses. The optimum-contrast reconstruction, the mixed reconstruction, and the monoenergetic reconstructions (ranging from 40 keV to 190 keV) were evaluated. Each image was analyzed for CT number accuracy, uniformity, noise, low-contrast visibility (LCV), spatial resolution and geometric distortion. For comparison purposes, all parameters were evaluated on 120 kVp single-energy CT (SECT) scans used for treatment planning, as well as, a sequential-scan DECT technique for corresponding doses. Results: For all DECT reconstructions no observable geometric distortion was found. Both the optimal-contrast and mixed images demonstrated slight improvements in LCV and noise when compared to the SECT, and slight reductions in CT number accuracy and spatial resolution. The CT numbers trended as expected for the monoenergetic reconstructions, with CT number accuracy within 50 HU for materials of density <2 g/cm3. Spatial resolution increased with energy, and for monoenergetic reconstructions >70 keV the spatial resolution exceeded that of the SECT. The noise in the monoenergetic reconstructions increased with decreasing energy. Thus, the image uniformity, signal-to-noise ratio and LCV were diminished at lower energies (70 keV). Applying iterative reconstruction techniques to the low-energy images reduced noise and improved LCV. The signal-to-noise ratio was stable for energies >100 keV. Conclusion: The initial commissioning of the novel split-filter DECT technology demonstrated favorable results for clinical implementation. The mixed reconstruction showed potential as a replacement for the treatment planning SECT. The image parameters for the monoenergetic reconstructions varied appropriately with energy. This work provides an initial understanding of the limitations and potential applications for monoenergetic imaging.

60 APPLIED LIFE SCIENCES↗

WE-DE-207B-11: Implementation of Size-Specific 3D Beam Modulation Filters On a Dedicated Breast CT Platform Using Breast Immobilization

Purpose: To implement a 3D beam modulation filter (3D-BMF) in dedicated breast CT (bCT) and develop a method for conforming the patient’s breast to a pre-defined shape, optimizing the effects of the filter. This work expands on previous work reporting the methodology for designing a 3D-BMF that can spare unnecessary dose and improve signal equalization at the detector by preferentially filtering the beam in the thinner anterior and peripheral breast regions. Methods: Effective diameter profiles were measured for 219 segmented bCT images, grouped into volume quintiles, and averaged within each group to represent the range of breast sizes found clinically. These profiles were then used to generate five size-specific computational phantoms and fabricate five size-specific UHMW phantoms. Each computational phantom was utilized for designing a size-specific 3D-BMF using previously reported methods. Glandular dose values and projection images were simulated in MCNP6 with and without the 3DBMF using the system specifications of our prototype bCT scanner “Doheny”. Lastly, thermoplastic was molded around each of the five phantom sizes and used to produce a series of breast immobilizers for use in conforming the patient’s breast during bCT acquisition. Results: After incorporating the 3D-BMF, MC simulations estimated an 80% average reduction in the detector dynamic range requirements across all phantom sizes. The glandular dose was reduced on average 57% after normalizing by the number of quanta reaching the detector under the thickest region of the breast. Conclusion: A series of bCT-derived breast phantoms were used to design size-specific 3D-BMFs and breast immobilizers that can be used on the bCT platform to conform the patient’s breast and therefore optimally exploit the benefits of the 3D-BMF. Current efforts are focused on fabricating several prototype 3D-BMFs and performing phantom scans on Doheny for MC simulation validation and image quality analysis. Research reported in this paper was supported in part by the National Cancer Institute of the National Institutes of Health under award R01CA181081. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institue of Health.

60 APPLIED LIFE SCIENCES↗

Accuracy of Some Approximate Gaussian Filters for the Navier--Stokes Equation in the Presence of Model Error

Bayesian state estimation of a dynamical system from a stream of noisy measurements is important in many geophysical and engineering applications where high dimensionality of the state space, sparse observations, and model error pose key challenges. In this work, three computationally feasible, approximate Gaussian data assimilation/filtering algorithms are considered in various regimes of turbulent 2D Navier--Stokes dynamics in the presence of model error. The first source of error arises from the necessary use of reduced models for the forward dynamics of the filters, while a particular type of representation error arises from the finite resolution of observations which mix up information about resolved and unresolved dynamics. Two stochastically parameterized filtering algorithms, referred to as cSPEKF and GCF, are compared with 3DVAR---a prototypical time-sequential algorithm known to be accurate for filtering dissipative systems for a suitably inflated “background” covariance. We provide the first evidence that the stochastically parameterized algorithms, which do not rely on detailed knowledge of the underlying dynamics and do not require covariance inflation, can compete with or outperform an optimally tuned 3DVAR algorithm, and they can overcome competing sources of error in a range of dynamical scenarios.

97 MATHEMATICS AND COMPUTING↗