Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shape from gradient”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗
DOE OSTI2019

Observing crystal nucleation in four dimensions using atomic electron tomography

Nucleation plays a critical role in many physical and biological phenomena that range from crystallization, melting and evaporation to the formation of clouds and the initiation of neurodegenerative diseases. However, nucleation is a challenging process to study in experiments, especially in the early stage when several atoms or molecules start to form a new phase from a parent phase. A number of experimental and computational methods have been used to investigate nucleation processes, but experimental determination of the three-dimensional atomic structure and the dynamics of early stage nuclei has been unachievable. Here we use atomic electron tomography to study early stage nucleation in four dimensions (4D: that is, including time) at atomic resolution. Using FePt nanoparticles as a model system, we find that early stage nuclei are irregularly shaped, each has a core of one to a few atoms with the maximum order parameter, and the order parameter gradient points from the core to the boundary of the nucleus. We capture the structure and dynamics of the same nuclei undergoing growth, fluctuation, dissolution, merging and/or division, which are regulated by the order parameter distribution and its gradient. These experimental observations are corroborated by molecular dynamics simulations of heterogeneous and homogeneous nucleation in liquid–solid phase transitions of Pt. Our experimental and molecular dynamics results differ from classical nucleation theory, indicating that a theory beyond this is needed to describe early stage nucleation at the atomic scale. Looking forward, we anticipate that the reported approach will open the door to the study of many fundamental problems in materials science, nanoscience, condensed matter physics and chemistry, such as phase transition, atomic diffusion, grain boundary dynamics, interface motion, defect dynamics and surface reconstruction with 4D atomic resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Single Shot Characterization of High Transformer Ratio Wakefields in Nonlinear Plasma Acceleration

We report that plasma wakefields can enable very high accelerating gradients for frontier high energy particle accelerators, in excess of 10 GeV/m. To overcome limits on single stage acceleration, specially shaped drive beams can be used in both linear and nonlinear plasma wakefield accelerators (PWFA), to increase the transformer ratio, implying that the drive beam deceleration is minimized relative to acceleration obtained in the wake. In this Letter, we report the results of a nonlinear PWFA, high transformer ratio experiment using high-charge, longitudinally asymmetric drive beams in a plasma cell. An emittance exchange process is used to generate variable drive current profiles, in conjunction with a long (multiple plasma wavelength) witness beam. The witness beam is energy modulated by the wakefield, yielding a response that contains detailed spectral information in a single-shot measurement. Using these methods, we generate a variety of beam profiles and characterize the wakefields, directly observing transformer ratios up to R=7.8. Lastly, a spectrally based reconstruction technique, validated by 3D particle-in-cell simulations, is introduced to obtain the drive beam current profile from the decelerating wake data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

WE-D-BRA-07: Analysis of ArcCHECK Diode Array Performance for ViewRay Quality Assurance

Purpose: Discrepancies in absolute dose values have been detected between the ViewRay treatment planning system and ArcCHECK readings when performing delivery quality assurance on the ViewRay system with the ArcCHECK-MR diode array (SunNuclear Corporation). In this work, we investigate whether these discrepancies are due to errors in the ViewRay planning and/or delivery system or due to errors in the ArcCHECK’s readings. Methods: Gamma analysis was performed on 19 ViewRay patient plans using the ArcCHECK. Frequency analysis on the dose differences was performed. To investigate whether discrepancies were due to measurement or delivery error, 10 diodes in low-gradient dose regions were chosen to compare with ion chamber measurements in a PMMA phantom with the same size and shape as the ArcCHECK, provided by SunNuclear. The diodes chosen all had significant discrepancies in absolute dose values compared to the ViewRay TPS. Absolute doses to PMMA were compared between the ViewRay TPS calculations, ArcCHECK measurements, and measurements in the PMMA phantom. Results: Three of the 19 patient plans had 3%/3mm gamma passing rates less than 95%, and ten of the 19 plans had 2%/2mm passing rates less than 95%. Frequency analysis implied a non-random error process. Out of the 10 diode locations measured, ion chamber measurements were all within 2.2% error relative to the TPS and had a mean error of 1.2%. ArcCHECK measurements ranged from 4.5% to over 15% error relative to the TPS and had a mean error of 8.0%. Conclusion: The ArcCHECK performs well for quality assurance on the ViewRay under most circumstances. However, under certain conditions the absolute dose readings are significantly higher compared to the planned doses. As the ion chamber measurements consistently agree with the TPS, it can be concluded that the discrepancies are due to ArcCHECK measurement error and not TPS or delivery system error. This work was funded by the Bhudatt Paliwal Professorship and the University of Wisconsin Medical Radiation Research Center.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Assessment of the Viability of Scaled Annular Pellet Fabrication Technologies

The DOE-NE AFC program broadly supports the creation of novel fuel forms including intermetallics, composites, and dopants for high uranium density compounds. While data can be obtained for key performance indicators such as thermal expansion and oxidation resistance in steam, fuel qualification remains a time intensive process. As such, accelerated fuel qualification methods are a recent focus of the program. By first looking at systems with known historical irradiation testing information, like UO 2 , the method of accelerated burnup testing can be validated, followed by use on novel materials with less information available, like silicides and nitrides. The biggest time delay in doing irradiation testing is reaching a “high burnup” fuel, to effectively model how the fuel will behave when it has been fully used. Reaching this point conventionally could take upwards of 30 years. Accelerated burnup would be approaching the threshold of “high burnup”, something on the order or 50GWd/MTU, more quickly, which although achievable in test reactors requires a decade to achieve using current irradiation testing schemes. To achieve this faster rate, the fuel is likely going to be subject to thermal effects commiserate with the increased heat by more fission. In a normal fuel pellet, this increased temperature would operate less efficiently because the thermal diffusivity of the pellet only allows it to transfer heat so quickly, thus causing a high centerline temperature and noticeable stress along the diameter of the pellet. Solutions could include smaller pellets, which would have a greater surface area giving off heat relative to their volume, or “shaped” fuel designs modifying the thermal behavior. The concept that drives the annular pellet design is reduction of the thermal gradient caused from the high “centerline” temperature at the center of a normal pellet to its surface. By removing the center of the pellet, the new centerline will instead be closer to an edge (whether that be the inner or outer edge), allowing cooling to wick away more heat than the traditional design. In an ideal world, the center of the pellet could also have a heat wick running through it allowing heat to be transported away from throughout the pellet. This increase in heat removal would allow the pellet to be subjected to higher energy output, and thus accelerate the burn-up process. UO 2 has been the staple fuel form of nuclear energy for decades, and replacement of this functionally understood material can only be done with sufficient fuel qualification. Increasing the burnup allows us to characterize new fuels with limited to no operational experience behind them, including carbides, nitrides, and silicide’s. Production of the unique geometries detailed herein will allow us to cut qualification times from multi-decadal experiments into the sub-decadal range, vastly increasing the quantity and quality of information available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Pragmatic Uncertainty Quantification and Propagation in Inverse Estimation of Structural Dynamics Parameters given Material Property Uncertainties and Limited Sensor Data

In this report we demonstrate some relatively simple and inexpensive methods to effectively account for various sources of epistemic lack-of-knowledge type uncertainty in inverse problems. The demonstration problem involves inverse estimation of six parameters of a bolted joint that attaches a kettlebell shaped object to a thick plate. The parameters are efficiently inverted in a modal-based model calibration using gradient-based optimization. Two material properties of the kettlebell are treated as uncertain to within given epistemic uncertainty bounds. We apply and test interval and sparse-sample probabilistic approaches to account for uncertainty in the estimated parameters (and various scalar functionals of the parameters as generic quantities of interest, QOIs) due to uncertainties in the material properties. We also investigate the error effects of limited numbers of vibration sensors (accelerometers) on the kettlebell and plate, and therefore abbreviated excitation/response information in the parameter inversions. We propose and demonstrate a Leave-K-Sensors-Out “cross-prediction” UQ approach to estimate related uncertainties on the parameters and QOI functionals. We indicate how uncertainties from material properties and limited sensors are treated in a combined manner. The economical combined UQ approach involves just three to five samples (i.e. three to five inverse simulations), with no added complication or error/uncertainty from use of surrogate models for affordability. Finally, we describe a related economical UQ approach for handling potential parameter solution non-uniqueness and numerical optimization related precision uncertainties in the estimated parameter values. Indicated further research is identified.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Experimental investigation of a throttlable 15 cm hollow cathode ion thruster

The use of dished high perveance grids on a 15 cm modified SERT 2 thruster is shown to facilitate throttled operation over a beam current range from 60 to 600 mA. Effects of increasing the radial component of the magnetic field in the main discharge chamber and decreasing the dimensions of the cathode discharge region are examined and found to degrade performance to the extent that primary electrons are forced in toward the center-line of the thruster. Studies of the baffle aperture region of two thrusters indicate that the electric potential gradient vector is perpendicular to the local magnetic field lines when the thruster is operating properly. The correlation between the shape of the ion beam current density and that of the ion density at the screen grid within the thruster is shown to be 94%. Additional experimental studies on maximum propellant utilization, plasma ion production cost, neutral density in the cathode discharge region, double ion production in hollow cathode thrusters and thermal flow meter performance are discussed.

Wilbur, P. J.↗
NASA NTRS1973

Effect of bottom topography, eddy diffusivity, and wind variation on circulation in a two-layer stratified lake

The steady-state, wind-driven circulation is calculated in a stratified lake composed of two layers having uniform but unequal densities and eddy diffusivities. The position of the thermocline and the velocities in both layers are calculated from an asymptotic solution of the shallow lake equations when the Ekman number in the epilimnion (upper layer) is of order one but the ratio of hypolimnion (lower layer) to epilimnion eddy diffusivities is much less than one. Large differences in the thermocline shape and the velocities occur between the solution for uniform wind stress and the one for unit order wind stress gradients. For the latter solution the hypolimnion eddy diffusivity magnitude and the bottom topography have a large and important effect.

Gedney, R. T.↗
NASA NTRS1995

Contribution to the optimal shape design of two-dimensional internal flows with embedded shocks

We explore the practicability of optimal shape design for flows modeled by the Euler equations. We define a functional whose minimum represents the optimality condition. The gradient of the functional with respect to the geometry is calculated with the Lagrange multipliers, which are determined by solving a co-state equation. The optimization problem is then examined by comparing the performance of several gradient-based optimization algorithms. In this formulation, the flow field can be computed to an arbitrary order of accuracy. Finally, some results for internal flows with embedded shocks are presented, including a case for which the solution to the inverse problem does not belong to the design space.

Iollo, Angelo↗
NASA NTRS1995

Electromagnetic scattering analysis of a three-dimensional-cavity-backed aperture in an infinite ground plane using a combined finite element method/method of moments approach

A combined finite element method/method of moments (FEM/MoM) approach is used to analyze the electromagnetic scattering properties of a three-dimensional-cavity-backed aperture in an infinite ground plane. The FEM is used to formulate the fields inside the cavity, and the MoM (with subdomain bases) in both spectral and spatial domains is used to formulate the fields above the ground plane. Fields in the aperture and the cavity are solved using a system of equations resulting from the combination of the FEM and the MoM. By virtue of the FEM, this combined approach is applicable to all arbitrarily shaped cavities with inhomogeneous material fillings, and because of the subdomain bases used in the MoM, the apertures can be of any arbitrary shape. This approach leads to a partly sparse and partly full symmetric matrix, which is efficiently solved using a biconjugate gradient algorithm. Numerical results are presented to validate the analysis.

Reddy, C. J.↗
NASA NTRS2021

Satellite Remote Sensing Observations of Trans-Atlantic Dust Transport and Deposition: A Multi-Sensor Analysis

We analyze the decade-long (2007-2016) record of aerosol measurements from four distinctive sensors, namely CALIOP, MODIS, MISR, and IASI, to quantify the trans-Atlantic dust transport and deposition. These satellite sensors use different techniques to characterize particle size and shape properties; and we have developed sensor-specific methods (broadly categorized into size-based and shape-based method) to derive dust optical depth (DOD). The size-based DOD from MODIS and IASI generally agrees better with AERONET-derived DOD than the shape-based DOD from CALIOP and MISR does. Overall, the shape-based DOD is smaller than the size-based DOD by about 25%, which is consistent with distinctive ways of accounting for coarse-spherical particles of dust-pollution internal mixture. While such dust-pollution mixtures are counted as dust in the size-based DOD, they are excluded in the shape-based DOD. DOD is not a good proxy for dust deposition. Instead, the dust deposition depends strongly on the gradient of DOD on a monthly basis and can be derived by calculating the meridional and zonal dust mass flux based on the three-dimensional distributions of dust. Among the remote sensing measurements, difference in dust deposition is smaller than that of DOD, suggesting that different satellites characterize the DOD gradient more consistently than DOD itself. Satellite measurements of dust deposition and DOD also provide an accurate estimate of the dust loss frequency (LF) that measures how efficient the dust is removed from the atmosphere. We found that these remote sensing measurements yield similar LF values of 0.078 –0.102 d-1, which however is factors of 2-5 smaller than model simulations. This analysis provides valuable insights into potential deficiencies in models’ emission and transport/removal processes and hence helps guide model improvement.

remote sensing↗
DOE OSTI2021

High flux boiling heat transfer enhancement using triangle shaped vertical walls in two-phase microchannel heat exchangers

Two-phase heat exchangers allow heat loads to be dissipated with lower fluid flow rates without inducing thermal gradients, and high surface area microchannels further reduce the thermal resistance of these heat exchangers. However, there are many practical limits that prevent widespread use of two-phase microchannel heat exchangers. Many authors have looked to address these practical limits by treating the channel floors or creating artificial nucleation sites. Few works look to enhance heat transfer in the walls of channels with hydraulic diameters below 200 μm. This work is aimed at increasing the efficacy of high aspect ratio silicon microchannels that use R134a as the working fluid. The 125 parallel channels are 40 μm wide, 200 μm deep, and 2 mm long and contain a 15 μm wide, 150 μm long restriction at each inlet. The thermal resistance and dryout characteristics of traditional, plain-walled channels are compared to channels with 2 area enhancement patterns each of which increases the area available for heat transfer by 41%. Furthermore, a resistance to dryout is seen in both area-enhanced test sections, and a 15.5% reduction in thermal resistance is shown for the area-enhanced part with a 6 μm triangular pitch. Only a 7.8% performance increase is seen in the test section with a 3 μm pitch. Both area enhancement patterns induce an average 27% increase in pressure drop for all tests. Heat transfer characteristics are solved for using a finite element analysis (FEA) model. By examining heat transfer coefficients and local heat fluxes, it is postulated that the area enhancement reduces heat transfer coefficients and the performance increase is caused by increasing the number of active nucleation sites.

42 ENGINEERING↗
NASA NTRSDate not supplied

Investigation of a Smooth Local Correlation-based Transition Model in a Discrete-Adjoint Aerodynamic Shape Optimization Algorithm

A smooth local correlation-based transition model is fully coupled to a RANS-based Newton-Krylov flow solver and discrete-adjoint gradient-based optimization algorithm. The free-transition optimization framework is evaluated using lift-constrained drag minimizations of airfoils at design conditions ranging from light to single-aisle aircraft and an infinite swept wing at design conditions representative of a transonic strut-braced wing aircraft. The impact of the streamwise grid resolution on the ability of the optimization algorithm to delay boundary-layer transition is investigated, with the results demonstrating that streamwise grid resolution requirements increase as the transition length decreases with increasing Reynolds number. The optimization problem at the light aircraft design conditions is demonstrated to be multi-modal, with the optimization algorithm producing two distinct designs: one with a thin, reflexed trailing edge and steep pressure recovery regions, the other with increased aft loading, with the latter design outperforming the former. A drag minimization of an airfoil at transonic design conditions demonstrates that the optimization algorithm successfully trades a decrease in viscous drag by delaying boundary-layer transition with an increase in wave drag, while the drag minimization of an infinite swept wing demonstrates the capability of the optimizational gorithm to delay both Tollmien-Schlichting and stationary crossflow instabilities.

AATT↗
NASA NTRS2000

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗
NASA NTRS1943

Determination of general relations for the behavior of turbulent boundary layers

An analysis has been made of a considerable amount of data for turbulent boundary layers along wings and bodies of various shapes in order to determine the fundamental variables that control the development of turbulent boundary layers. It was found that the type of velocity distribution in the boundary layer could be expressed in terms of a single parameter. This parameter was chosen as the ratio of the displacement thickness to the momentum thickness of the boundary layer. The variables that control the development of the turbulent boundary layer apparently are: (1) the ratio of the nondimensional pressure gradient, expressed in terms of the local dynamic pressure outside the boundary layer and boundary-layer thickness, to the local skin-friction coefficient and (2) the shape of the boundary layer. An empirical equation has been developed in terms of these variables that, when used with the momentum equation and the skin-friction relation, makes it possible to trace the development of the turbulent boundary layer to the separation point.

Von Doenhoff, Albert E↗
DOE OSTI2016

EME 192 Report

I spent the quarter working in Lawrence Livermore National Laboratory’s (LLNL) Materials Engineering Division. The group I have been working with (I’ve been here for two summers already) focuses on advanced manufacturing techniques such as stereolithography, electrophoretic deposition, and the printing of silicon based inks. Part of the goal of what is done in our group is to create designer materials not by altering the composition but by altering the micro-architecture. Our technology can create shapes that are not possible with traditional manufacturing techniques. This allows us to create structures that are light, yet very strong and stiff. It also allows us to create materials with property gradients. In other words, we can make structures and parts that are stronger in some locations than others. I have been working with electrophoretic deposition for the duration of my stay and have focused on advancing the technology from a thin-film technique to a true additive manufacturing paradigm. Put succinctly, electrophoretic deposition is the deposition of particles in suspension with electric fields. Particles have a potential on the surface which allows them to be driven to an electrode using an electric field. The particles then deposit onto the conductive regions of the substrate, traditionally, the entire surface. Electrophoretic deposition is powerful in that it can handle a wide variety of materials (ceramics, metals, bacteria), create material gradients in the deposits, and create layered deposition of multiple materials. A drawback of traditional electrophoretic deposition is that patterned deposits are only possible with a non-reconfigurable patterned electrode. A technique was developed at LLNL that allows for the arbitrary patterning of the electric field using photoconductive electrodes and light. This way, you can create interesting shapes and reconfigure the pattern of the deposit using the same electrode. A photoconductive electrode is made by hydrothermally growing titania nanorods onto a transparent current collector. A photomask is used to block incoming some light and only allow the desired pattern of light through. The photoconductive electrode then activates when and where the light hits, once an electric field is applied. Particles will migrate to the areas of illumation and deposit.

36 MATERIALS SCIENCE↗
NASA NTRS2007

Simple and Efficient Numerical Evaluation of Near-Hypersingular Integrals

Recently, significant progress has been made in the handling of singular and nearly-singular potential integrals that commonly arise in the Boundary Element Method (BEM). To facilitate object-oriented programming and handling of higher order basis functions, cancellation techniques are favored over techniques involving singularity subtraction. However, gradients of the Newton-type potentials, which produce hypersingular kernels, are also frequently required in BEM formulations. As is the case with the potentials, treatment of the near-hypersingular integrals has proven more challenging than treating the limiting case in which the observation point approaches the surface. Historically, numerical evaluation of these near-hypersingularities has often involved a two-step procedure: a singularity subtraction to reduce the order of the singularity, followed by a boundary contour integral evaluation of the extracted part. Since this evaluation necessarily links basis function, Green s function, and the integration domain (element shape), the approach ill fits object-oriented programming concepts. Thus, there is a need for cancellation-type techniques for efficient numerical evaluation of the gradient of the potential. Progress in the development of efficient cancellation-type procedures for the gradient potentials was recently presented. To the extent possible, a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. However, since the gradient kernel involves singularities of different orders, we also require that the transformation leaves remaining terms that are analytic. The terms "normal" and "tangential" are used herein with reference to the source element. Also, since computational formulations often involve the numerical evaluation of both potentials and their gradients, it is highly desirable that a single integration procedure efficiently handles both.

Fink, Patrick W.↗