Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level Set Method”

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,261 records · Page 70

Optimal Strategies for Validating TRMM Rainfall Estimates with Rain-Gauge Data

In spite of all their problems, rain gauges measure rainfall in such a direct way when compared with other methods of estimating rainfall that comparing their totals to satellite estimates remains an essential tool in the validation of satellite products. Some disagreement between averages of satellite data and rain-gauge data is expected because of the very different sampling patterns of the two systems--the satellite provides only occasional snapshots of large areas, whereas rain gauges provide continuous measurements over very small areas. A quantitative estimate of the expected level of disagreement due to the differences in sampling is needed in order for comparisons of the two kinds of averages to be informative. As part of an earlier effort to determine the sampling error of satellite averages, a space-time model for rainfall statistics was developed. Sets of parameter values for the model have been obtained for which the model statistics are found to agree quite well with those of radar data taken over the Atlantic during GATE and the western Pacific during TOGA COARE. Although the model was intended to represent the statistics of relatively large scale fluctuations of rain, it is surprisingly consistent with the very different scales on which rain gauges observe. Implications for maximizing the benefit in comparisons of rain gauge averages and TRMM estimates are derived using the model predictions.

Bell, Thomas L.↗

A Reference Model for Software and System Inspections. White Paper

Software Quality Assurance (SQA) is an important component of the software development process. SQA processes provide assurance that the software products and processes in the project life cycle conform to their specified requirements by planning, enacting, and performing a set of activities to provide adequate confidence that quality is being built into the software. Typical techniques include: (1) Testing (2) Simulation (3) Model checking (4) Symbolic execution (5) Management reviews (6) Technical reviews (7) Inspections (8) Walk-throughs (9) Audits (10) Analysis (complexity analysis, control flow analysis, algorithmic analysis) (11) Formal method Our work over the last few years has resulted in substantial knowledge about SQA techniques, especially the areas of technical reviews and inspections. But can we apply the same QA techniques to the system development process? If yes, what kind of tailoring do we need before applying them in the system engineering context? If not, what types of QA techniques are actually used at system level? And, is there any room for improvement.) After a brief examination of the system engineering literature (especially focused on NASA and DoD guidance) we found that: (1) System and software development process interact with each other at different phases through development life cycle (2) Reviews are emphasized in both system and software development. (Figl.3). For some reviews (e.g. SRR, PDR, CDR), there are both system versions and software versions. (3) Analysis techniques are emphasized (e.g. Fault Tree Analysis, Preliminary Hazard Analysis) and some details are given about how to apply them. (4) Reviews are expected to use the outputs of the analysis techniques. In other words, these particular analyses are usually conducted in preparation for (before) reviews. The goal of our work is to explore the interaction between the Quality Assurance (QA) techniques at the system level and the software level.

He, Lulu↗

An Explainable Machine-Learning Model for Compensatory Reserve Measurement: Methods for Feature Selection and the Effects of Subject Variability

Tracking vital signs accurately is critical for triaging a patient and ensuring timely therapeutic intervention. The patient’s status is often clouded by compensatory mechanisms that can mask injury severity. The compensatory reserve measurement (CRM) is a triaging tool derived from an arterial waveform that has been shown to allow for earlier detection of hemorrhagic shock. However, the deep-learning artificial neural networks developed for its estimation do not explain how specific arterial waveform elements lead to predicting CRM due to the large number of parameters needed to tune these models. Alternatively, we investigate how classical machine-learning models driven by specific features extracted from the arterial waveform can be used to estimate CRM. More than 50 features were extracted from human arterial blood pressure data sets collected during simulated hypovolemic shock resulting from exposure to progressive levels of lower body negative pressure. A bagged decision tree design using the ten most significant features was selected as optimal for CRM estimation. This resulted in an average root mean squared error in all test data of 0.171, similar to the error for a deep-learning CRM algorithm at 0.159. By separating the dataset into sub-groups based on the severity of simulated hypovolemic shock withstood, large subject variability was observed, and the key features identified for these sub-groups differed. This methodology could allow for the identification of unique features and machine-learning models to differentiate individuals with good compensatory mechanisms against hypovolemia from those that might be poor compensators, leading to improved triage of trauma patients and ultimately enhancing military and emergency medicine.

60 APPLIED LIFE SCIENCES↗

Network-Level Traffic Signal Cooperation: A Higher-Order Conflict Graph Approach

Traffic signal control and cooperation are extremely important to alleviate traffic congestion in a large traffic network. This study develops a higher-order conflict graph approach for network-wide traffic signal control and cooperation. A conflict graph is applied to model the traffic signal configurations, which identifies the conflict and unconflicted movements for each intersection. In conflict graph, the node represents each movement. The weight of each node can be defined as traffic volume, queue length, fuel consumption, or any weighted combinations of these measurements. The calculation of the optimal green light duration and green light sequence (for different movements) is equivalent to sequentially finding the maximum weight independent set (MWIS) in the conflict graph. The conflict graph also provides a uniform and efficient way to connect traffic signal operations among nearby intersections spatially. Then, we introduced the concept of the k -th order neighborhood to model the degree of connectivity between each movement to the movements at upstream or downstream intersections. The weight of each node in the higher-order conflict graph not only represents its own congestion level, but also relates to the traffic conditions of nearby intersections. Through this approach, the cooperation of multiple intersections can be realized by incorporating their spatial connectivity into conflict graph and solving the MWIS problem. A simulation network is built in SUMO to test the effectiveness of the proposed method. Results suggested that the proposed model outperformed other state-of-the-art signal control methods. Also, the scheme maintains good performance under varying traffic demands.

42 ENGINEERING↗

Linear quadratic servo control of a reusable rocket engine

The paper deals with the development of a design method for a servo component in the frequency domain using singular values and its application to a reusable rocket engine. A general methodology used to design a class of linear multivariable controllers for intelligent control systems is presented. Focus is placed on performance and robustness characteristics, and an estimator design performed in the framework of the Kalman-filter formalism with emphasis on using a sensor set different from the commanded values is discussed. It is noted that loop transfer recovery modifies the nominal plant noise intensities in order to obtain the desired degree of robustness to uncertainty reflected at the plant input. Simulation results demonstrating the performance of the linear design on a nonlinear engine model over all power levels during mainstage operation are discussed.

Musgrave, Jeffrey L.↗

A Gold Standards Approach to Training Instructors to Evaluate Crew Performance

The Advanced Qualification Program requires that airlines evaluate crew performance in Line Oriented Simulation. For this evaluation to be meaningful, instructors must observe relevant crew behaviors and evaluate those behaviors consistently and accurately against standards established by the airline. The airline industry has largely settled on an approach in which instructors evaluate crew performance on a series of event sets, using standardized grade sheets on which behaviors specific to event set are listed. Typically, new instructors are given a class in which they learn to use the grade sheets and practice evaluating crew performance observed on videotapes. These classes emphasize reliability, providing detailed instruction and practice in scoring so that all instructors within a given class will give similar scores to similar performance. This approach has value but also has important limitations; (1) ratings within one class of new instructors may differ from those of other classes; (2) ratings may not be driven primarily by the specific behaviors on which the company wanted the crews to be scored; and (3) ratings may not be calibrated to company standards for level of performance skill required. In this paper we provide a method to extend the existing method of training instructors to address these three limitations. We call this method the "gold standards" approach because it uses ratings from the company's most experienced instructors as the basis for training rater accuracy. This approach ties the training to the specific behaviors on which the experienced instructors based their ratings.

Baker, David P.↗

Locomotion in simulated microgravity: gravity replacement loads

BACKGROUND: When an astronaut walks or runs on a treadmill in microgravity, a subject load device (SLD) is used to return him or her back to the treadmill belt. The gravity replacement load (GRL) in the SLD is transferred, via a harness, to the pelvis and/or the shoulders. This research compared comfort and ground reaction forces during treadmill running in a microgravity locomotion simulator at GRLs of 60%, 80%, and 100% of body weight (BW). Two harness designs (shoulder springs only (SSO) and waist and shoulder springs (WSS)) were used. HYPOTHESES: 1) The 100% BW gravity replacement load conditions would be comfortably tolerated and would result in larger ground reaction forces and loading rates than the lower load conditions, and 2) the WSS harness would be more comfortable than the SSO harness. METHODS: Using the Penn State Zero Gravity Locomotion Simulator (ZLS), 8 subjects ran at 2.0 m x s(-1) (4.5 mph) for 3 min at each GRL setting in each harness. Subjective ratings of harness comfort, ground reaction forces, and GRL data were collected during the final minute of exercise. RESULTS: The 100% BW loading conditions were comfortably tolerated (2.3 on a scale of 0-10), although discomfort increased as the GRL increased. There were no overall differences in perceived comfort between the two harnesses. The loading rates (27.1, 33.8, 39.1 BW x s(-1)) and the magnitudes of the first (1.0, 1.4, 1.6 BW) and second (1.3, 1.7, 1.9 BW) peaks of the ground reaction force increased with increasing levels (60, 80, 100% BW respectively) of GRL. CONCLUSIONS: Subjects were able to tolerate a GRL of 100% BW well. The magnitude of the ground reaction force peaks and the loading rate is directly related to the magnitude of the GRL.

NASA Discipline Musculoskeletal↗

Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays

Single-slope analog-to-digital converters (ADCs) are particularly useful for onchip digitization in focal plane arrays (FPAs) because of their inherent monotonicity, relative simplicity, and efficiency for column-parallel applications, but they are comparatively slow. Squareroot encoding can allow the number of code values to be reduced without loss of signal-to-noise ratio (SNR) by keeping the quantization noise just below the signal shot noise. This encoding can be implemented directly by using a quadratic ramp. The reduction in the number of code values can substantially increase the quantization speed. However, in an FPA, the fixed pattern noise (FPN) limits the use of small quantization steps at low signal levels. If the zero-point is adjusted so that the lowest column is onscale, the other columns, including those at the center of the distribution, will be pushed up the ramp where the quantization noise is higher. Additionally, the finite frequency response of the ramp buffer amplifier and the comparator distort the shape of the ramp, so that the effective ramp value at the time the comparator trips differs from the intended value, resulting in errors. Allowing increased settling time decreases the quantization speed, while increasing the bandwidth increases the noise. The FPN problem is solved by breaking the ramp into two portions, with some fraction of the available code values allocated to a linear ramp and the remainder to a quadratic ramp. To avoid large transients, both the value and the slope of the linear and quadratic portions should be equal where they join. The span of the linear portion must cover the minimum offset, but not necessarily the maximum, since the fraction of the pixels above the upper limit will still be correctly quantized, albeit with increased quantization noise. The required linear span, maximum signal and ratio of quantization noise to shot noise at high signal, along with the continuity requirement, determines the number of code values that must be allocated to each portion. The distortion problem is solved by using a lookup table to convert captured code values back to signal levels. The values in this table will be similar to the intended ramp value, but with a correction for the finite bandwidth effects. Continuous-time comparators are used, and their bandwidth is set below the step rate, which smoothes the ramp and reduces the noise. No settling time is needed, as would be the case for clocked comparators, but the low bandwidth enhances the distortion of the non-linear portion. This is corrected by use of a return lookup table, which differs from the one used to generate the ramp. The return lookup table is obtained by calibrating against a stepped precision DC reference. This results in a residual non-linearity well below the quantization noise. This method can also compensate for differential non-linearity (DNL) in the DAC used to generate the ramp. The use of a ramp with a combination of linear and quadratic portions for a single-slope ADC is novel. The number of steps is minimized by keeping the step size just below the photon shot noise. This in turn maximizes the speed of the conversion. High resolution is maintained by keeping small quantization steps at low signals, and noise is minimized by allowing the lowest analog bandwidth, all without increasing the quantization noise. A calibrated return lookup table allows the system to maintain excellent linearity.

Wrigley, Chris J.↗

Excited-state uncertainties in lattice-QCD calculations of multi-hadron systems

Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.

Detmold, William [MIT, Cambridge, CTP]↗

Operator-level quantum acceleration of non-logconcave sampling

Sampling from probability distributions of the form 𝝈 ∝ e −𝜷V , where V is a continuous potential, is a fundamental task across physics, chemistry, biology, computer science, and statistics. However, when V is nonconvex, the resulting distribution becomes non-logconcave, and classical methods such as Langevin dynamics often exhibit poor performance. We introduce a quantum algorithm that provably accelerates a broad class of continuous-time sampling dynamics. For Langevin dynamics, our method encodes the target Gibbs measure into the amplitudes of aquantum state, identified as the kernel of a block matrix derived from a factorization of the Witten Laplacian operator. This connection enables Gibbs sampling via singular value thresholding and yields up to a quartic quantum speedup over best-knownclassical Langevin-based methods in the non-logconcave setting. Building on this framework, we further develop the first quantum algorithm that accelerates replica exchange Langevin diffusion, a widely used method for sampling from complex, rugged energy landscapes.

97 MATHEMATICS AND COMPUTING↗

A Contextually Supervised Optimization-Based HVAC Load Disaggregation Methodology

This paper presents a novel contextually supervised optimization-based approach for disaggregating heating, ventilation, and air-conditioning (HVAC) loads using smart meter or Supervisory Control and Data Acquisition data. To disaggregate the load into HVAC loads, large and infrequently used loads (LIUL), and base loads, we formulate an optimization problem to minimize a set of five loss terms, consisting of the reconstruction errors of the overall load profile, the ramp rate losses, and three distinct loss functions linked with the HVAC load, base load, and LIUL, respectively. To enhance accuracy, we incorporate two forms of contextual information into the problem formulation. First, we utilize mutual information to estimate HVAC energy consumption. Second, we employ a base load dictionary to constrain HVAC load estimation errors. The obtained HVAC load profiles are fine-tuned by abnormal ramp detection followed by binary hypothesis testing. Here, the proposed method is developed and tested using sub-metered residential and commercial building data. Simulation results show that the proposed method outperforms existing methods across various data resolutions and load aggregation levels, showing excellent transferability and generalizability.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Variable mixer propulsion cycle

A design technique, method and apparatus are delineated for controlling the bypass gas stream pressure and varying the bypass ratio of a mixed flow gas turbine engine in order to achieve improved performance. The disclosed embodiments each include a mixing device for combining the core and bypass gas streams. The variable area mixing device permits the static pressures of the core and bypass streams to be balanced prior to mixing at widely varying bypass stream pressure levels. The mixed flow gas turbine engine therefore operates efficiently over a wide range of bypass ratios and the dynamic pressure of the bypass stream is maintained at a level which will keep the engine inlet airflow matched to an optimum design level throughout a wide range of engine thrust settings.

Rundell, D. J.↗

Gas temperature measurements using the dual-line detection Rayleigh scattering technique

A new laser-induced Rayleigh scattering method is presented for the improved temperature diagnostics of gas flows. In the present technique, the two lines of a copper vapor laser are used to obtain the time and space resolved temperature. A single set of optics is used to form the optical probe and to collect the signal simultaneously from both the 510 nm and the 578 nm lines. The dual-line detection allows for the determination and removal of surface-scattered laser light from a Rayleigh signal thereby improving the applicability of Rayleigh scattering to near wall flows with a high degree of glare. An optical system using the dual-line detection technique is built, calibrated and tested in a hot air jet under various levels of background contamination. The results indicate that highly accurate temperature measurements are possible even when the laser-line background intensity, captured by the collecting optics, is five times that of the Rayleigh signal.

Otugen, M. Volkan↗

Gas temperature measurements using the dual-line detection Rayleigh scattering technique

A new laser-induced Rayleigh scattering method is presented for the improved temperature diagnostics of gas flows. In the present technique, the two lines of a copper vapor laser are used to obtain the time and space resolved temperature. A single set of optics is used to form the optical probe and to collect the signal simultaneously from both the 510 nm and the 578 nm lines. The dual-line detection allows for the determination and removal of surface-scattered laser light from a Rayleigh signal thereby improving the applicability of Rayleigh scattering to near wall flows with a high degree of glare. An optical system using the dual-line detection technique is built, calibrated and tested in a hot air jet under various levels of background contamination. The results indicate that highly accurate temperature measurements are possible even when the laser-line background intensity, captured by the collecting optics, is five times that of the Rayleigh signal.

Otugen, M. V.↗

Relationship Between Ecosystem Productivity and Photosynthetically Active Radiation for Northern Peatlands

We analyzed the relationship between net ecosystem exchange of carbon dioxide (NEE) and irradiance (as photosynthetic photon flux density or PPFD), using published and unpublished data that have been collected during midgrowing season for carbon balance studies at seven peatlands in North America and Europe, NEE measurements included both eddy-correlation tower and clear, static chamber methods, which gave very similar results. Data were analyzed by site, as aggregated data sets by peatland type (bog, poor fen, rich fen, and all fens) and as a single aggregated data set for all peatlands. In all cases, a fit with a rectangular hyperbola (NEE = alpha PPFD P(sub max)/(alpha PPFD + P(sub max) + R) better described the NEE-PPFD relationship than did a linear fit (NEE = beta PPFD + R). Poor and rich fens generally had similar NEE-PPFD relationships, while bogs had lower respiration rates (R = -2.0 micro mol m(exp -2) s(exp -1) for bogs and -2.7 micro mol m(exp -2) s(exp -1)) for fens) and lower NEE at moderate and high light levels (P(sub max)= 5.2 micro mol m(exp -2) s(exp -1) for bogs and 10.8 micro mol m(exp -2) s(exp -1) for fens). As a single class, northern peatlands had much smaller ecosystem respiration (R = -2.4 micro mol m(exp -2) s(exp -1)) and NEE rates (alpha = 0.020 and P(sub max)= 9.2 micro mol m(exp -2) s(exp -1)) than the upland ecosystems (closed canopy forest, grassland, and cropland). Despite this low productivity, northern peatland soil carbon pools are generally 5-50 times larger than upland ecosystems because of slow rates of decomposition caused by litter quality and anaerobic, cold soils.

Frolking, S. E.↗

Determination of muscle-specific glucose flux using radioactive stereoisomers and microdialysis

The purpose of the present study was to evaluate a novel approach for determining skeletal muscle-specific glucose flux using radioactive stereoisomers and the microdialysis technique. Microdialysis probes were inserted into the vastus lateralis muscle of human subjects and perfused (4 microl/min) with a Ringer solution containing small amounts of radioactive D- and L-glucose as the internal reference markers for determining probe recovery as well as varying concentrations of insulin (0-10 microM). The rationale behind this approach was that both stereoisomers would be equally affected by the factors that determine probe recovery, with the exception of L-glucose, which is nonmetabolizable and would not be influenced by tissue uptake. Therefore, any differences in the probe recovery ratios between the D- and L-stereoisomers represent changes in skeletal muscle glucose uptake directly at the tissue level. There were no differences in probe recovery between the D- (42.3 +/- 3.5%) and L- (41.2 +/- 3.5) stereoisomers during the control period (no insulin), which resulted in a D/L ratio of 1.04 +/- 0.03. However, during insulin perfusion (1 microM), The D/L ratio increased to 1.62 +/- 0.08 and 1.58 +/- 0.07 (P < 0.05) during the two collection (0-15 and 15-30 min) periods, respectively. This was accomplished solely by an increase (P < 0.05) in D-glucose probe recovery, as L-glucose probe recovery remained unchanged. In a second set of experiments, the perfusion of 10 microM insulin did not increase the D/L ratio (1.40 +/- 0.11) above that observed during 1.0 microM (1.41 +/- 0.07) insulin perfusion. These data suggest that this method is sufficiently sensitive to detect differences in insulin-stimulated glucose uptake; thus the use of radioactive stereoisomers in conjunction with the microdialysis technique provides a novel and useful technique for determining tissue-specific glucose flux and insulin sensitivity.

Non-NASA Center↗

First application of a liquid argon time projection chamber for the search for intranuclear neutron-antineutron transitions and annihilation in 40 Ar using the MicroBooNE detector

We present a novel methodology to search for intranuclear neutron-antineutron transition (n→n̅) followed by n̅-nucleon annihilation within an 40 Ar nucleus, using the MicroBooNE liquid argon time projection chamber (LArTPC) detector. A discovery of n→n̅ transition or a new best limit on the lifetime of this process would either constitute physics beyond the Standard Model or greatly constrain theories of baryogenesis, respectively. The approach presented in this paper makes use of deep learning methods to select n→n̅ events based on their unique features and differentiate them from cosmogenic backgrounds. The achieved signal and background efficiencies are (70.22 ± 6.04)% and (0.0020 ± 0.0003)%, respectively. A demonstration of a search is performed with a data set corresponding to an exposure of 3.32 ×10 26 neutron-years, and where the background rate is constrained through direct measurement, assuming the presence of a negligible signal. With this approach, no excess of events over the background prediction is observed, setting a demonstrative lower bound on the n→n̅ lifetime in 40 Ar of τ m ≳ 1.1×10 26 years, and on the free n→n̅ transition time of τ n→n̅ ≳ 2.6×10 5 s, each at the 90% confidence level. This analysis represents a first-ever proof-of-principle demonstration of the ability to search for this rare process in LArTPCs with high efficiency and low background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A consensus-based ensemble approach to improve transcriptome assembly

Systems-level analyses, such as differential gene expression analysis, co-expression analysis, and metabolic pathway reconstruction, depend on the accuracy of the transcriptome. Multiple tools exist to perform transcriptome assembly from RNAseq data. However, assembling high quality transcriptomes is still not a trivial problem. This is especially the case for non-model organisms where adequate reference genomes are often not available. Different methods produce different transcriptome models and there is no easy way to determine which are more accurate. Furthermore, having alternative-splicing events exacerbates such difficult assembly problems. While benchmarking transcriptome assemblies is critical, this is also not trivial due to the general lack of true reference transcriptomes. In this study, we first provide a pipeline to generate a set of the simulated benchmark transcriptome and corresponding RNAseq data. Using the simulated benchmarking datasets, we compared the performance of various transcriptome assembly approaches including both de novo and genome-guided methods. The results showed that the assembly performance deteriorates significantly when alternative transcripts (isoforms) exist or for genome-guided methods when the reference is not available from the same genome. To improve the transcriptome assembly performance, leveraging the overlapping predictions between different assemblies, we present a new consensus-based ensemble transcriptome assembly approach, ConSemble. Without using a reference genome, ConSemble using four de novo assemblers achieved an accuracy up to twice as high as any de novo assemblers we compared. When a reference genome is available, ConSemble using four genome-guided assemblies removed many incorrectly assembled contigs with minimal impact on correctly assembled contigs, achieving higher precision and accuracy than individual genome-guided methods. Furthermore, ConSemble using de novo assemblers matched or exceeded the best performing genome-guided assemblers even when the transcriptomes included isoforms. We thus demonstrated that the ConSemble consensus strategy both for de novo and genome-guided assemblers can improve transcriptome assembly. The RNAseq simulation pipeline, the benchmark transcriptome datasets, and the script to perform the ConSemble assembly are all freely available from: http://bioinfolab.unl.edu/emlab/consemble/.

59 BASIC BIOLOGICAL SCIENCES↗