Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic boundaries”

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,369 records · Page 76

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗

Overview of the Spirit Mars Exploration Rover Mission to Gusev Crater: Landing Site to Backstay Rock in the Columbia Hills

Spirit landed on the floor of Gusev Crater and conducted initial operations on soil covered, rock-strewn cratered plains underlain by olivine-bearing basalts. Plains surface rocks are covered by wind-blown dust and show evidence for surface enrichment of soluble species as vein and void-filling materials and coatings. The surface enrichment is the result of a minor amount of transport and deposition by aqueous processes. Layered granular deposits were discovered in the Columbia Hills, with outcrops that tend to dip conformably with the topography. The granular rocks are interpreted to be volcanic ash and/or impact ejecta deposits that have been modified by aqueous fluids during and/or after emplacement. Soils consist of basaltic deposits that are weakly cohesive, relatively poorly sorted, and covered by a veneer of wind blown dust. The soils have been homogenized by wind transport over at least the several kilometer length scale traversed by the rover. Mobilization of soluble species has occurred within at least two soil deposits examined. The presence of mono-layers of coarse sand on wind-blown bedforms, together with even spacing of granule-sized surface clasts, suggest that some of the soil surfaces encountered by Spirit have not been modified by wind for some time. On the other hand, dust deposits on the surface and rover deck have changed during the course of the mission. Detection of dust devils, monitoring of the dust opacity and lower boundary layer, and coordinated experiments with orbiters provided new insights into atmosphere-surface dynamics.

Arvidson, R. E.↗

Circulation Control Model Experimental Database for CFD Validation

A 2D circulation control wing was tested in the Basic Aerodynamic Research Tunnel at the NASA Langley Research Center. A traditional circulation control wing employs tangential blowing along the span over a trailing-edge Coanda surface for the purpose of lift augmentation. This model has been tested extensively at the Georgia Tech Research Institute for the purpose of performance documentation at various blowing rates. The current study seeks to expand on the previous work by documenting additional flow-field data needed for validation of computational fluid dynamics. Two jet momentum coefficients were tested during this entry: 0.047 and 0.114. Boundary-layer transition was investigated and turbulent boundary layers were established on both the upper and lower surfaces of the model. Chordwise and spanwise pressure measurements were made, and tunnel sidewall pressure footprints were documented. Laser Doppler Velocimetry measurements were made on both the upper and lower surface of the model at two chordwise locations (x/c = 0.8 and 0.9) to document the state of the boundary layers near the spanwise blowing slot.

Paschal, Keith B.↗

Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications

The need to maintain optimal energy efficiency is critical during the drilling operations performed on future and current planetary rover missions (see figure). Specifically, this innovation seeks to solve the following problem. Given a spring-loaded percussive drill driven by a voice-coil motor, one needs to determine the optimal input voltage waveform (periodic function) and the optimal hammering period that minimizes the dissipated energy, while ensuring that the hammer-to-rock impacts are made with sufficient (user-defined) impact velocity (or impact energy). To solve this problem, it was first observed that when voice-coil-actuated percussive drills are driven at high power, it is of paramount importance to ensure that the electrical current of the device remains in phase with the velocity of the hammer. Otherwise, negative work is performed and the drill experiences a loss of performance (i.e., reduced impact energy) and an increase in Joule heating (i.e., reduction in energy efficiency). This observation has motivated many drilling products to incorporate the standard bang-bang control approach for driving their percussive drills. However, the bang-bang control approach is significantly less efficient than the optimal energy-efficient control approach solved herein. To obtain this solution, the standard tools of classical optimal control theory were applied. It is worth noting that these tools inherently require the solution of a two-point boundary value problem (TPBVP), i.e., a system of differential equations where half the equations have unknown boundary conditions. Typically, the TPBVP is impossible to solve analytically for high-dimensional dynamic systems. However, for the case of the spring-loaded vibro-impactor, this approach yields the exact optimal control solution as the sum of four analytic functions whose coefficients are determined using a simple, easy-to-implement algorithm. Once the optimal control waveform is determined, it can be used optimally in the context of both open-loop and closed-loop control modes (using standard realtime control hardware).

Aldrich, Jack B.↗

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗

Dynamic processes in transient phases during self-assembly of organic semiconductor thin films

Understanding and manipulating crystallization processes has been an important challenge for solution-processed organic thin films, both for fundamental studies and for fabricating thin films with near-intrinsic charge transport properties. We report an in situ X-ray scattering study of the crystallization of 2-decyl-7-phenyl-[1]benzothieno[3,2-b][1]benzothiophene (Ph-BTBT-C 10 ) during its deposition from solution. At temperatures modestly below the smectic-E/crystalline phase boundary, the crystallization goes through a transient liquid crystal state before reaching the final stable crystalline phase. Significant dynamics occur in the first few seconds of the transition, which are observed through fluctuations in the X-ray scattering intensity, and are correlated with the time interval that the transient thin film coexists with the evaporating solvent. The transition to the stable crystalline phase takes minutes or even hours under these conditions, which may be a result of the asymmetry of the molecule. Transient phases are of potential interest for applications, since they can act as a route to self-assembly of organic thin films. However, our observations show that the long-lived monolayer-stacked intermediate state does not act as a template for the bilayer-stacked crystalline phase. Rather, the grain structure is replaced through nucleation, where the nucleation free-energy barrier is related to a potential barrier that prevents molecules to flip their long axis by 180°.

36 MATERIALS SCIENCE↗

Hypersonic Turbulent Boundary-Layer and Free Sheer Database Datasets

A critical assessment and compilation of data are presented on attached hypersonic turbulent boundary layers in pressure gradients and compressible turbulent mixing layers. Extensive searches were conducted to identify candidate experiments, which were subjected to a rigorous set of acceptance criteria. Accepted datasets are both tabulated and provided in machine-readable form. The purpose of this database effort is to make existing high quality data available in detailed form for the turbulence-modeling and computational fluid dynamics communities. While significant recent data were found on the subject of compressible turbulent mixing, the available boundary-layer/pressure-gradient experiments are all older ones of which no acceptable data were found at hypersonic Mach numbers.

Settles, Gary S.↗

Möbius mirrors

An accelerating boundary (mirror) acts as a horizon and black hole analog, radiating energy with some particle spectrum. In this work, we demonstrate that a Möbius transformation on the null coordinate advanced time mirror trajectory uniquely keeps invariant not only the energy flux but the particle spectrum. We clarify how the geometric entanglement entropy is also invariant. The transform allows generation of families of dynamically distinct trajectories, including $\mathcal{PT}$ -symmetric ones, mapping from the eternally thermal mirror to the de Sitter horizon, and different boundary motions corresponding to Kerr or Schwarzschild black holes

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗

Nonlinear Nano‐Imaging of Interlayer Coupling in 2D Graphene‐Semiconductor Heterostructures

The emergent electronic, spin, and other quantum properties of 2D heterostructures of graphene and transition metal dichalcogenides are controlled by the underlying interlayer coupling and associated charge and energy transfer dynamics. However, these processes are sensitive to interlayer distance and crystallographic orientation, which are in turn affected by defects, grain boundaries, or other nanoscale heterogeneities. This obfuscates the distinction between interlayer charge and energy transfer. Here, nanoscale imaging in coherent four‐wave mixing (FWM) and incoherent two‐photon photoluminescence (2PPL) is combined with a tip distance‐dependent coupled rate equation model to resolve the underlying intra‐ and inter‐layer dynamics while avoiding the influence of structural heterogeneities in mono‐ to multi‐layer graphene/WSe 2 heterostructures. With selective insertion of hBN spacer layers, it is shown that energy, as opposed to charge transfer, dominates the interlayer‐coupled optical response. From the distinct nano‐FWM and ‐2PPL tip‐sample distance‐dependent modification of interlayer and intralayer relaxation by tip‐induced enhancement and quenching, an interlayer energy transfer time of τ ET ≈ ( 0 . 35 − 0.15 + 0.65 ) $\tau _{\rm ET} \approx (0.35^{+0.65}_{-0.15})$ ps consistent with recent reports is derived. As a local probe technique, this approach highlights the ability to determine intrinsic sample properties even in the presence of large sample heterogeneity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗

Perfectly Matched Layers and Characteristic Boundaries in Lattice Boltzmann: Accuracy vs Cost

Artificial boundary conditions (BCs) play a ubiquitous role in numerical simulations of transport phenomena in several diverse fields, such as fluid dynamics, electromagnetism, acoustics, geophysics, and many more. They are essential for accurately capturing the behavior of physical systems whenever the simulation domain is truncated for computational efficiency purposes. Ideally, an artificial BC would allow relevant information to enter or leave the computational domain without introducing artifacts or unphysical effects. Boundary conditions designed to control spurious wave reflections are referred to as nonreflective boundary conditions (NRBCs). Another approach is given by the perfectly matched layers (PMLs), in which the computational domain is extended with multiple dampening layers, where outgoing waves are absorbed exponentially in time. Here, in this work, the definition of PML is revised in the context of the lattice Boltzmann method. The impact of adopting different types of BCs at the edge of the dampening zone is evaluated and compared, in terms of both accuracy and computational costs. It is shown that for sufficiently large buffer zones, PMLs allow stable and accurate simulations even when using a simple zeroth-order extrapolation BC. Moreover, employing PMLs in combination with NRBCs potentially offers significant gains in accuracy at a modest computational overhead, provided the parameters of the BC are properly tuned to match the properties of the underlying fluid flow.

97 MATHEMATICS AND COMPUTING↗

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗

Aerothermodynamic Measurement and Prediction for Modified Orbiter at Mach 6 and 10

Detailed heat-transfer rate distributions measured laterally over the windward surface of an orbiter-like configuration using thin-film resistance heat-transfer gauges and globally using the newly developed relative intensity, two-color thermographic phosphor technique are presented for Mach 6 and 10 in air. The angle of attack was varied from 0 to 40 deg, and the freestream Reynolds number based on the model length was varied from 4 x 10(exp 5) to 6 x 10(exp 6) at Mach 6, corresponding to laminar, transitional, and turbulent boundary layers; the Reynolds number at Mach 10 was 4 x 10(exp 5), corresponding to laminar flow. The primary objective of the present study was to provide detailed benchmark heat-transfer data for the calibration of computational fluid-dynamics codes. Predictions from a Navier-Stokes solver referred to as the Langley aerothermodynamic upwind relaxation algorithm and an approximate boundary-layer solving method known as the axisymmetric analog three-dimensional boundary layer code are compared with measurement. In general, predicted laminar heat-transfer rates are in good agreement with measurements.

Micol, John R.↗

Riverbed Temperature and 4D ERT Monitoring Reveals Heterogenous Horizontal and Vertical Groundwater-Surface Water Exchange Flows Under Dynamic Stage Conditions

Groundwater surface water exchange plays a critical role in physical, biological, and geochemical function of coastal and riverine systems. Observing exchange flow behavior in heterogeneous systems is a primary challenge, particularly when flows are governed by dynamic river stage or tidal variations. In this paper we demonstrate a novel application of time-lapse 3D electrical resistivity tomography and temperature monitoring where an array of thermistors installed beneath a riverbed double as resistivity electrodes. We use the array to monitor stage driven exchange flows over a 6-day period in a dynamic, stage-driven high order stream. We present a method for addressing the otherwise confounding effects of the moving river-surface boundary on the raw resistivity data, thereby enabling successful tomographic imaging. Temperature time-series at each thermistor location and time-lapse 3D images of changes in bulk electrical conductivity together provide a detailed description of exchange dynamics over a 10-meter by 45-meter section of the riverbed, to a depth of approximately 5 m. Results reveal highly variable flux behavior throughout the monitoring domain including both horizontal and vertical exchange flows.

54 ENVIRONMENTAL SCIENCES↗

Phase Synchronization and Desynchronization of Structural Response Induced by Turbulent and External Sound

Acoustic and turbulent boundary layer flow loadings over a flexible structure are used to study the spatial-temporal dynamics of the response of the structure. The stability of the spatial synchronization and desynchronization by an active external force is investigated with an array of coupled transducers on the structure. In the synchronous state, the structural phase is locked, which leads to the formation of spatial patterns while the amplitude peaks exhibit chaotic behaviors. Large amplitude, spatially symmetric loading is superimposed on broadband, but in the desynchronized state, the spectrum broadens and the phase space is lost. The resulting pattern bears a striking resemblance to phase turbulence. The transition is achieved by using a low power external actuator to trigger broadband behaviors from the knowledge of the external acoustic load inducing synchronization. The changes are made favorably and efficiently to alter the frequency distribution of power, not the total power level. Before synchronization effects are seen, the panel response to the turbulent boundary layer loading is discontinuously spatio-temporally correlated. The stability develops from different competing wavelengths; the spatial scale is significantly shorter than when forced with the superimposed external sound. When the external sound level decreases and the synchronized phases are lost, changes in the character of the spectra can be linked to the occurrence of spatial phase transition. These changes can develop broadband response. Synchronized responses of fuselage structure panels have been observed in subsonic and supersonic aircraft; results from two flights tests are discussed.

Maestrello, Lucio↗

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗

A gyrokinetic simulation model for 2D equilibrium potential in the scrape-off layer of a field-reversed configuration

The equilibrium potential structure in the scrape-off layer (SOL) of the field-reversed configuration (FRC) can be affected by the penetration of edge biasing applied at the divertor ends. The primary focus of the paper is to establish a formulation that accurately captures both parallel and radial variations of the two-dimensional (2D) potential in SOL. The formulation mainly describes a quasi-neutral plasma with a logical sheath boundary. A full-f gyrokinetic ion model and a massless electron model are implemented in the GTC-X code to solve for the self-consistent equilibrium potential, given fixed radial potential profiles at the boundaries. The first essential point of this 2D model lies in its ability to couple radial and parallel dynamics stemming from resistive currents and drag force on ions. The model successfully recovers the fluid force balance and continuity equations. These collisional effects on 2D potential mainly appear through the density profile changes, modifying the potential through electron pressure gradient. This means an accurate prescription of electron density and temperature profiles is important in predicting the potential structure in the FRC SOL. The Debye sheath potential and the potential profiles applied at the boundaries can be additional factors contributing to the 2D variations in SOL. This comprehensive full-f scheme holds promise for future investigations into turbulent transport in the presence of the self-consistent 2D potential together with the non-Maxwellian distributions and open boundary conditions in the FRC SOL.

Physics↗