Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local rational modeling”

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 19 records

A dynamic localization model with stochastic backscatter

The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.

Carati, Daniele↗

Structure of Evaporating and Combusting Sprays: Measurements and Predictions

Complete measurements of the structure of nonevaporating, evaporating and combusting sprays for sufficiently well defined boundary conditions to allow evaluation of models of these processes were obtained. The development of rational design methods for aircraft combustion chambers and other devices involving spray combustion were investigated. Three methods for treating the discrete phase are being considered: a locally homogeneous flow (LHF) model, a deterministic separated flow (DSF) model, and a stochastic separated flow (SSF) model. The main properties of these models are summarized.

Shuen, J. S.↗

The origin of the galactic emission in IRAS data

Using the high resolution IRAS data and the radial distribution of gaseous material and those of the interstellar radiation field (ISRF), a model of the Galactic infrared emission was built. The first step consisted of a separation of the diffuse emission in IRAS data from that of the well-defined strong Galactic sources. A well accepted idea is that IR emission comes from dust heated either by the ISRF and/or by internal cloud heating sources. Thus the IR galactic emission was modelled from radial distributions of gas and ISRF and the following three main hypothesis: (1) the dust-to-gas ration is the same in the whole Galaxy; (2) IR emission is proportional to local dust density; and (3) IR emission is also proportional to the local ISRF. The IR emission as modelled was integrated over each line of sight and compared with observed IRAS data. The results show that the IR diffuse component comes from dust associated with H1 and heated by the general ISRF. For the dust embedded in cold H2 component, the heating source is also the general ISRF while the warm component is explained by dust embedded in molecular clouds and heated by high-mass stars born in the close vicinity of the clouds and by disc population ISRF.

Caux, Emanuel↗

Crustal deformation along the San Andreas, California

The goal is to achieve a better understanding of the regional and local deformation and crustal straining processes in western North America, particularly the effects of the San Andreas and nearby faults on the spatial and temporal crustal deformation behavior. Construction of theoretical models based on the mechanics of coupled elastic plate, viscoelastic foundation and large scale crack mechanics provide a rational basis for the interpretation of seismic and aseismic anomalies and expedite efforts in forecasting the stability of plate boundary deformation. Special focus is placed on the three dimensional time dependent surface deformation due to localized slippage in a elastic layer coupled to a visco-elastic substrate. The numerical analysis is based on a 3-D boundary element technique. Extension to visco-elastic coupling demands the derivation of 3-D time dependent Green's function. This method was applied to analyze the viscoelastic surface displacements due to a dislocated embedded patch. Surface uplift as a function of time and position are obtained. Comparisons between surface uplift for long and short dislocated patches are made.

Li, Victor C.↗

Towards a Rational Model for the Triple Velocity Correlations of Turbulence

This paper presents a rational approach to modelling the triple velocity correlations that appear in the transport equations for the Reynolds stresses. All existing models of these correlations have largely been formulated on phenomenological grounds and are defective in one important aspect: they all neglect to allow for the dependence of these correlations on the local gradients of mean velocity. The mathematical necessity for this dependence will be demonstrated in the paper. The present contribution lies in the novel use of Group Representation Theory to determine the most general tensorial form of these correlations in terms of all the second- and third-order tensor quantities that appear in the exact equations that govern their evolution. The requisite representation did not exist in the literature and therefore had to be developed specifically for this purpose by Professor G. F. Smith. The outcome of this work is a mathematical framework for the construction of algebraic, explicit, and rational models for the triple velocity correlations that are theoretically consistent and include all the correct dependencies. Previous models are reviewed, and all are shown to be an incomplete subset of this new representation, even to lowest order.

Younis, B. A.↗

A Generalized, Compactly-Supported Correlation Function for Data Assimilation Applications

Correlation functions play an essential role in modern data assimilation, where they are used to model covariances given a set of tunable parameters or applied as tapering functions to localize covariances in ensemble-based schemes. One of the most widely-used correlation functions in data assimilation is the Gaspari and Cohn (1999) piecewise-rational, compactly-supported parametric correlation function (hereafter referred to as GC99). The GC99 correlation function is useful due to its tunable cut-off parameter c and Gaussian-like shape achieved when the parameter a is set to one-half. These properties are attractive for tapering functions in data assimilation applications. However, the GC99 correlation function is homogeneous over Euclidean 3-space and isotropic when restricted to the sphere, properties that may be less than ideal for some geophysical applications. GC99 is also compactly-supported on a sphere of fixed radius, which requires tuning of the cut-off parameter c that can depend on the specific application. This work presents a generalization of the GC99 correlation function that allows the cut-off parameter c and shape parameter a to vary over space to gain more flexibility in shape while maintaining its compact support property. The function, which we call the Generalized Gaspari Cohn (GenGC) correlation function, introduces inhomogeneity in Euclidean 3-space and anisotropy when restricted to the sphere by allowing both parameters c and a to vary, as functions, over the spatial domain. The GC99 correlation function is a special case of GenGC where the functions c and a are held constant, as fixed parameters rather than functions. The GenGC correlation function also generalizes the follow-on to the work of Gaspari and Cohn (1999) presented in Gaspari et al. (2006), which allowed a to vary while keeping c fixed. We illustrate through simple one- and two-dimensional examples the variety of inhomogeneous and anisotropic correlation functions GenGC can produce by varying c and a over space, and suggest applications where they may be useful in data assimilation, such as covariance modeling or localization. In particular, we describe how the GenGC correlation function can be used to construct covariances using correlation length and variance fields derived from dynamics. For example, the correlation length field for advective dynamics is governed by a partial differential equation (PDE) in N spatial dimensions, where N is the number of space dimensions of the state. Correlation length fields can be determined from this PDE and used with GenGC to construct the corresponding correlations. We can then approximate the full covariance by rescaling by the variance, which also satisfies a PDE in N spatial dimensions for advective dynamics. Thus we can approximate the full covariance without solving the covariance PDE, which is in 2N spatial dimensions, by solving just two PDEs each in only N spatial dimensions. This approach to evolving the correlation length and variance fields, then reconstructing the correlations using GenGC, is suggested as an alternative to current methods of covariance modeling in data assimilation algorithms.

GC99↗

Distributed digital signal processors for multi-body flexible structures

Multi-body flexible structures, such as those currently under investigation in spacecraft design, are large scale (high-order) dimensional systems. Controlling and filtering such structures is a computationally complex problem. This is particularly important when many sensors and actuators are located along the structure and need to be processed in real time. This report summarizes research activity focused on solving the signal processing (that is, information processing) issues of multi-body structures. A distributed architecture is developed in which single loop processors are employed for local filtering and control. By implementing such a philosophy with an embedded controller configuration, a supervising controller may be used to process global data and make global decisions as the local devices are processing local information. A hardware testbed, a position controller system for a servo motor, is employed to illustrate the capabilities of the embedded controller structure. Several filtering and control structures which can be modeled as rational functions can be implemented on the system developed in this research effort. Thus the results of the study provide a support tool for many Control/Structure Interaction (CSI) NASA testbeds such as the Evolutionary model and the nine-bay truss structure.

Lee, Gordon K. F.↗

Spatial Growth Modeling and High Resolution Remote Sensing Data Coupled with Air Quality Modeling to Assess the Impact of Atlanta, Georgia on the Local and Regional Environment

The growth of cities, both in population and areal extent, appears as an inexorable process. Urbanization continues at a rapid rate, and it is estimated that by the year 2025, 60 percent of the world s population will live in cities. Urban expansion has profound impacts on a host of biophysical, environmental, and atmospheric processes within an urban ecosystems perspective. A reduction in air quality over cities is a major result of these impacts. Because of its complexity, the urban landscape is not adequately captured in air quality models such as the Community Multiscale Air Quality (CMAQ) model that is used to assess whether urban areas are in attainment of EPA air quality standards, primarily for ground level ozone. This inadequacy of the CMAQ model to sufficiently respond to the heterogeneous nature of the urban landscape can impact how well the model predicts ozone levels over metropolitan areas and ultimately, whether cities exceed EPA ozone air quality standards. We are exploring the utility of high-resolution remote sensing data and urban spatial growth modeling (SGM) projections as improved inputs to a meteorological/air quality modeling system focusing on the Atlanta, Georgia metropolitan area as a case study. These growth projections include business as usual and smart growth scenarios out to 2030. The growth projections illustrate the effects of employing urban heat island mitigation strategies, such as increasing tree canopy and albedo across the Atlanta metro area, which in turn, are used to model how air temperature can potentially be moderated as impacts on elevating ground-level ozone, as opposed to not utilizing heat island mitigation strategies. The National Land Cover Dataset at 30m resolution is being used as the land use/land cover input and aggregated to the 4km scale for the MM5 mesoscale meteorological model and the CMAQ modeling schemes. Use of these data has been found to better characterize low density/suburban development as compared with USGS lkm land use/land cover data that have traditionally been used in modeling. Air quality prediction for future scenarios to 2030 is being facilitated by land use projections using a spatial growth model. Land use projections were developed using the 2030 Regional Transportation Plan developed by the Atlanta Regional Commission, the regional planning agency for the area. This allows the Georgia Environmental Protection Division to evaluate how these transportation plans will affect future air quality. The coupled SGM and air quality modeling approach provides insight on what the impacts of Atlanta s growth will be on the local and regional environment and exists as a mechanism that can be used by policy makers to make rational decisions on urban growth and sustainability for the metropolitan area in the future.

Quattrochi, Dale A.↗

Continuum models for beam-like and plate-like lattice structures

A simple, rational approach is presented for developing continuum models for large repetitive beam-like and plate-like lattices with arbitrary configurations subjected to static, thermal and dynamic loadings. The continuum models for these structures are shear flexible beams and plates. They account for local effects in the repeating element of the actual structure and are characterized by their thermoelastic strain and kinetic energies from which the equations of motion and constitutive relations can be derived. The procedure for developing the expressions for thermoelastic strain and kinetic energies of the continuum involves introducing basic assumptions regarding the variation of the temperature, displacement and strain components in one or two directions (for plate-like and beam-like lattices) and obtaining effective thermoelastic and dynamic coefficients of the continuum in terms of material properties and geometry of the original lattice structure. The high accuracy of the continuum models developed is demonstrated by means of numerical examples.

Noor, A. K.↗

Strain accumulation and surface deformation along the San Andreas, California

The goal of this project remains to be the achievement of a better understanding of the regional and local deformation and crustal straining processes in western North America, particularly the effect of the San Andreas and nearby faults on the spatial and temporal crustal deformation behavior. Construction of theoretical models based on the mechanics of coupled elastic plate/viscoelastic foundation and large scale crack mechanics provide a rational basis for the interpretation of seismic and aseismic anomalies and expedite efforts in forecasting the stability of plate boundary deformations. In the present period, special focus is placed on the 3-D effect of irregular fault locked patches on the ground measured deformation fields. Specifically, use is made of a newly developed 3-D boundary element program to analyze the fault slip and vertical ground motion in the Parkfield area on the San Andreas.

Li, Victor C.↗

MER Field Geologic Traverse in Gusev Crater, Mars: Initial Results From the Perspective of Spirit

This report casts the initial results of the traverse and science investigations by the Mars Exploration Rover (MER) Spirit at Gusev crater [1] in terms of data sets commonly used in field geologic investigations: Local mapping of geologic features, analyses of selected samples, and their location within the local map, and the regional context of the field traverse in terms of the larger geologic and physiographic region. These elements of the field method are represented in the MER characterization of the Gusev traverse by perspective-based geologic/morphologic maps, the placement of the results from Mossbauer, APXS, Microscopic Imager, Mini-TES and Pancam multispectral studies in context within this geologic/ morphologic map, and the placement of the overall traverse in the context of narrow-angle MOC (Mars Orbiter Camera) and descent images. A major campaign over a significance fraction of the mission will be the first robotic traverse of the ejecta from a Martian impact crater along an approximate radial from the crater center. The Mars Exploration Rovers have been conceptually described as 'robotic field geologists', that is, a suite of instruments with mobility that enables far-field traverses to multiple sites located within a regional map/image base at which in situ analyses may be done. Initial results from MER, where the field geologic method has been used throughout the initial course of the investigation, confirm that this field geologic model is applicable for remote planetary surface exploration. The field geologic method makes use of near-field geologic characteristics ('outcrops') to develop an understanding of the larger geologic context through continuous loop of rational steps focused on real-time hypothesis identification and testing. This poster equates 'outcrops' with the locations of in situ investigations and 'regional context' with the geology over distance of several kilometers. Using this fundamental field geologic method, we have identified the basic local geologic materials on the floor of Gusev at this site, their compositions and likely lithologies, origins, processes that have modified these materials, and their potential significance in the interpretation of the regional geology both spatially and temporally.

Crumpler, L.↗

Dynamic testing of a two-dimensional box truss beam

Testing to determine the effects of joint freeplay and pretensioning of diagonal members on the dynamic characteristics of a two-dimensional box truss beam was conducted. The test article was ten bays of planar truss suspended by long wires at each joint. Each bay measured 2 meters per side. Pins of varying size were used to simulate various joint freeplay conditions. Single-point random excitation was the primary method of test. The rational fraction polynomial method was used to extract modal characteristics from test data. A finite element model of the test article was generated from which modal characteristics were predicted. These were compared with those obtained from tests. With the exception of the fundamental mode, correlation of theoretical and experimental results was poor, caused by the resonant coupling of local truss member bending modes with global truss beam modes. This coupling introduced many modes in the frequency range of interest whose frequencies were sensitive to joint boundary conditions. It was concluded that local/global coupling must be avoided in the frequency range where accurate modal characteristics are required.

White, Charles W.↗

Theoretical studies on the mechanical behavior of granular materials under very low intergranular stresses

The salient aspects of the theoretical modeling of a conventional triaxial test (CTC) of a cohesionless granular medium with stress and strain rate loading are described. Included are a controllable gravitational body force and provision for low confining pressure and/or very low intergranular stress. The modeling includes rational, analytic, and numerical phases, all in various stages of development. The numerical evolutions of theoretical models will be used in final design stages and in the analysis of the experimental data. In this the experimental design stage, it is of special interest to include in the candidate considerations every anomaly found in preliminary terrestrial experimentation. Most of the anomalies will be eliminated by design or enhanced for measurement as the project progresses. The main aspect of design being not the physical apparatus but the type and trajectories of loading elected. The major considerations that have been treated are: appearance and growth of local surface aberrations, stress-power coefficients, strain types, optical strain, radial bead migration, and measures of rotation for the proper stress flux.

French, Kenneth W., Jr.↗

Asymptotic theory of two-dimensional trailing-edge flows

Problems of laminar and turbulent viscous interaction near trailing edges of streamlined bodies are considered. Asymptotic expansions of the Navier-Stokes equations in the limit of large Reynolds numbers are used to describe the local solution near the trailing edge of cusped or nearly cusped airfoils at small angles of attack in compressible flow. A complicated inverse iterative procedure, involving finite-difference solutions of the triple-deck equations coupled with asymptotic solutions of the boundary values, is used to accurately solve the viscous interaction problem. Results are given for the correction to the boundary-layer solution for drag of a finite flat plate at zero angle of attack and for the viscous correction to the lift of an airfoil at incidence. A rational asymptotic theory is developed for treating turbulent interactions near trailing edges and is shown to lead to a multilayer structure of turbulent boundary layers. The flow over most of the boundary layer is described by a Lighthill model of inviscid rotational flow. The main features of the model are discussed and a sample solution for the skin friction is obtained and compared with the data of Schubauer and Klebanoff for a turbulent flow in a moderately large adverse pressure gradient.

Melnik, R. E.↗

On the breakup of viscous liquid threads

A one-dimensional model evolution equation is used to describe the nonlinear dynamics that can lead to the breakup of a cylindrical thread of Newtonian fluid when capillary forces drive the motion. The model is derived from the Stokes equations by use of rational asymptotic expansions and under a slender jet approximation. The equations are solved numerically and the jet radius is found to vanish after a finite time yielding breakup. The slender jet approximation is valid throughout the evolution leading to pinching. The model admits self-similar pinching solutions which yield symmetric shapes at breakup. These solutions are shown to be the ones selected by the initial boundary value problem, for general initial conditions. Further more, the terminal state of the model equation is shown to be identical to that predicted by a theory which looks for singular pinching solutions directly from the Stokes equations without invoking the slender jet approximation throughout the evolution. It is shown quantitatively, therefore, that the one-dimensional model gives a consistent terminal state with the jet shape being locally symmetric at breakup. The asymptotic expansion scheme is also extended to include unsteady and inerticial forces in the momentum equations to derive an evolution system modelling the breakup of Navier-Stokes jets. The model is employed in extensive simulations to compute breakup times for different initial conditions; satellite drop formation is also supported by the model and the dependence of satellite drop volumes on initial conditions is studied.

Papageorgiou, Demetrios T.↗

Portfolios in Stochastic Local Search: Efficiently Computing Most Probable Explanations in Bayesian Networks

Portfolio methods support the combination of different algorithms and heuristics, including stochastic local search (SLS) heuristics, and have been identified as a promising approach to solve computationally hard problems. While successful in experiments, theoretical foundations and analytical results for portfolio-based SLS heuristics are less developed. This article aims to improve the understanding of the role of portfolios of heuristics in SLS. We emphasize the problem of computing most probable explanations (MPEs) in Bayesian networks (BNs). Algorithmically, we discuss a portfolio-based SLS algorithm for MPE computation, Stochastic Greedy Search (SGS). SGS supports the integration of different initialization operators (or initialization heuristics) and different search operators (greedy and noisy heuristics), thereby enabling new analytical and experimental results. Analytically, we introduce a novel Markov chain model tailored to portfolio-based SLS algorithms including SGS, thereby enabling us to analytically form expected hitting time results that explain empirical run time results. For a specific BN, we show the benefit of using a homogenous initialization portfolio. To further illustrate the portfolio approach, we consider novel additive search heuristics for handling determinism in the form of zero entries in conditional probability tables in BNs. Our additive approach adds rather than multiplies probabilities when computing the utility of an explanation. We motivate the additive measure by studying the dramatic impact of zero entries in conditional probability tables on the number of zero-probability explanations, which again complicates the search process. We consider the relationship between MAXSAT and MPE, and show that additive utility (or gain) is a generalization, to the probabilistic setting, of MAXSAT utility (or gain) used in the celebrated GSAT and WalkSAT algorithms and their descendants. Utilizing our Markov chain framework, we show that expected hitting time is a rational function - i.e. a ratio of two polynomials - of the probability of applying an additive search operator. Experimentally, we report on synthetically generated BNs as well as BNs from applications, and compare SGSs performance to that of Hugin, which performs BN inference by compilation to and propagation in clique trees. On synthetic networks, SGS speeds up computation by approximately two orders of magnitude compared to Hugin. In application networks, our approach is highly competitive in Bayesian networks with a high degree of determinism. In addition to showing that stochastic local search can be competitive with clique tree clustering, our empirical results provide an improved understanding of the circumstances under which portfolio-based SLS outperforms clique tree clustering and vice versa.

Mengshoel, Ole J.↗

Five-Year Wilkinson Microwave Anisotropy Probe Observations: Data Processing, Sky Maps, and Basic Results

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the Wilkinson Microwave Anisotropy Probe (WMAP) sky survey. The new maps are consistent with previous maps and are more sensitive. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that W-band polarization data is not yet suitable for cosmological studies, but we suggest directions for further study. We do find that Ka-band data is suitable for use; in conjunction with the additional years of data, the addition of Ka band to the previously used Q- and V-band channels significantly reduces the uncertainty in the optical depth parameter, tau. Further scientific results from the five-year data analysis are presented in six companion papers and are summarized in Section 7 of this paper. With the five-year WMAP data, we detect no convincing deviations from the minimal six-parameter ACDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations in the galaxy distribution, we find (68% CL uncertainties): OMEGA(sub b)h(sup 2) = 0.02267(sup +0.00058)(sub -0.00059), OMEGA(sub c)h(sup 2) = 0.1131 plus or minus 0.0034, OMEGA(sub logical and) = 0.726 plus or minus 0.015, ns = .960 plus or minus 0.013, tau = 0.84 plus or minus 0.016, and DELTA(sup 2)(sub R) = (22.445 plus or minus 0.096) x 10(exp -9) at k = 0.002 Mpc(exp -1). From these we derive sigma(sub 8) = 0.812 plus or minus 0.026, H(sub 0) = 70.5 plus or minus 1.3 kilometers per second Mpc(exp -1), OMEGA(sub b) = 0.0456 plus or minus 0.0015, OMEGA(sub c) = .228 plus or minus 0.013, OMEGA(sub m)h(sup 2) = 0.1358(sup +0.0037)(sub -0.0036), z reion = 10.9 plus or minus 1.4, and t(sub 0) = 13.72 plus or minus 0.12 Gyr. The new limit on the tensor-to-scalar ration is r less than 0.22 (95% CL), while the evidence for a running spectral index is insignificant, dn(sub s)/d ln k = -0.028 plus or minus 0.020 (68% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and the spatial curvature of the universe: -0.14 less than 1 + w less than 0.12 (95% CL) and -0.0179 less than OMEGA(sub k) less than 0.0081 (95% CL). The number of relativistic degrees of freedom, expressed in units of the effective number of neutrino species, is found to be N(sub eff) = 4.4 plus or minus 1.5 (69% CL), consistent with the standard value of 3.04. Models with N(sub eff) = 0 are disfavored at greater than 99% confidence. Finally, new limits on physically motivated primordial non-Gaussianity parameters are -9 less than f(sup local)(sub NL) less than 111 (95% CL) and -151 less than f(sup equal)(sub NL) less than 253 (95% CL) for the local and equilateral models, respectively.

Hinshaw, G.↗

Rational Design of Nanoplasmonic Array Geometries for Biosensing

Background: Molecular diagnostics provide early and accurate diagnosis, which is essential for the prevention and treatment of infectious as well as chronic diseases. These tests are designed to detect disease-specific bioanalytes such as nucleic acid (DNA or RNA) or protein (antigens, antibodies) biomarkers. In the context of infectious disease diagnosis, nucleic acid-based detection methods are known to provide more specific and sensitive results. Here, the presence of a unique sequence belonging to the pathogenic genomic material is targeted to identify species, organism, genera and/or antimicrobial resistant gene markers. The majority of the common nucleic acid based diagnostic techniques require amplification (polymerase chain reaction, isothermal amplification etc.) of the pathogenic genetic material prior to detection impacting diagnostic speed, complexity, and cost thereby limiting ease of use. Thus, the development of simplified nucleic acid-based diagnostics that can be even used in resource-poor settings may hugely benefit patients across the globe. Nanopath is a molecular diagnostics company utilizing a solid-state nanosensor to enable sequence-specific detection of target nucleic acids without the need of amplification. These nanostructures enable ultra-sensitive biomarker detection using geometric, feature-dependent properties highly dependent on the local dielectric environment, allowing them to be sensitive to low concentration binding events. This paper describes an application of this approach to provide highly relevant clinical information within a single doctor’s office visit. Intro: The Nanopath team is in collaboration with NASA (National Aeronautics and Space Administration) and NIST (National Institute of Standards and Technology) to push the bounds of the fundamental physics associated with their biosensing platform. The ability of metals to support electromagnetic surface waves gives rise to surface plasmons when optically illuminated. This property, and its strong sensitivity to changes in the local refractive index, allows for the use of metal nanoparticles as ultra-sensitive transducers. In prior work by members of this team, ensembles of randomly oriented nanoparticles (i.e., colloidal nanorods dispersed on chip) were employed for sequence-specific nucleic acid sensing (1-3). While these particle sensors have the advantage of rapid fabrication, they suffer from low sensitivity and quality factor due to the random particle dispersity. In contrast, in this study we employ ordered array nanoparticle ensembles which can be used to improve sensor sensitivity and figure-of-merit. Study Methods Overview: In this talk, we detail the results of sensing experiments and computational simulations to outline a rational design of the structure of these plasmonic nanoparticle arrays for biomolecular sensing. Through simulation and experiment, we iteratively tailor nanostructure dimension to provide high quality signal and large resonance shifts upon modeled nucleic acid binding. In particular, full-wave electromagnetic simulations were conducted using Lumerical photonic simulation software in which periodic boundary conditions were applied in the x- and y- dimensions for each of the nanoplasmonic sensor geometries. To simulate the resonance response to changes in the bulk solution in contact with the sensor surface, the refractive index of the surrounding media was changed appropriately. Nucleic acid hybridization events were modeled using either using spherical structures approximating the relevant radius of genomic material as estimated by polymer models, or as conformal layers with the known refractive indices for nucleic acids. On the basis of initial simulations, nanosensors were fabricated using traditional electron-beam lithography protocols at NIST. To evaluate consensus between simulations and experiments, bulk sensing experiments were carried out in which the resonance peaks were obtained by submerging the sensors in refractive index standards. Key nanosensor characteristics including resonance peak locations, resonance peak shifts as a function of refractive index, and figure of merit (FOM) of extinction curves were examined between the experimental and simulation results prior to proceeding with simulations on additional geometries and more complex solution conditions, and further device fabrication. This iterative process is repeated toward a rational design of nanoplasmonic array geometries for biosensing optimizing response for targeted disease detection. In summary, this study puts forth a methodology for rational design and characterization of regularly spaced nanoparticle arrays for optics-based biosensing. The results of this study will allow for more informed design of nanostructure geometries towards sequence-specific nucleic acid detection. These improved designs have the potential to improve clinical sensitivity and limit-of-detection across disease indication.

sensor↗