Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parametric processes”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

True Concurrent Thermal Engineering Integrating CAD Model Building with Finite Element and Finite Difference Methods

Thermal engineering has long been left out of the concurrent engineering environment dominated by CAD (computer aided design) and FEM (finite element method) software. Current tools attempt to force the thermal design process into an environment primarily created to support structural analysis, which results in inappropriate thermal models. As a result, many thermal engineers either build models "by hand" or use geometric user interfaces that are separate from and have little useful connection, if any, to CAD and FEM systems. This paper describes the development of a new thermal design environment called the Thermal Desktop. This system, while fully integrated into a neutral, low cost CAD system, and which utilizes both FEM and FD methods, does not compromise the needs of the thermal engineer. Rather, the features needed for concurrent thermal analysis are specifically addressed by combining traditional parametric surface based radiation and FD based conduction modeling with CAD and FEM methods. The use of flexible and familiar temperature solvers such as SINDA/FLUINT (Systems Improved Numerical Differencing Analyzer/Fluid Integrator) is retained.

Panczak, Tim↗
NASA NTRS2003

Polar Heat Flow on Io

Recently, Galileo spacecraft data have revealed Io's polar regions to be much warmer than previously expected. This unexpected development came from Photo-Polarimeter Radiometer (PPR) data which show that the minimum night temperatures are in the range of 90-95 K virtually everywhere on Io. The minimum night temperatures show no dependence upon latitude and, when away from the sunset terminator, they show no dependence upon time of night. This is indeed bizarre behavior for surface units which generally had been assumed to be passive with respect to Io's pervasive volcanism. Night temperatures of 90-95 K at high, polar latitudes are particularly hard to explain. Even assuming infinite thermal inertia, at these latitudes there is insufficient sunlight to support these warm night temperatures. Thus, through the process of elimination of other possibilities, we come to the conclusion that these surfaces are volcanically heated. Taking previously passive units and turning them into new sources of heat flow is a radical departure from previous thermophysical model paradigms. However, the geological interpretation is straight forward. We are simply seeing the effect of old, cool lava flows which cover most of the surface of Io but yet have some heat to radiate. Under these new constraints, we have taken on the challenge of formulating a physical model which quantitatively reproduces all of the observations of Io's thermal emission. In the following we introduce a new parametric model which suffices to identify a previously unrecognized polar component of Io's heat flow.

Veeder, G. J.↗
NASA NTRS1995

A probabilistic model of a porous heat exchanger

This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Agrawal, O. P.↗
NASA NTRS2006

Shuttle Damage/Repair from the Perspective of Hypersonic Boundary Layer Transition - Experimental Results

An overview is provided of the experimental wind tunnel program conducted at the NASA Langley Research Center Aerothermodynamics Laboratory in support of an agency-wide effort to prepare the Shuttle Orbiter for Return-to-Flight. The effect of an isolated protuberance and an isolated rectangular cavity on hypersonic boundary layer transition onset on the windward surface of the Shuttle Orbiter has been experimentally characterized. These experimental studies were initiated to provide a protuberance and cavity effects database for developing hypersonic transition criteria to support on-orbit disposition of thermal protection system damage or repair. In addition, a synergistic experimental investigation was undertaken to assess the impact of an isolated mass-flow entrainment source (simulating pyrolysis/outgassing from a proposed tile repair material) on boundary layer transition. A brief review of the relevant literature regarding hypersonic boundary layer transition induced from cavities and localized mass addition from ablation is presented. Boundary layer transition results were obtained using 0.0075-scale Orbiter models with simulated tile damage (rectangular cavities) of varying length, width, and depth and simulated tile damage or repair (protuberances) of varying height. Cavity and mass addition effects were assessed at a fixed location (x/L = 0.3) along the model centerline in a region of near zero pressure gradient. Cavity length-to-depth ratio was systematically varied from 2.5 to 17.7 and length-to-width ratio of 1 to 8.5. Cavity depth-to-local boundary layer thickness ranged from 0.5 to 4.8. Protuberances were located at several sites along the centerline and port/starboard attachment lines along the chine and wing leading edge. Protuberance height-to-boundary layer thickness was varied from approximately 0.2 to 1.1. Global heat transfer images and heating distributions of the Orbiter windward surface using phosphor thermography were used to infer the state of the boundary layer (laminar, transitional, or turbulent). Test parametrics include angles-of-attack of 30 deg and 40 deg, sideslip angle of 0 deg, freestream Reynolds numbers from 0.02x106 to 7.3x106 per foot, edge-to-wall temperature ratio from 0.4 to 0.8, and normal shock density ratios of approximately 5.3, 6.0, and 12 in Mach 6 air, Mach 10 air, and Mach 6 CF4, respectively. Testing to simulate the effects of ablation from a proposed tile repair concept indicated that transition was not a concern. The experimental protuberance and cavity databases highlighted in this report were used to formulate boundary layer transition correlations that were an integral part of an analytical process to disposition observed Orbiter TPS damage during STS- 114.

Horvath, Thomas J.↗
NASA NTRS2009

NASA Tech Briefs, March 2009

Topics covered include: Improved Instrument for Detecting Water and Ice in Soil; Real-Time Detection of Dust Devils from Pressure Readings; Determining Surface Roughness in Urban Areas Using Lidar Data; DSN Data Visualization Suite; Hamming and Accumulator Codes Concatenated with MPSK or QAM; Wide-Angle-Scanning Reflectarray Antennas Actuated by MEMS; Biasable Subharmonic Membrane Mixer for 520 to 600 GHz; Hardware Implementation of Serially Concatenated PPM Decoder; Symbolic Processing Combined with Model-Based Reasoning; Presentation Extensions of the SOAP; Spreadsheets for Analyzing and Optimizing Space Missions; Processing Ocean Images to Detect Large Drift Nets; Alternative Packaging for Back-Illuminated Imagers; Diamond Machining of an Off-Axis Biconic Aspherical Mirror; Laser Ablation Increases PEM/Catalyst Interfacial Area; Damage Detection and Self-Repair in Inflatable/Deployable Structures; Polyimide/Glass Composite High-Temperature Insulation; Nanocomposite Strain Gauges Having Small TCRs; Quick-Connect Windowed Non-Stick Penetrator Tips for Rapid Sampling; Modeling Unsteady Cavitation and Dynamic Loads in Turbopumps; Continuous-Flow System Produces Medical-Grade Water; Discrimination of Spore-Forming Bacilli Using spoIVA; nBn Infrared Detector Containing Graded Absorption Layer; Atomic References for Measuring Small Accelerations; Ultra-Broad-Band Optical Parametric Amplifier or Oscillator; Particle-Image Velocimeter Having Large Depth of Field; Enhancing SERS by Means of Supramolecular Charge Transfer; Improving 3D Wavelet-Based Compression of Hyperspectral Images; Improved Signal Chains for Readout of CMOS Imagers; SOI CMOS Imager with Suppression of Cross-Talk; Error-Rate Bounds for Coded PPM on a Poisson Channel; Biomorphic Multi-Agent Architecture for Persistent Computing; and Using Covariance Analysis to Assess Pointing Performance.

Source record↗
NASA NTRS1994

Physical processes in eclipsing pulsars: Eclipse mechanisms and diagnostics

We investigate how the radio emission of a pulsar interacts with plasma derived from a stellar companion. Various physical mechanisms that can cause radio pulse eclipse are discussed, and predictions are made for the polarization properties of the emergent radio wave. We consider eclipses by a wind from the stellar companion, by a stellar magnetosphere, or by material entrained in the pulsar wind. Eclipses due to refraction require either a relatively high plasma density or a sharp edge to the plasma distribution. The conditions that must prevail for free-free absorption to be effective in eclipsing a radio beam are also outlined. Pulse smearing may be important at higher frequencies; related eclipse mechanisms include pulse spreading due to a rapidly changing electron column, and scattering by Langmuir turbulence. The high brightness temperature radio beam can generate its own plasma turbulence via a number of nonlinear parametric instabilities, such as the instability associated with stimulated Raman scattering. When the plasma turbulence is heavily damped, the radio bean can still undergo induced Compton scattering. Stimulated scattering effects such as these are very sensitive to the presence of narrow-band substructure in the pulsar radio emission. Finally, we consider the possibility that plasma derived from a stellar companion may mix with the relativistic pulsar wind and cause cyclotron absorption at low radio frequencies. Even if the cyclotron optical depth is small, fluctuations in the emergent polarization of the radio beam on the timescale of a few seconds are a very sensitive probe of the spatial structure of the magnetic field in the pulsar wind. The current observational properties of two known eclipsing pulsar systems, PSR 1957+20 and PSR 1744-24A, are used to construct tentative eclipse models. The favored model for PSR 1957+20 is cyclotron or synchrotron absorption by plasma embedded in the pulsar wind combined with pulse smearing at high frequency, and the favored model for PSR 1744-24A is backscattering off plasma turbulence generated by the stimulated Raman scattering parametric instability. Pulsar eclipses promise to provide a good diagnostic of pulsar winds and possible of the pulse emission mechanism.

Thompson, C.↗
NASA NTRS1996

Signal Processing and Calibration of Continuous-Wave Focused CO2 Doppler Lidars for Atmospheric Backscatter Measurement

Two continuous-wave(CW)focused C02 Doppler lidars (9.1 and 10.6 micrometers) were developed for airborne in situ aerosol backscatter measurements. The complex path of reliably calibrating these systems, with different signal processors, for accurate derivation of atmospheric backscatter coefficients is documented. Lidar calibration for absolute backscatter measurement for both lidars is based on range response over the lidar sample volume, not solely at focus. Both lidars were calibrated with a new technique using well-characterized aerosols as radiometric standard targets and related to conventional hard-target calibration. A digital signal processor (DSP), a surface acoustic and spectrum analyzer and manually tuned spectrum analyzer signal analyzers were used. The DSP signals were analyzed with an innovative method of correcting for systematic noise fluctuation; the noise statistics exhibit the chi-square distribution predicted by theory. System parametric studies and detailed calibration improved the accuracy of conversion from the measured signal-to-noise ratio to absolute backscatter. The minimum backscatter sensitivity is approximately 3 x 10(exp -12)/m/sr at 9.1 micrometers and approximately 9 x 10(exp -12)/m/sr at 10.6 micrometers. Sample measurements are shown for a flight over the remote Pacific Ocean in 1990 as part of the NASA Global Backscatter Experiment (GLOBE) survey missions, the first time to our knowledge that 9.1-10.6 micrometer lidar intercomparisons were made. Measurements at 9.1 micrometers, a potential wavelength for space-based lidar remote-sensing applications, are to our knowledge the first based on the rare isotope C-12 O(2)-18 gas.

Rothermel, Jeffry↗
NASA NTRSDate not supplied

Characterization of Reactant Refill and Detonation Wave Dynamics in a GOx/Natural-gas RDRE Using Simultaneous High Repetition-Rate OH-PLIF and Chemiluminescence

The potential application of rotating detonation engines (RDEs) in rocket combustors hinges on a fundamental understanding of detonation wave structure and injector characteristics with fuel and oxidizer compositions relevant in rocket systems. Simultaneous 300 kHz-rate broadband OH* chemiluminescence and OH-PLIF imaging is employed in a fully optically accessible Natural Gas-GOx rotating detonation rocket engine (RDRE) to visualize reactant refill dynamics and detonation wave structure. A custom-built KTP-type optical parametric oscillator (OPO) is coupled with a nanosecond high-repetition-rate burst-mode laser to output284 nm light and target excitation of the Q1(9) transition in the OH radical. Significant deflagrative burning is observed throughout the chamber as a consequence of the oxygen-rich environment. Trailing Azimuthal Reflected Shock Combustion (ARSC) system, similar to those in a H2-air RDE are observed, burning unburned reactants in the region immediately following the primary detonation wave. Contact burning, as indicated in this study, does not seem to be a primary loss mechanism. The simultaneous measurement of OH and OH* show that axial locations exist in the refill process where OH radicals are present, and produced due to shear layer induced deflagration, however, these zones do not produce excited state OH*. While a deeper understanding of the underlying physics in RDRE systems requires further investigation, this work highlights a first-of-its-kind visualization of the turbulent combustion product field and reactant refill characteristics in this highly unsteady environment.

Propulsion↗
NASA NTRS2005

Integrated System-Level Optimization for Concurrent Engineering With Parametric Subsystem Modeling

The introduction of concurrent design practices to the aerospace industry has greatly increased the productivity of engineers and teams during design sessions as demonstrated by JPL's Team X. Simultaneously, advances in computing power have given rise to a host of potent numerical optimization methods capable of solving complex multidisciplinary optimization problems containing hundreds of variables, constraints, and governing equations. Unfortunately, such methods are tedious to set up and require significant amounts of time and processor power to execute, thus making them unsuitable for rapid concurrent engineering use. This paper proposes a framework for Integration of System-Level Optimization with Concurrent Engineering (ISLOCE). It uses parametric neural-network approximations of the subsystem models. These approximations are then linked to a system-level optimizer that is capable of reaching a solution quickly due to the reduced complexity of the approximations. The integration structure is described in detail and applied to the multiobjective design of a simplified Space Shuttle external fuel tank model. Further, a comparison is made between the new framework and traditional concurrent engineering (without system optimization) through an experimental trial with two groups of engineers. Each method is evaluated in terms of optimizer accuracy, time to solution, and ease of use. The results suggest that system-level optimization, running as a background process during integrated concurrent engineering sessions, is potentially advantageous as long as it is judiciously implemented.

Schuman, Todd↗
NASA NTRS2002

A geometric deformable model for echocardiographic image segmentation

Gradient vector flow (GVF), an elegant external force for parametric deformable models, can capture object boundaries from both sides. A new geometric deformable model is proposed that combines GVF and the geodesic active contour model. The level set method is used as the numerical method of this model. The model is applied for echocardiographic image segmentation.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS1990

The application of CFD to rotary wing flow problems

Rotorcraft aerodynamics is especially rich in unsolved problems, and for this reason the need for independent computational and experimental studies is great. Three-dimensional unsteady, nonlinear potential methods are becoming fast enough to enable their use in parametric design studies. At present, combined CAMRAD/FPR analyses for a complete trimmed rotor soltution can be performed in about an hour on a CRAY Y-MP (or ten minutes, with multiple processors). These computational speeds indicate that in the near future many of the large CFD problems will no longer require a supercomputer. The ability to convect circulation is routine for integral methods, but only recently was it discovered how to do the same with differential methods. It is clear that the differential CFD rotor analyses are poised to enter the engineering workplace. Integral methods already constitute a mainstay. Ultimately, it is the users who will integrate CFD into the entire engineering process and provide a new measure of confidence in design and analysis. It should be recognized that the above classes of analyses do not include several major limiting phenomena which will continue to require empirical treatment because of computational time constraints and limited physical understanding. Such empirical treatment should be included, however, into the developing CFD, engineering level analyses. It is likely that properly constructed flow models containing corrections from physical testing will be able to fill in unavoidable gaps in the experimental data base, both for basic studies and for specific configuration testing. For these kinds of applications, computational cost is not an issue. Finally, it should be recognized that although rotorcraft are probably the most complex of aircraft, the rotorcraft engineering community is very small compared to the fixed-wing community. Likewise, rotorcraft CFD resources can never achieve fixed-wing proportions and must be used wisely. Therefore the fixed-wing work must be gleaned for many of the basic methods.

Caradonna, F. X.↗
NASA NTRS1992

Theoretical studies of the physics of the solar atmosphere

Significant advances in our theoretical basis for understanding several physical processes related to dynamical phenomena on the sun were achieved. We have advanced a new model for spicules and fibrils. We have provided a simple physical view of resonance absorption of MHD surface waves; this allowed an approximate mathematical procedure for obtaining a wealth of new analytical results which we applied to coronal heating and p-mode absorption at magnetic regions. We provided the first comprehensive models for the heating and acceleration of the transition region, corona, and solar wind. We provided a new view of viscosity under coronal conditions. We provided new insights into Alfven wave propagation in the solar atmosphere. And recently we have begun work in a new direction: parametric instabilities of Alfven waves.

Hollweg, Joseph V.↗
NASA NTRS1994

Towards a Comprehensive Computational Simulation System for Turbomachinery

The objective of this work is to develop algorithms associated with a comprehensive computational simulation system for turbomachinery flow fields. This development is accomplished in a modular fashion. These modules includes grid generation, visualization, network, simulation, toolbox, and flow modules. An interactive grid generation module is customized to facilitate the grid generation process associated with complicated turbomachinery configurations. With its user-friendly graphical user interface, the user may interactively manipulate the default settings to obtain a quality grid within a fraction of time that is usually required for building a grid about the same geometry with a general-purpose grid generation code. Non-Uniform Rational B-Spline formulations are utilized in the algorithm to maintain geometry fidelity while redistributing grid points on the solid surfaces. Bezier curve formulation is used to allow interactive construction of inner boundaries. It is also utilized to allow interactive point distribution. Cascade surfaces are transformed from three-dimensional surfaces of revolution into two-dimensional parametric planes for easy manipulation. Such a transformation allows these manipulated plane grids to be mapped to surfaces of revolution by any generatrix definition. A sophisticated visualization module is developed to al-low visualization for both grid and flow solution, steady or unsteady. A network module is built to allow data transferring in the heterogeneous environment. A flow module is integrated into this system, using an existing turbomachinery flow code. A simulation module is developed to combine the network, flow, and visualization module to achieve near real-time flow simulation about turbomachinery geometries. A toolbox module is developed to support the overall task. A batch version of the grid generation module is developed to allow portability and has been extended to allow dynamic grid generation for pitch changing turbomachinery configurations. Various applications with different characteristics are presented to demonstrate the success of this system.

Shih, Ming-Hsin↗
NASA NTRSDate not supplied

Lunar Water Extraction Via Lunar Auger Dryer ISRU (LADI)

In 2009, the Lunar Reconnaissance Orbiter (LRO) and Lunar Crater Observation and Sensing Satellite (LCROSS) provided definitive proof of water in the Lunar’s southern permanently shadowed region (PSR)1. Both the 2020 NASA Technology Taxonomy2 and the Lunar Surface Innovation Initiative (LSII) team identified capability gaps in icy regolith transfer and reactor processing in Permanently Shadowed Region (PSR) environmental conditions. A water processing plant operating inside the PSR can continuously process water (and volatiles) for both breathable air and propellant. NASA’s Johnson Space Center (JSC) fabricated a unique breadboard test stand for validating the feasibility of this concept,developed a lunar auger dryer sizing tool, and identified the physical flow and thermal models required to design an Engineering Development Unit (EDU) for environmental testing. Over the two years this project was active, models were developed including: a) a sizing tool that scales the auger parameters based on production targets, b) accurate thermal models in both lab and lunar conditions, c) a Discreet Element Model (DEM) for observing particle effects in low-g, and d) incorporated these subsystem models into a MATLAB parametric Modular Open System Approach (MOSA) system model. These models were validated by testing the auger in both a mechanical and thermal system configuration. The results of this testing increased the Technology Readiness Level (TRL) of the Lunar Auger Dryer ISRU (LADI) system to TRL 4.

Lunar↗
NASA NTRS2010

Excavator Design Validation

The Excavator Design Validation tool verifies excavator designs by automatically generating control systems and modeling their performance in an accurate simulation of their expected environment. Part of this software design includes interfacing with human operations that can be included in simulation-based studies and validation. This is essential for assessing productivity, versatility, and reliability. This software combines automatic control system generation from CAD (computer-aided design) models, rapid validation of complex mechanism designs, and detailed models of the environment including soil, dust, temperature, remote supervision, and communication latency to create a system of high value. Unique algorithms have been created for controlling and simulating complex robotic mechanisms automatically from just a CAD description. These algorithms are implemented as a commercial cross-platform C++ software toolkit that is configurable using the Extensible Markup Language (XML). The algorithms work with virtually any mobile robotic mechanisms using module descriptions that adhere to the XML standard. In addition, high-fidelity, real-time physics-based simulation algorithms have also been developed that include models of internal forces and the forces produced when a mechanism interacts with the outside world. This capability is combined with an innovative organization for simulation algorithms, new regolith simulation methods, and a unique control and study architecture to make powerful tools with the potential to transform the way NASA verifies and compares excavator designs. Energid's Actin software has been leveraged for this design validation. The architecture includes parametric and Monte Carlo studies tailored for validation of excavator designs and their control by remote human operators. It also includes the ability to interface with third-party software and human-input devices. Two types of simulation models have been adapted: high-fidelity discrete element models and fast analytical models. By using the first to establish parameters for the second, a system has been created that can be executed in real time, or faster than real time, on a desktop PC. This allows Monte Carlo simulations to be performed on a computer platform available to all researchers, and it allows human interaction to be included in a real-time simulation process. Metrics on excavator performance are established that work with the simulation architecture. Both static and dynamic metrics are included.

Pholsiri, Chalongrath↗
NASA NTRS2014

Analyzing the Impacts of Natural Environments on Launch and Landing Availability for NASA's Exploration Systems Development Programs

The National Aeronautics and Space Administration (NASA) is developing new capabilities for human and scientific exploration beyond Earth orbit. Natural environments information is an important asset for NASA's development of the next generation space transportation system as part of the Exploration Systems Development (ESD) Programs, which includes the Space Launch System (SLS) and Multi-Purpose Crew Vehicle (MPCV) Programs. Natural terrestrial environment conditions - such as wind, lightning and sea states - can affect vehicle safety and performance during multiple mission phases ranging from pre-launch ground processing to landing and recovery operations, including all potential abort scenarios. Space vehicles are particularly sensitive to these environments during the launch/ascent and the entry/landing phases of mission operations. The Marshall Space Flight Center (MSFC) Natural Environments Branch provides engineering design support for NASA space vehicle projects and programs by providing design engineers and mission planners with natural environments definitions as well as performing custom analyses to help characterize the impacts the natural environment may have on vehicle performance. One such analysis involves assessing the impact of natural environments to operational availability. Climatological time series of operational surface weather observations are used to calculate probabilities of meeting/exceeding various sets of hypothetical vehicle-specific parametric constraint thresholds. Outputs are tabulated by month and hour of day to show both seasonal and diurnal variation. This paper will discuss how climate analyses are performed by the MSFC Natural Environments Branch to support the ESD Launch Availability (LA) Technical Performance Measure (TPM), the SLS Launch Availability due to Natural Environments TPM, and several MPCV (Orion) launch and landing availability analyses - including the 2014 Orion Exploration Flight Test 1 (EFT-1) mission.

Altino, Karen M.↗
NASA NTRS1994

Optimizations on Designing High-Resolution Finite-Difference Schemes

We describe a general optimization procedure for both maximizing the resolution characteristics of existing finite differencing schemes as well as designing finite difference schemes that will meet the error tolerance requirements of numerical solutions. The procedure is based on an optimization process. This is a generalization of the compact scheme introduced by Lele in which the resolution is improved for single, one-dimensional spatial derivative, whereas in the present approach the complete scheme, after spatial and temporal discretizations, is optimized on a range of parameters of the scheme and the governing equations. The approach is to linearize and Fourier analyze the discretized equations to check the resolving power of the scheme for various wave number ranges in the solution and optimize the resolution to satisfy the requirements of the problem. This represents a constrained nonlinear optimization problem which can be solved to obtain the nodal weights of discretization. An objective function is defined in the parametric space of wave numbers, Courant number, Mach number and other quantities of interest. Typical criterion for defining the objective function include the maximization of the resolution of high wave numbers for acoustic and electromagnetic wave propagations and turbulence calculations. The procedure is being tested on off-design conditions of non-uniform mesh, non-periodic boundary conditions, and non-constant wave speeds for scalar and system of equations. This includes the solution of wave equations and Euler equations using a conventional scheme with and without optimization and the design of an optimum scheme for the specified error tolerance.

Liu, Yen↗
NASA NTRS1999

Brain segmentation and the generation of cortical surfaces

This paper describes methods for white matter segmentation in brain images and the generation of cortical surfaces from the segmentations. We have developed a system that allows a user to start with a brain volume, obtained by modalities such as MRI or cryosection, and constructs a complete digital representation of the cortical surface. The methodology consists of three basic components: local parametric modeling and Bayesian segmentation; surface generation and local quadratic coordinate fitting; and surface editing. Segmentations are computed by parametrically fitting known density functions to the histogram of the image using the expectation maximization algorithm [DLR77]. The parametric fits are obtained locally rather than globally over the whole volume to overcome local variations in gray levels. To represent the boundary of the gray and white matter we use triangulated meshes generated using isosurface generation algorithms [GH95]. A complete system of local parametric quadratic charts [JWM+95] is superimposed on the triangulated graph to facilitate smoothing and geodesic curve tracking. Algorithms for surface editing include extraction of the largest closed surface. Results for several macaque brains are presented comparing automated and hand surface generation. Copyright 1999 Academic Press.

NASA Discipline Neuroscience↗