Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MINIMIZATION”

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,765 records · Page 98

Novel Integrated Radial and Axial Magnetic Bearing

Typically, fully active magnetically suspended systems require one axial and two radial magnetic bearings. Combining radial and axial functions into a single device allows for more compact and elegant packaging. Furthermore, in the case of high-speed devices such as energy storage flywheels, it is beneficial to minimize shaft length to keep rotor mode frequencies as high as possible. Attempts have been made to combine radial and axial functionality, but with certain drawbacks. One approach requires magnetic control flux to flow through a bias magnet reducing control effectiveness, thus resulting in increased resistive losses. This approach also requires axial force producing magnetic flux to flow in a direction into the rotor laminate that is undesirable for minimizing eddy-current losses resulting in rotational losses. Another approach applies a conical rotor shape to what otherwise would be a radial heteropolar magnetic bearing configuration. However, positional non-linear effects are introduced with this scheme and the same windings are used for bias, radial, and axial control adding complexity to the controller and electronics. For this approach, the amount of axial capability must be limited. It would be desirable for an integrated radial and axial magnetic bearing to have the following characteristics, separate inputs for radial and axial control for electronics and control simplicity, all magnetic control fluxes should only flow through their respective air gaps and should not flow through any bias magnets for minimal resistive losses, be of a homopolar design to minimize rotational losses, position related non-linear effects should be minimized, and dependent upon the design parameters, be able to achieve any radial/axial force or power ratio as desired. The integrated radial and axial magnetic bearing described in this paper exhibits all these characteristics. Magnetic circuit design, design equations, and analysis results will be presented.

Blumenstock, Kenneth A.↗

Novel Integration Radial and Axial Magnetic Bearing

Typically, fully active magnetically suspended systems require one axial and two radial magnetic bearings. Combining radial and axial functions into a single device allows for more compact and elegant packaging. Furthermore, in the case of high-speed devices such as energy storage flywheels, it is beneficial to minimize shaft length to keep rotor mode frequencies as high as possible. Attempts have been made to combine radial and axial functionality, but with certain drawbacks. One approach requires magnetic control flux to flow through a bias magnet reducing control effectiveness, thus resulting in increased resistive losses. This approach also requires axial force producing magnetic flux to flow in a direction into the rotor laminate that is undesirable for minimizing eddy-current losses resulting in rotational losses. Another approach applies a conical rotor shape to what otherwise would be a radial heteropolar magnetic bearing configuration. However, positional non-linear effects are introduced with this scheme and the same windings are used for bias, radial, and axial control adding complexity to the controller and electronics. For this approach, the amount of axial capability must be limited. It would be desirable for an integrated radial and axial magnetic bearing to have the following characteristics; separate inputs for radial and axial control for electronics and control simplicity, all magnetic control fluxes should only flow through their respective air gaps and should not flow through any bias magnets for minimal resistive losses, be of a homopolar design to minimize rotational losses, position related non-linear effects should be minimized, and dependent upon the design parameters, be able to achieve any radial/axial force or power ratio as desired. The integrated radial and axial magnetic bearing described in this paper exhibits all these characteristics. Magnetic circuit design, design equations, and magnetic field modeling results will be presented.

Blumenstock, Kenneth↗

Noise Considerations in the Design and Operation of V/STOL Aircraft

Available propulsion-system noise data have been applied to the problems of design and operation of V/STOL aircraft. In particular, considerations have been given to minimizing adverse community reaction for operations between airports and to minimizing detection due to noise for special missions. For minimizing adverse community reaction, configurations incorporating low-blade-loading rotors, low-tip-speed propellers, or turbofan-type engines are judged to be most satisfactory. For minimizing detection, consideration must be given to minimizing the noise of the generating airplane, having maximum background noise in the vicinity of the observer, and operating the aircraft at minimum altitude.

Maglieri, Domenic J.↗

Optimized suspension culture: the rotating-wall vessel

Suspension culture remains a popular modality, which manipulates mechanical culture conditions to maintain the specialized features of cultured cells. The rotating-wall vessel is a suspension culture vessel optimized to produce laminar flow and minimize the mechanical stresses on cell aggregates in culture. This review summarizes the engineering principles, which allow optimal suspension culture conditions to be established, and the boundary conditions, which limit this process. We suggest that to minimize mechanical damage and optimize differentiation of cultured cells, suspension culture should be performed in a solid-body rotation Couette-flow, zero-headspace culture vessel such as the rotating-wall vessel. This provides fluid dynamic operating principles characterized by 1) solid body rotation about a horizontal axis, characterized by colocalization of cells and aggregates of different sedimentation rates, optimally reduced fluid shear and turbulence, and three-dimensional spatial freedom; and 2) oxygenation by diffusion. Optimization of suspension culture is achieved by applying three tradeoffs. First, terminal velocity should be minimized by choosing microcarrier beads and culture media as close in density as possible. Next, rotation in the rotating-wall vessel induces both Coriolis and centrifugal forces, directly dependent on terminal velocity and minimized as terminal velocity is minimized. Last, mass transport of nutrients to a cell in suspension culture depends on both terminal velocity and diffusion of nutrients. In the transduction of mechanical culture conditions into cellular effects, several lines of evidence support a role for multiple molecular mechanisms. These include effects of shear stress, changes in cell cycle and cell death pathways, and upstream regulation of secondary messengers such as protein kinase C. The discipline of suspension culture needs a systematic analysis of the relationship between mechanical culture conditions and biological effects, emphasizing cellular processes important for the industrial production of biological pharmaceuticals and devices.

Non-NASA Center↗

Attitude-Independent Magnetometer Calibration for Spin-Stabilized Spacecraft

The paper describes a three-step estimator to calibrate a Three-Axis Magnetometer (TAM) using TAM and slit Sun or star sensor measurements. In the first step, the Calibration Utility forms a loss function from the residuals of the magnitude of the geomagnetic field. This loss function is minimized with respect to biases, scale factors, and nonorthogonality corrections. The second step minimizes residuals of the projection of the geomagnetic field onto the spin axis under the assumption that spacecraft nutation has been suppressed by a nutation damper. Minimization is done with respect to various directions of the body spin axis in the TAM frame. The direction of the spin axis in the inertial coordinate system required for the residual computation is assumed to be unchanged with time. It is either determined independently using other sensors or included in the estimation parameters. In both cases all estimation parameters can be found using simple analytical formulas derived in the paper. The last step is to minimize a third loss function formed by residuals of the dot product between the geomagnetic field and Sun or star vector with respect to the misalignment angle about the body spin axis. The method is illustrated by calibrating TAM for the Fast Auroral Snapshot Explorer (FAST) using in-flight TAM and Sun sensor data. The estimated parameters include magnetic biases, scale factors, and misalignment angles of the spin axis in the TAM frame. Estimation of the misalignment angle about the spin axis was inconclusive since (at least for the selected time interval) the Sun vector was about 15 degrees from the direction of the spin axis; as a result residuals of the dot product between the geomagnetic field and Sun vectors were to a large extent minimized as a by-product of the second step.

Natanson, Gregory↗

Efficient Model-Based Diagnosis Engine

An efficient diagnosis engine - a combination of mathematical models and algorithms - has been developed for identifying faulty components in a possibly complex engineering system. This model-based diagnosis engine embodies a twofold approach to reducing, relative to prior model-based diagnosis engines, the amount of computation needed to perform a thorough, accurate diagnosis. The first part of the approach involves a reconstruction of the general diagnostic engine to reduce the complexity of the mathematical-model calculations and of the software needed to perform them. The second part of the approach involves algorithms for computing a minimal diagnosis (the term "minimal diagnosis" is defined below). A somewhat lengthy background discussion is prerequisite to a meaningful summary of the innovative aspects of the present efficient model-based diagnosis engine. In model-based diagnosis, the function of each component and the relationships among all the components of the engineering system to be diagnosed are represented as a logical system denoted the system description (SD). Hence, the expected normal behavior of the engineering system is the set of logical consequences of the SD. Faulty components lead to inconsistencies between the observed behaviors of the system and the SD (see figure). Diagnosis - the task of finding faulty components - is reduced to finding those components, the abnormalities of which could explain all the inconsistencies. The solution of the diagnosis problem should be a minimal diagnosis, which is a minimal set of faulty components. A minimal diagnosis stands in contradistinction to the trivial solution, in which all components are deemed to be faulty, and which, therefore, always explains all inconsistencies.

Fijany, Amir↗

Low-Friction, High-Stiffness Joint for Uniaxial Load Cell

A universal-joint assembly has been devised for transferring axial tension or compression to a load cell. To maximize measurement accuracy, the assembly is required to minimize any moments and non-axial forces on the load cell and to exhibit little or no hysteresis. The requirement to minimize hysteresis translates to a requirement to maximize axial stiffness (including minimizing backlash) and a simultaneous requirement to minimize friction. In practice, these are competing requirements, encountered repeatedly in efforts to design universal joints. Often, universal-joint designs represent compromises between these requirements. The improved universal-joint assembly contains two universal joints, each containing two adjustable pairs of angular-contact ball bearings. One might be tempted to ask why one could not use simple ball-and-socket joints rather than something as complex as universal joints containing adjustable pairs of angularcontact ball bearings. The answer is that ball-and-socket joints do not offer sufficient latitude to trade stiffness versus friction: the inevitable result of an attempt to make such a trade in a ball-and-socket joint is either too much backlash or too much friction. The universal joints are located at opposite ends of an axial subassembly that contains the load cell. The axial subassembly includes an axial shaft, an axial housing, and a fifth adjustable pair of angular-contact ball bearings that allows rotation of the axial housing relative to the shaft. The preload on each pair of angular-contact ball bearings can be adjusted to obtain the required stiffness with minimal friction, tailored for a specific application. The universal joint at each end affords two degrees of freedom, allowing only axial force to reach the load cell regardless of application of moments and non-axial forces. The rotational joint on the axial subassembly affords a fifth degree of freedom, preventing application of a torsion load to the load cell.

Lewis, James L.↗

Miniature L-Band Radar Transceiver

A miniature L-band transceiver that operates at a carrier frequency of 1.25 GHz has been developed as part of a generic radar electronics module (REM) that would constitute one unit in an array of many identical units in a very-large-aperture phased-array antenna. NASA and the Department of Defense are considering the deployment of such antennas in outer space; the underlying principles of operation, and some of those of design, also are applicable on Earth. The large dimensions of the antennas make it advantageous to distribute radio-frequency electronic circuitry into elements of the arrays. The design of the REM is intended to implement the distribution. The design also reflects a requirement to minimize the size and weight of the circuitry in order to minimize the weight of any such antenna. Other requirements include making the transceiver robust and radiation-hard and minimizing power demand. Figure 1 depicts the functional blocks of the REM, including the L-band transceiver. The key functions of the REM include signal generation, frequency translation, amplification, detection, handling of data, and radar control and timing. An arbitrary-waveform generator that includes logic circuitry and a digital-to-analog converter (DAC) generates a linear-frequency-modulation chirp waveform. A frequency synthesizer produces local-oscillator signals used for frequency conversion and clock signals for the arbitrary-waveform generator, for a digitizer [that is, an analog-to-digital converter (ADC)], and for a control and timing unit. Digital functions include command, timing, telemetry, filtering, and high-rate framing and serialization of data for a high-speed scientific-data interface. The aforementioned digital implementation of filtering is a key feature of the REM architecture. Digital filters, in contradistinction to analog ones, provide consistent and temperature-independent performance, which is particularly important when REMs are distributed throughout a large array. Digital filtering also enables selection among multiple filter parameters as required for different radar operating modes. After digital filtering, data are decimated appropriately in order to minimize the data rate out of an antenna panel. The L-band transceiver (see Figure 2) includes a radio-frequency (RF)-to-baseband down-converter chain and an intermediate- frequency (IF)-to-RF up-converter chain. Transmit/receive (T/R) switches enable the use of a single feed to the antenna for both transmission and reception. The T/R switches also afford a built-in test capability by enabling injection of a calibration signal into the receiver chain. In order of decreasing priority, components of the transceiver were selected according to requirements of radiation hardness, then compactness, then low power. All of the RF components are radiation-hard. The noise figure (NF) was optimized to the extent that (1) a low-noise amplifier (LNA) (characterized by NF < 2 dB) was selected but (2) the receiver front-end T/R switches were selected for a high degree of isolation and acceptably low loss, regardless of the requirement to minimize noise.

McWatters, Dalia↗

Resource Balancing Control Allocation

Next generation aircraft with a large number of actuators will require advanced control allocation methods to compute the actuator commands needed to follow desired trajectories while respecting system constraints. Previously, algorithms were proposed to minimize the l1 or l2 norms of the tracking error and of the control effort. The paper discusses the alternative choice of using the l1 norm for minimization of the tracking error and a normalized l(infinity) norm, or sup norm, for minimization of the control effort. The algorithm computes the norm of the actuator deflections scaled by the actuator limits. Minimization of the control effort then translates into the minimization of the maximum actuator deflection as a percentage of its range of motion. The paper shows how the problem can be solved effectively by converting it into a linear program and solving it using a simplex algorithm. Properties of the algorithm are investigated through examples. In particular, the min-max criterion results in a type of resource balancing, where the resources are the control surfaces and the algorithm balances these resources to achieve the desired command. A study of the sensitivity of the algorithms to the data is presented, which shows that the normalized l(infinity) algorithm has the lowest sensitivity, although high sensitivities are observed whenever the limits of performance are reached.

Frost, Susan A.↗

Passive Water Assurance Delivery System

This report provides an analysis and evaluation of proposed design improvements to the internal water delivery system of NASA’s passive vegetable growing system. The VEGGIE system is designed to provide fresh vegetables and psychological benefits for astronauts aboard the International Space Station. Preliminary testing conducted by NASA for the VEGGIE system revealed flaws in the water delivery system. In these tests, the interface between the plants and the water reservoir failed to provide passive water delivery through the system both in microgravity experiments and experiments conducted on Earth. Our team’s goal was to fix this interface such that water can be passively delivered from the water reservoir to the plant rooting pillows to grow vegetative crops from seed to harvest. The system requirements outlined by NASA include: minimal total mass, on-demand passive water delivery, same dimensions as existing fixtures, minimal pressure on system to prevent leaking, zero mold growth, minimal swelling or clogging for non-flammable capillary materials, pressure stabilization between the water reservoir and plant pillow bags, even dispersal of water during initial priming, minimal bubble obstruction of capillary interface, maximum gas availability to plant roots, growth in a mixed artificial media (50:50 Arcillite:Fafard #2), and to avoid overwatering and drought conditions for plants. The major proposed design improvements are to use a single-interface capillary cord design to directly connect the water reservoir to the plant-rooting pillows, and to alter the water reservoir to model a propellant management device (PMD) in order to ensure consistent and long term watering for the VEGGIE system. The plant rooting pillow required minimal changes outside of replacing the capillary mat on the bottom of the pillow with O-ring insertion points the single-interface capillary system. - Research was conducted to determine which materials are able to uptake water through capillary action to grow Outredgeous Lettuce plants from seed to harvest. The primary requirements to be met by the team's design and evaluation of the capillary interface were: continuous passive watering for 90 days using non-flammable capillary materials and a peak water delivery rate of 30 mL/hr./0.15 m 2 . The results of the team's experimentation showed that the capillary material Nomex displayed the highest capillary water delivery potential with a maximum flow rate of 3.6 mL/hr. The Nomex capillary systems were the only capillary material to consistently grow Outredgeous Lettuce plants from seed to harvest, and displayed the highest average flow rate for multiple experiment sets. The Nomex material was previously incorporated into the VEGGIE system using a matted version of the capillary material, though the team recommends using a cord configuration in the single-interface capillary system for greater system stability. Nomex has proven to be a promising material, as it has passed both health and fire standards for use aboard the International Space Station. PMDs are made of materials that utilize surface tension and adhesive forces to overcome adverse accelerations to improve stability and ensure fluid delivery. PMDs are typically used in fuel tanks to ensure fuel delivery. The team recommends using a sponge PMD in order to mitigate bubble obstruction, decrease system weight, and ensure reliable water delivery to the capillary interface. The PMD water reservoir requires a rigid water reservoir and a vent tube for pressure stabilization. Since PMD’s cannot be tested in 1-G (Earth conditions), further theoretical modeling and testing is required for the proposed water reservoir design. The intent of this proposed system is to passively water plants in microgravity, though the technology is not limited to microgravity applications. The testing at The Ohio State University has proven that the design is highly effective on Earth, demonstrating that it could serve as a simple water delivery system in home and office applications. This would make vegetative crops more accessible in all indoor applications, thereby improving indoor air quality and occupant comfort. Additionally, this technology has great potential to be utilized in greenhouse plant production, cutting back on more sophisticated watering system energy and time requirements.

Alexandria Jensen↗

Development of 3-D Geological Model of Tuscarora Sandstone for Feasibility of Deep Direct-Use Geothermal at West Virginia University’s Main Campus

The Morgantown campus of West Virginia University (WVU) is uniquely positioned to host the first geothermal deep direct-use district heating system in the eastern United States. While much of the eastern United States is not blessed with extremely high heat flow and elevated temperatures, the northeastern part of West Virginia is unique in having a basin that is expected to support the achievable flowrate of geofluid through target formations, with sufficient temperatures at depth. These two factors were identified in the 2006 MIT Future of Geothermal Energy Report to be the two most critical factors in minimizing cost of geothermal energy. Our overall project objective is to determine the feasibility of designing a Geothermal District Heating and Cooling (GDHC) system for the West Virginia University campus utilizing Geothermal Deep Direct-Use (DDU) by 1) minimization of the uncertainty and risk associated with developing the geothermal resource for use on campus at WVU and 2) completion of an optimized design for the geothermal system, minimizing the delivered Levelized Cost of Heat (LCOH). Our first goal, to minimize the risk of project development, will be achieved by decreasing the uncertainty in both the subsurface geothermal system as well the surface distribution system. The subsurface uncertainty is dominated by the uncertainty in the projections of geofluid flowrate in the target formation, the Tuscarora Sandstone. In this paper, three cores from the heterogeneous reservoir, available through West Virginia Geologic and Economic Survey, are analyzed by performing core analysis using CT scanning and permeability measurements via minipermeameter. Additional geological data are collected through cores, published literature, seismic data, and nearby, existing wells to estimate thickness, fracture network configuration and geothermal gradient to minimize the uncertainty of well deliverability. Using these estimated reservoir properties; a 3D conceptual model for the proposed geothermal site is developed.

15 GEOTHERMAL ENERGY↗

In-situ Detection of Glovebox Glove Degradation Prior to Glove Failure; Effects on Safety, Waste, and Operational Costs - 20402

Glovebox glove pressure decay leak testing and analysis of the data can greatly improve glovebox safety, minimize Transuranic (TRU) glove waste, and minimize operational cost. It is well known, the weakest point of containment on a glovebox is the glove. Mitigating unplanned openings in gloves is critical to minimizing operational and safety costs in glovebox operations. Mitigating unplanned glove openings due to glove failures will be discussed in this paper. The lack of an engineered solution to determine a gloves operational life in nuclear operations, has lead to use of theoretical and manual subjective methods to determine the life and safety of a glovebox glove. Whether the glove is inspected visually prior to use, or gloves are changed regularly to theoretically avoid failure, both the operators PPE and the room are being exposed to contamination, or un-necessary amounts of TRU glove waste are being generated. Millions of dollars yearly in operational costs are spent on glove failures due to loss of production, incident analysis, regulatory audits, clean up, and paperwork. In a nuclear application, a glove failure could contaminate an operator and laboratory, shut down the operation for weeks, potentially create a regulatory audit, and potentially create a media frenzy; the damages to the organization can be astronomical. Elimination of glove failures is possible by tracking glove material degradation and setting limits to allow glove change when it becomes necessary, prior to glove failure. Analytical pressure decay leak testing allows the detection of material degradation, and in turn the ability to limit operational glove failures. A German based company called MK Versuchsanlagen, e.K., has developed an advanced Glove Integrity Testing System. The system, with its use of proprietary RFID and special software technology, is unique in its capabilities. The system is capable to perform highly accurate regular glove testing on any size glovebox line in minimal time. The system records and can analyze a tremendous amount of system data that in turn can determine, for example, a safe glove change time prior to a glove material failure. This paper will show how an advanced Glove Integrity Testing System can detect glove material degradation and help determine a safe change point prior to glove failure. Pressure decay curves of new gloves and recorded effects of accelerated aging on the glove over time will be shown. The use of this technology in nuclear facilities can greatly improve glovebox operational safety, minimize glovebox glove TRU waste, and save organizations tremendous costs in downtime, clean up, documentation, and potential regulatory review. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Cardiac C-arm computed tomography using a 3D + time ROI reconstruction method with spatial and temporal regularization

Purpose: Reconstruction of the beating heart in 3D + time in the catheter laboratory using only the available C-arm system would improve diagnosis, guidance, device sizing, and outcome control for intracardiac interventions, e.g., electrophysiology, valvular disease treatment, structural or congenital heart disease. To obtain such a reconstruction, the patient's electrocardiogram (ECG) must be recorded during the acquisition and used in the reconstruction. In this paper, the authors present a 4D reconstruction method aiming to reconstruct the heart from a single sweep 10 s acquisition. Methods: The authors introduce the 4D RecOnstructiOn using Spatial and TEmporal Regularization (short 4D ROOSTER) method, which reconstructs all cardiac phases at once, as a 3D + time volume. The algorithm alternates between a reconstruction step based on conjugate gradient and four regularization steps: enforcing positivity, averaging along time outside a motion mask that contains the heart and vessels, 3D spatial total variation minimization, and 1D temporal total variation minimization. Results: 4D ROOSTER recovers the different temporal representations of a moving Shepp and Logan phantom, and outperforms both ECG-gated simultaneous algebraic reconstruction technique and prior image constrained compressed sensing on a clinical case. It generates 3D + time reconstructions with sharp edges which can be used, for example, to estimate the patient's left ventricular ejection fraction. Conclusions: 4D ROOSTER can be applied for human cardiac C-arm CT, and potentially in other dynamic tomography areas. It can easily be adapted to other problems as regularization is decoupled from projection and back projection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Theoretically informed Monte Carlo simulation of liquid crystals by sampling of alignment-tensor fields

A theoretically informed coarse-grained Monte Carlo method is proposed for studying liquid crystals. The free energy functional of the system is described in the framework of the Landau-de Gennes formalism. The alignment field and its gradients are approximated by finite differences, and the free energy is minimized through a stochastic sampling technique. The validity of the proposed method is established by comparing the results of the proposed approach to those of traditional free energy minimization techniques. Its usefulness is illustrated in the context of three systems, namely, a nematic liquid crystal confined in a slit channel, a nematic liquid crystal droplet, and a chiral liquid crystal in the bulk. It is found that for systems that exhibit multiple metastable morphologies, the proposed Monte Carlo method is generally able to identify lower free energy states that are often missed by traditional approaches. Importantly, the Monte Carlo method identifies such states from random initial configurations, thereby obviating the need for educated initial guesses that can be difficult to formulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optimized mirror shape tuning using beam weightings based on distance, angle of incidence, reflectivity, and power

For glancing-incidence optical systems, such as short-wavelength optics used for nano-focusing, incorporating physical factors in the calculations used for shape optimization can improve performance. Wavefront metrology, including the measurement of a mirror's shape or slope, is routinely used as input for mirror figure optimization on mirrors that can be bent, actuated, positioned, or aligned. Modeling shows that when the incident power distribution, distance from focus, angle of incidence, and the spatially varying reflectivity are included in the optimization, higher Strehl ratios can be achieved. Following the works of Maréchal and Mahajan, optimization of the Strehl ratio (for peak intensity with a coherently illuminated system) occurs when the expectation value of the phase error's variance is minimized. We describe an optimization procedure based on regression analysis that incorporates these physical parameters. This approach is suitable for coherently illuminated systems of nearly diffraction-limited quality. Mathematically, this work is an enhancement of the methods commonly applied for ex situ alignment based on uniform weighting of all points on the surface (or a sub-region of the surface). It follows a similar approach to the optimization of apodized and non-uniformly illuminated optical systems. Significantly, it reaches a different conclusion than a more recent approach based on minimization of focal plane ray errors.

47 OTHER INSTRUMENTATION↗

Simultaneous motion estimation and image reconstruction (SMEIR) for 4D cone-beam CT

Purpose: Image reconstruction and motion model estimation in four-dimensional cone-beam CT (4D-CBCT) are conventionally handled as two sequential steps. Due to the limited number of projections at each phase, the image quality of 4D-CBCT is degraded by view aliasing artifacts, and the accuracy of subsequent motion modeling is decreased by the inferior 4D-CBCT. The objective of this work is to enhance both the image quality of 4D-CBCT and the accuracy of motion model estimation with a novel strategy enabling simultaneous motion estimation and image reconstruction (SMEIR).Methods: The proposed SMEIR algorithm consists of two alternating steps: (1) model-based iterative image reconstruction to obtain a motion-compensated primary CBCT (m-pCBCT) and (2) motion model estimation to obtain an optimal set of deformation vector fields (DVFs) between the m-pCBCT and other 4D-CBCT phases. The motion-compensated image reconstruction is based on the simultaneous algebraic reconstruction technique (SART) coupled with total variation minimization. During the forward- and backprojection of SART, measured projections from an entire set of 4D-CBCT are used for reconstruction of the m-pCBCT by utilizing the updated DVF. The DVF is estimated by matching the forward projection of the deformed m-pCBCT and measured projections of other phases of 4D-CBCT. The performance of the SMEIR algorithm is quantitatively evaluated on a 4D NCAT phantom. The quality of reconstructed 4D images and the accuracy of tumor motion trajectory are assessed by comparing with those resulting from conventional sequential 4D-CBCT reconstructions (FDK and total variation minimization) and motion estimation (demons algorithm). The performance of the SMEIR algorithm is further evaluated by reconstructing a lung cancer patient 4D-CBCT.Results: Image quality of 4D-CBCT is greatly improved by the SMEIR algorithm in both phantom and patient studies. When all projections are used to reconstruct a 3D-CBCT by FDK, motion-blurring artifacts are present, leading to a 24.4% relative reconstruction error in the NACT phantom. View aliasing artifacts are present in 4D-CBCT reconstructed by FDK from 20 projections, with a relative error of 32.1%. When total variation minimization is used to reconstruct 4D-CBCT, the relative error is 18.9%. Image quality of 4D-CBCT is substantially improved by using the SMEIR algorithm and relative error is reduced to 7.6%. The maximum error (MaxE) of tumor motion determined from the DVF obtained by demons registration on a FDK-reconstructed 4D-CBCT is 3.0, 2.3, and 7.1 mm along left–right (L-R), anterior–posterior (A-P), and superior–inferior (S-I) directions, respectively. From the DVF obtained by demons registration on 4D-CBCT reconstructed by total variation minimization, the MaxE of tumor motion is reduced to 1.5, 0.5, and 5.5 mm along L-R, A-P, and S-I directions. From the DVF estimated by SMEIR algorithm, the MaxE of tumor motion is further reduced to 0.8, 0.4, and 1.5 mm along L-R, A-P, and S-I directions, respectively.Conclusions: The proposed SMEIR algorithm is able to estimate a motion model and reconstruct motion-compensated 4D-CBCT. The SMEIR algorithm improves image reconstruction accuracy of 4D-CBCT and tumor motion trajectory estimation accuracy as compared to conventional sequential 4D-CBCT reconstruction and motion estimation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Electrode Edge Effects and the Failure Mechanism of Lithium-Metal Batteries

Lithium (Li) metal with very high specific capacity makes it one of the ideal anode for high energy batteries. However, Li dendrite growth and the formation of isolated (or “dead”) Li during repeated Li plating/stripping processes lead to its low Coulombic efficiency (CE). In this work, we discovered for the first time that part of Li CE loss can be attributed to the edge effect of the copper (Cu) substrate electrode in a conventional Li||Cu cell configuration. The “dead” Li formed on the edge of Cu substrate was systematically investigated through scanning electron microscopy, energy dispersive X-ray spectroscopy, and two-dimensional X-ray photoelectron spectroscopy. In order to minimize the Li loss at the edge of the Cu exposed to pressure-free space, modified Li||Cu cell configuration with a Cu electrode smaller than Li metal is preferred. It was clearly demonstrated that using an electrode configuration with minimal open space or “pressure-free space” across electrodes can reduce “dead” Li accumulated during cycling and increase Li CE. This phenomena was also verified in Li metal batteries (Li||LiNi1/3Mn1/3Co1/3O2) and should be considered in the design of practical Li metal batteries.

Lee, Hongkyung↗

Parallel FETI‐DP algorithm for efficient simulation of large‐scale EM problems

Summary An efficient parallelization of the dual‐primal finite‐element tearing and interconnecting (FETI‐DP) algorithm is presented for large‐scale electromagnetic simulations. As a nonoverlapping domain decomposition method, the FETI‐DP algorithm formulates a global interface problem, whose iterative solution is accelerated with a solution of a global corner problem. To achieve a good load balance for parallel computation, the original computational domain is decomposed into subdomains with similar sizes and shapes. The subdomains are then distributed to processors based on their close proximity to minimize inter‐processor communication. The parallel generalized minimal residual method, enhanced with the iterative classical Gram‐Schmidt orthogonalization scheme to reduce global communication, is adopted to solve the global interface problem with a fast convergence rate. The global corner‐related coarse problem is solved iteratively with a parallel communication‐avoiding biconjugate gradient stabilized method to minimize global communication, and its convergence is accelerated by a diagonal preconditioner constructed from the coarse system matrix. To alleviate neighboring communication overhead, the non‐blocking communication approach is employed in both generalized minimal residual and communication‐avoiding biconjugate gradient stabilized iterative solutions. Three numerical examples are presented to demonstrate the accuracy, scalability, and capability of the proposed parallel FETI‐DP algorithm for electromagnetic modeling of general objects and antenna arrays. Copyright © 2016 John Wiley & Sons, Ltd.

Zhang, Kedi↗