SEARCH · Engineering Papers
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.
A shock tube's driver gas can apparently provide flow with turbulent bursts. The fluctuations are interpreted using a boundary layer model of contact surface flow and results form a kinetic theory of turbulence. With this, a lower limit of 4 on the ratio of maximum to minimum turbulent intensities in contact surface instabilities has been estimated.
Collision Detection for Spacecraft Proximity Operations This thesis describes the development of a new collision detection algorithm to be used when two spacecraft are operating in the same vicinity. The two spacecraft are modelled as unions of convex polyhedra, where the polyhedron resulting from the union may be either convex or nonconvex. The relative motion of the two spacecraft is assumed to be such that one vehicle is moving with constant linear and angular velocity with respect to the other. The algorithm determines if a collision is possible and, if so, predicts the time when the collision will take place. The theoretical basis for the new collision detection algorithm is the C-function formulation of the configuration space approach recently introduced by researchers in robotics. Three different types of C-functions are defined that model the contacts between the vertices, edges, and faces of the polyhedra representing the two spacecraft. These C-functions are used to formulate three "collision" conditions. The first of these conditions limits the points representing potential collisions to the zeros of the C-functions. The new algorithm is fundamentally a search for the smallest zero of any C-function that satisfies the second and third collision conditions. The C-functions are shown to be transcendental functions of time for the assumed trajectory of the moving spacecraft. The zeros of these functions cannot be expressed in dosed form. Therefore, numerical search procedures are developed to find aLl of the zeros of a C-function in specified bounded intervals of time. These bounded intervals of time are found by examining the second and third collision conditions. The capabilities of the new algorithm are demonstrated for several example cases. These include examples of collisions determined by zeros of each of the three different types of C-functions. In addition to predicting the time of first contact of the polyhedra, the algorithm identifies the features of the two polyhedra that are touching at this time. The new collision detection algorithm is the first such algorithm that is capable of solving the collision detection problem exactly for the case where the moving object has constant linear and angular velocities. This is a significant improvement on previous collision detection algorithms described in the literature. In particular, the ability to handle constant angular velocity represents a more realistic type of rotational motion than those which have been used in other algorithms.
In this manuscript we address the challenges associated with the ability to predict radiation sensitivity associated with exposure to either cosmic radiation or X-rays in a population study, by monitoring DNA damage sensing protein 53BP1 forming small nuclear radiation-induced foci (RIF) as a surrogate biomarker of DNA double strand breaks (DSB). 76 primary skin fibroblasts were isolated from 10 collaborative cross strains and five reference inbred mice (C57Bl/6, BALB/CByJ, B6C3, C3H and CBA/CaJ) and exposed to three different charged nuclei of increasing LET (350 MeV/n Si, 350 MeV/n Ar and 600 MeV/n Fe) and X-ray. Our data brings strong evidence against the classic "contact-first" model where DSBs are assumed to be immobile and repaired at the lesion site. In contrast, our model suggests nearby DSBs move into single repair unit characterized by large RIF before the repair machinery kicks in. Such model has the advantage of being much more efficient molecularly but is poorly suited to deal with cosmic radiation, where energy is concentrated along the particle trajectory, inducing a large density of DSBs along each particle track. In accordance with this model, RIF quantification after X-ray exposition showed a saturated dose response for early time points post-irradiation for all strains. Similarly, the high-LET response showed that RIF number matched the number of track per cell, not the number of expected DSB per cell (1). At the temporal level, we noted that the percentage of unrepaired high-LET tracks over a 48 hour time-course increased with LET, confirming that the DNA repair process becomes more difficult as more DSB coalesce into single RIF. There was also good agreement between persistent RIF levels measured in-vitro in the primary skin cultures and survival levels of T-cells and B-cells collected in blood samples from 10 CC strains 24 hours after 0.1 Gy whole-body dose of X-ray. This suggests that persistent RIF 24 hour post-IR is a good surrogate in-vitro biomarker for in-vivo radiation toxicity. Finally, at the genomic level, large differences in repair rates between strains for high-LET allowed us to identify suggestive genetic loci associated with radiation sensitivity. Interestingly, the two highest LETs provided the most strain variation with a common locus on Chromosome 10 highly enriched for DNA repair associated genes we discussed in detail.
The Point-Contact Solar Cell (PCSC) geometry has proven very effective for silicon cells. To date, it has not been implemented in III-V materials. In addition, modeling such a geometry is very difficult because of its three-dimensional nature. We have developed a three-dimensional finite element modeling code (FIESTA ROC). In this paper, we present results from a three-dimensional modeling study of InP point-contact solar cells.
The paper proposes a methodology to fill a gap in the context of mechanistic particle resuspension models related to the adhesion and lift force lever arms ( a ): to date, only one experimental observation exists from a centrifuge experiment with a specific particle/surface pair which found r P / a = 100 and assumed that the contact points were symmetric with respect to the particle center. However, the Reeks and Hall’s (2001) model is highly sensitive to this parameter, which therefore cannot be used as-is for different particle diameters or materials. This paper proposes a methodology to compute the adhesion and lift force lever arms ( a A and a L , respectively) based on Atomic Force Microcopy (AFM) surface roughness measurements and contact mechanics arguments. In practice, a spherical, perfectly smooth particle is brought in contact (computationally) with the measured surface roughness and is then rotated by the point of contact to determine the two or more asperities over which the particle would sit on the surface. A method to compute the asperity radius of curvature is devised to compute the asperity deformation and the correct partitioning of the total deformation between the particle and the asperity is presented. For the glass surface under consideration, the asperity deformation was negligible, i.e., it was below the van der Waals radius, which is the limit of applicability of contact mechanics. The Reeks and Hall’s (2001) model is then extended to account for the joint probability of a A and a L . Significant differences in particle resuspended after a day in common atmospheric conditions are found between the approach that uses the lever-arm distributions or their mean values. Also, important differences are found when comparing the lever-arm distribution and r P / a = 100. Hence, one should be cautious of using r P / a = 100 for surface/particle pairs differing from those in the original study. Finally, limits of applicability of the theories used are presented together with considerations regarding how to extend the method to non-spherical particles and different materials.
Integrated SO x and NO x removal technology using a direct contact cooler is a promising alternative for cleaning coal combustion flue gas, especially for pressurized combustion systems. Past investigations of integrated SO x and NO x removal have been limited to low-temperature processes, where the flue gas has been cooled upstream of the removal unit, and the moisture in the flue gas has been condensed out. Here this work presents a novel design that recovers the flue gas heat while simultaneously removing SO x and NO x in a single reactive-absorption column – direct contact cooler. A validated model was used to evaluate different direct contact cooler designs. Modeling results suggested that a counter-current column with a single water inlet can scrub 83% of the NO x and 96.5% of the SO x in 115 s, while an optimized design with multiple water inlets enhances the scrubbing of NO x by 9 percentage points and SO x by 3.5 percentage points, due to a 34% increase in residence time. It is also observed that the liquid-phase reaction between absorbed SO x and NO x plays a significant role in the removal of SO 2 from the flue gas for high-temperature scrubbing compared with low-temperature scrubbing. The design reinforces the potential of high-temperature SO x -NO x removal and heat recovery and suggests a means of reducing capital costs for such columns.
We develop a model showing how reflection from both the front and back contacts of thermophotovoltaic (TPV) cells affects the system efficiency. From this foundation we determine the design rules for front grids for TPV cells, and assess the consequences of those rules on cell design and system performance. We show that these design rules are different and more complex than the corresponding design rules for solar photovoltaic (PV) cells, and can lead to substantially different grid design and system performance than if the simpler solar PV design rules were used.
This research aims to study the effects of femtosecond laser treatment on surface texture and wettability properties. The first portion of this research studies 2024-T3 aluminum alloy concerning surface texture (e.g., roughness), wettability, and surface free energy. SEM images help aid the researchers in detecting minuscule changes in surface texture. Profilometry extracts topographical data from the surface. A contact angle goniometer measures contact angle, which follows the sessile drop method. Owens, Wendt, Rabel, and Kaelble's method calculates the surface free energy in conjunction with the contact angle. Femtosecond laser modification produces relatively uniform surface morphology at the macroscale. When looking at the microscale, surface morphology differs depending on the interplay between material melting and self-assembly, where the surface becomes rougher after laser treatment. Contact angle results show that all but the control sample and sample 7 exhibit hydrophobic behavior, i.e., contact angle >90°, with Sample 1 having the most significant angle of 115°. The effects of laser parameters on surface energy is apparent. More work is needed to establish correlations among these three quantities: contact angle, surface energy, and laser parameters. The second portion of this research studies aluminum alloy 6061 as the subject material with the primary goal of creating a superhydrophobic surface. A central composite design creates the laser surface texturing design of experiments. SEM images help aid the researchers in detecting minuscule changes in surface texture. A 3D optical profiler extracts topographical data from the surface. The sessile drop method incorporates a contact angle goniometer to measure the contact angle. The response surface method builds a second-order polynomial model for the contact angle and obtains optimized parameters to maximize the contact angle. This research shows that laser surface texturing can generate a wide range of surface profiles and roughness values with geometric features ranging from hundreds of μm to submicron. All three laser parameters (pulse energy, pulse duration, repetition rate) affect surface roughness and contact angle to some degree. To quantify the relationship between the contact angle, pulse energy, and pulse repetition rate, a response surface model for Al 6061 is identified and used to find the optimal conditions of E=214 μJ, tp=10 ps, and fp=2427 Hz with a predicted maximum contact angle of 161°. A confirmation experiment produces a contact angle of 168°, in good agreement with the predicted value.
A face gear drive with a spur involute pinion is considered. The generation of the face gear is based on application of a grinding or cutting worm whereas the conventional method of generation is based on application of an involute shaper. An analytical approach is proposed for the determination of: (1) the worm thread surface; (2) avoidance of singularities of the worm thread surface, (air) dressing of the worm; and (3) determination of stresses of the face-gear drive. A computer program for simulation of meshing and contact of the pinion and face-gear has been developed. Correction of machine-tool settings is proposed for reduction of the shift of the bearing contact caused by misalignment. An automatic development of the model of five contacting teeth has been proposed for stress analysis. Numerical examples for illustration of the developed theory are provided.
A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.
The upstream of bioenergy industry has suffered from unreliable operations of granular biomass feedstocks in handling equipment. Computational modeling, including continuum-mechanics models and discrete-particle models, offers insightful understandings and predictive capabilities on the flow of milled biomass and can assist equipment design and optimization. This paper presents a benchmark study on the fidelity of the continuum and discrete modeling approaches for predicting granular biomass flow. We first introduce the constitutive law of the continuum-mechanics model and the contact law of the coarse-grained discrete-particle model, with model parameters calibrated against laboratory characterization tests of the milled loblolly pine. Three classical granular material flow systems (i.e., a lab-scale rotating drum, a pilot-scale hopper, and a full-scale inclined plane) are then simulated using the two models with the same initial and boundary conditions as the physical experiments. The close agreement of the numerical predictions with the experimental measurements on the hopper mass flow rate, the hopper critical outlet width, the material stopping thickness on the inclined plane, and the dynamic angle of repose, clearly indicates that the two methods can capture the critical flow behavior of granular biomass. The qualitative comparison shows that the continuum-mechanics model outperforms in parameterization of materials and wall friction, and large-scale systems, while the discrete-particle model is more preferred for discontinuous flow systems at smaller scales. Industry stakeholders can use these findings as guidance for choosing appropriate numerical tools to model biomass material flow in part of the optimization of material handling equipment in biorefineries.
An adaptive force control algorithm for small celestial body sampling for a variety of surface properties is developed. The control algorithm consists of an adaptive controller combined with feedback linearization. When a spacecraft makes contact with the surface, it must maintain a desired contact force in order to capture a sample. The properties of the surface are unknown or uncertain before contact with the surface is made. The adaptive controller performs system identification online to create an input-output model of the feedback linearized system. From the input-output model a block observable canonical form is realized and the control input is determined by model predictive control (MPC) to maintain a desired contact force in spite of the unknown surface properties. The approach is applied to a variety of surface properties with linear and nonlinear contact models.
This study numerically investigated the effects of various contact strategies on the thermal hydraulic behavior within a structured bed of 100 explicitly modeled pebbles. Four contact strategies and two thermal hydraulic conditions were considered. The strategies to avoid contact singularities include decreasing the pebble diameter, increasing the pebble diameter, bridging the pebble surfaces near the contact region, and capping the pebble surfaces near the contact region. One strategy, Strategy 3a, which involves bridging with a cylinder equal to 10% of the pebble diameter, was selected as the baseline strategy because it addressed the contact singularity while minimizing the geometric changes that affect the bed porosity. The two thermal hydraulic conditions were full-power operation (Case 1) and pressurized loss of forced cooling or PLOFC (Case 2). Simulations of the conjugate heat transfer within the structured bed were performed using the Reynolds-averaged Navier–Stokes approach with the realizable k-ϵ turbulence model and two-layer all y + wall treatment. The thermal-fluid quantities of interest were compared between the contact strategies for each case. In Case 1, the hydraulic behavior was sensitive to the contact strategy, with large differences in the pressure drop (30%) and volume-average velocity (4%). The thermal behavior was not sensitive, with less than a 0.5% difference across the strategies. To better understand the separate effects of each heat transfer mode, Case 2 was divided into the following subcases: conduction (Case 2a); conduction/radiation (Case 2b); and conduction/radiation/convection (Case 2c). Case 2a represents an early phase of the PLOFC transient. Case 2b represents an intermediate phase of the PLOFC transient, with the pebble temperatures sufficiently high for the radiative heat transfer to be non-negligible. Case 2c represents a late phase of the PLOFC transient after the establishment of the natural circulation of the heat transfer fluid. For Case 2, large differences in the contact strategy were observed only in Case 2a with only conduction. The difference in the maximum pebble temperature was 23% in Case 2a, 2% in Case 2b, and 0.3% in Case 2c.
In this paper, we propose that the manipulator impact control problem be approached from a stochastic optimal control perspective. The reason is that not only is such approach be able to model uncertainties in contact environment, force sensing, as well as manipulator dynamics, the controllers obtained is optimally robust in terms of performance. This result is verified by analyses and simulations.
Device modeling is performed to propose > 25% efficient industry-compatible rear junction double-side passivated contacts solar cell structure with full area p-TOPCon on the rear and selective area n-TOPCon under the front grid pattern (selective TOPCon). Here, this design enables the use of thicker TOPCon (>100nm) on the front for traditional screen-printed contacts without incurring metal-induced damage, high parasitic absorption loss, and compromise in lateral transport or carrier collection on the front. Rear junction design with appropriate bulk lifetime and resistivity combination eliminate the need for heavy doping in the front field region because carriers can flow through the bulk Si without appreciable FF loss. High VOC is maintained because high-quality Si surface passivation in the field region by Al 2 O 3 /SiN gives J 0 comparable to the TOPCon. Our device modeling specifies the practically achievable properties and parameters for each region, including full area rear p-TOPCon, selective area front n-TOPCon, bulk and contacts, to achieve 25.4% efficiency screen-printed bifacial rear junction selective TOPCon cells.
A quasi-two-dimensional analytical model is developed to account for vertical and horizontal current flow in and adjacent to a square ohmic contact between a metal and a thin semiconducting strip which is wider than the contact. The model includes side taps to the contact area for voltage probing and relates the 'apparent' interfacial resistivity to the (true) interfacial resistivity, the sheet resistance of the semiconducting layer, the contact size, and the width of the 'flange' around the contact. This relation is checked against numerical simulations. With the help of the model, interfacial resistivities of ohmic contacts to GaAs were extracted and found independent of contact size in the range of 1.5-10 microns.
The hydrodynamic characteristics of a flying-boat incorporating a low-drag, planing-tail hull were determined from model tests made in Langley tank number 2 and compared with tests of the same flying boat incorporating a conventional-type hull. The planing-tail model, with which stable take-offs were possible for a large range of elevator positions at all center-of-gravity locations tested, had more take-off stability than the conventional model. No upper-limit porpoising was encountered by the planing-tail model. The maximum changes in rise during landings were lower for the planing-tail model than for the conventional model at most contact trims, an indication of improved landing stability for the planing-tail model. The hydrodynamic resistance of the planing-tail hull was lower than the conventional hull at all speeds, and the load-resistance ratio was higher for the planing-tail hull, being especially high at the hump. The static trim of the planing-tail hull was much higher than the conventional hull, but the variation of trim with speed during take-off was smaller.