Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative dynamics”

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 757 records · Page 42

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. II. Entropy production and irreversibility of RG flows

Herein we demonstrate that the reformulation of renormalization group (RG) flow equations as nonlinear heat equations has severe implications on the understanding of RG flows in general. We demonstrate by explicitly constructing an entropy function for a zero-dimensional Z 2 -symmetric model that the dissipative character of generic nonlinear diffusion equations is also hard-coded in the functional RG equation. This renders RG flows manifestly irreversible, revealing the semigroup property of RG transformations on the level of the flow equation itself. Additionally, we argue that the dissipative character of RG flows, its irreversibility and the entropy production during the RG flow may be linked to the existence of a so-called C– / A-function. In total, this introduces an asymmetry in the so-called RG time—in complete analogy to the thermodynamic arrow of time—and allows for an interpretation of infrared actions as equilibrium solutions of dissipative RG flows equations. The impossibility of resolving microphysics from macrophysics is evident in this framework. Furthermore, we directly link the irreversibility and the entropy production in RG flows to an explicit numerical entropy production, which is manifest in diffusive and non-linear partial differential equations (PDEs) and a standard mathematical tool for the analysis of PDEs. Using exactly solvable zero-dimensional Z 2 -symmetric models, we explicitly compute the (numerical) entropy production related to the total variation nonincreasing property of the PDE during RG flows toward the infrared limit. Finally, we discuss generalizations of our findings and relations to the C– / A-theorem as well as how our work may help to construct truncations of RG flow equations in the future, including numerically stable schemes for solving the corresponding PDEs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Protein Adhesion on Semi-Fluorinated Polystyrene Surfaces in Static and Dynamic Measurements

Reducing protein adhesion is a critical strategy in fouling-resistant material innovation, with broad applications spanning biomedical and healthcare devices, biosensors, industrial and environmental systems, and other important technological domains. Here, in this study, we elucidated protein adhesion behavior on polystyrene-based thin films by neutron reflectometry (NR) and quartz crystal microbalance with dissipation (QCM-D), using both lysozyme and bovine serum albumin (BSA) as model proteins. To this end, semifluorinated polystyrene thin films with gradient wettability and surface energy were fabricated through dry processing using plasma oxidation and gas-phase deposition. Although it is believed that a fully fluorinated alkyl chain offers extremely low surface energy, thus rejecting foulants, and has been used in many fouling-resistant surface designs, enhanced protein–surface interactions were observed consistently in NR and QCM-D results, due to the combined effects of surface morphology and chemistry. On the contrary, depositing shorter fluorinated silane onto a hydrophilic PS surface contributed to a more homogeneous nanoscale fluorine coating, resulting in less initial protein adsorption and improved surface recovery. Comparative analysis of proteins with different sizes on the nanopatterned semifluorinated surface revealed the influence of molecular characteristics on surface interactions. Lysozyme, being smaller and more compact, showed faster adsorption kinetics and higher surface coverage but largely reversible binding, whereas BSA, with its larger and more flexible structure, formed broader and more stable interfacial layers. This study fills the gap in understanding protein adhesion within the range of hydrophobicity (water contact angle ∼90°), as current strategies often associate with extreme hydrophilic and superhydrophobic surfaces due to hydration or low-surface-energy rejection mechanisms, respectively. It also provides in-depth insights into current combinatorial fouling-resistant surface design.

Yuan, Yue [Oak Ridge National Laboratory (ORNL), O↗

Dynamical models for a spacecraft idealized as a set of multi-hinged rigid bodies

A brief description is presented of a canonical set of equations which governs the behavior of an n-body spacecraft. General results are given for the case in which the spacecraft is modeled in terms of n rigid bodies connected by dissipative elastic joints. The final equations are free from constraint torques and involve only r variables (r is the number of degrees of freedom of the system). An advantage which accompanies the elimination of the constraint torques is a decrease in the computer run time (especially when n is large).

Larson, V.↗

Double-beam optical method and apparatus for measuring thermal diffusivity and other molecular dynamic processes in utilizing the transient thermal lens effect

A sample material was irradiated by relatively high power, short pulses from a dye laser. Energy from the pulses was absorbed by the sample material, thereby forming a thermal lens in the area of absorption. The pulse repetition rate was chosen so that the thermal lens is substantially dissipated by the time the next pulse reaches the sample material. A probe light beam, which in a specific embodiment is a relatively low power, continuous wave (cw) laser beam, irradiated the thermal lens formed in the sample material. The intensity characteristics of the probe light beam subsequent to irradiation of the thermal lens is related to changes in the refractive index of the sample material as the thermal lens is formed and dissipated.

Moacanin, J.↗

Shockwave dissipation by interface-dominated porous structures

The advent of additive manufacturing (AM) has enabled topological control of structures at the micrometer scale, transforming the properties of polymers for a variety of applications. Examples include tailored mechanical responses, acoustic properties, and thermal properties. Porous polymer materials are a class of materials used for shock and blast mitigation, yet they frequently possess a lack of structural order and are largely developed and evaluated via trial-and-error. Here, we demonstrate control of shockwave dissipation through interface-dominated structures prepared by AM using 2-photon polymerization. A fractal structure with voids, or free surfaces, arranged less than 100 μm apart, allows for rarefaction interactions on the timescale of the shockwave loading. Simulations and dynamic x-ray phase contrast imaging experiments show that fractal structures with interfaces assembled within a “critical” volume reduce shockwave stress and wave velocity by over an order of magnitude within the first unit cell.

36 MATERIALS SCIENCE↗

Development of a precision, wide-dynamic-range actuator for use in active optical systems

The design, operation, and performance of a wide-dynamic-range optical-quality actuator are discussed. The actuator uses a closed-loop control system to maintain accurate positioning and has an rms noise performance of 20 nm. A unique force offloading mechanism allows the actuator coil to dissipate less than 3 mW under quiescent conditions. The operation of an experimental segmented optical system that uses 18 of the actuators is examined to show how they are integrated into an actual system.

Lorell, K. R.↗

Prominence formation in a coronal loop

A model is presented which depends on the preferential deposition of heating in the legs of a coronal loop and which produces a stable prominence-scale condensation at the loop top. Dynamic stability is attained by the subsequent adjustment of local parallel gravity by a magnetic inversion at the loop (or arcade) apex. A nonlinear numerical simulation of this process, which includes a deep chromosphere, a heating rate with a fixed dissipation length, and full solar gravity is described.

Mok, Y.↗

Calculating Shocks In Chemically Reacting Flows

Total-variation-diminishing dissipation operator results in crisp shocks. Report discusses method for numerical simulation of shocks in hypersonic, chemically reacting flows. Introduces mathematical model of gas of as many as three species for which species-conservation equations coupled with gas-dynamic equations.

Eberhardt, Scott↗

Multistage Schemes with Multigrid for Euler and Navier-Strokes Equations: Components and Analysis

A class of explicit multistage time-stepping schemes with centered spatial differencing and multigrids are considered for the compressible Euler and Navier-Stokes equations. These schemes are the basis for a family of computer programs (flow codes with multigrid (FLOMG) series) currently used to solve a wide range of fluid dynamics problems, including internal and external flows. In this paper, the components of these multistage time-stepping schemes are defined, discussed, and in many cases analyzed to provide additional insight into their behavior. Special emphasis is given to numerical dissipation, stability of Runge-Kutta schemes, and the convergence acceleration techniques of multigrid and implicit residual smoothing. Both the Baldwin and Lomax algebraic equilibrium model and the Johnson and King one-half equation nonequilibrium model are used to establish turbulence closure. Implementation of these models is described.

Swanson, R. C.↗

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing Stress-Adaptive Dense Suspensions Using Dynamic Covalent Chemistry

The non-Newtonian behaviors of dense suspensions are central to their use in technological and industrial applications and arise from a network of particle–particle contacts that dynamically adapt to imposed shear. Reported herein are studies aimed at exploring how dynamic covalent chemistry between particles and the polymeric solvent can be used to tailor such stress-adaptive contact networks, leading to their unusual rheological behaviors. Specifically, a room temperature dynamic thia-Michael bond is employed to rationally tune the equilibrium constant (K eq ) of the polymeric solvent to the particle interface. It is demonstrated that low Keq leads to shear thinning, while high K eq produces antithixotropy, a rare phenomenon where the viscosity increases with shearing time. It is proposed that an increase in K eq increases the polymer graft density at the particle surface and that antithixotropy primarily arises from partial debonding of the polymeric graft/solvent from the particle surface and the formation of polymer bridges between particles. Thus, the implementation of dynamic covalent chemistry provides a new molecular handle with which to tailor the macroscopic rheology of suspensions by introducing programmable time dependence. These studies open the door to energy-absorbing materials that not only sense mechanical inputs and adjust their dissipation as a function of time or shear rate but also can switch between these two modalities on demand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

A compressible near-wall turbulence model for boundary layer calculations

A compressible near-wall two-equation model is derived by relaxing the assumption of dynamical field similarity between compressible and incompressible flows. This requires justifications for extending the incompressible models to compressible flows and the formulation of the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilational part, which is directly affected by these changes. This approach isolates terms with explicit dependence on compressibility so that they can be modeled accordingly. An equation that governs the transport of the solenoidal dissipation rate with additional terms that are explicitly dependent on the compressibility effects is derived similarly. A model with an explicit dependence on the turbulent Mach number is proposed for the dilational dissipation rate. Thus formulated, all near-wall incompressible flow models could be expressed in terms of the solenoidal dissipation rate and straight-forwardly extended to compressible flows. Therefore, the incompressible equations are recovered correctly in the limit of constant density. The two-equation model and the assumption of constant turbulent Prandtl number are used to calculate compressible boundary layers on a flat plate with different wall thermal boundary conditions and free-stream Mach numbers. The calculated results, including the near-wall distributions of turbulence statistics and their limiting behavior, are in good agreement with measurements. In particular, the near-wall asymptotic properties are found to be consistent with incompressible behavior; thus suggesting that turbulent flows in the viscous sublayer are not much affected by compressibility effects.

So, R. M. C.↗

A Modeling Investigation of Thermal and Strain Induced Recovery and Nonlinear Hardening in Potential Based Viscoplasticity

Specific forms for both the Gibb's and the complementary dissipation potentials were chosen such that a complete potential based multiaxial, isothermal, viscoplastic model was obtained. This model in general possesses three internal state variables (two scalars associated with dislocation density and one tensor associated with dislocation motion) both thermal and dynamic recovery mechanisms, and nonlinear kinematic hardening. This general model, although possessing associated flow and evolutionary laws, is shown to emulate three distinct classes of theories found in the literature, by modification of the driving threshold function F. A parametric study was performed on a specialized nondimensional multiaxial form containing only a single tensorial internal state variable (i.e., internal stress). The study was conducted with the idea of examining the impact of including a strain-induced recovery mechanism and the compliance operator, derived from the Gibb's potential, on the uniaxial and multiaxial response. One important finding was that inclusion of strain recovery provided the needed flexibility in modeling stress-strain and creep response of metals at low homologous temperatures, without adversely affecting the high temperature response. Furthermore, for nonproportional loading paths, the inclusion of the compliance operator had a significant influence on the multiaxial response, but had no influence on either uniaxial or proportional load histories.

Arnold, S. M.↗

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗

Theoretical, observational, and isotopic estimates of the lifetime of the solar nebula

There are a variety of isotopic data for meteorites which suggest that the protostellar nebula existed and was involved in making planetary materials for some 10(exp 7) yr or more. Many cosmochemists, however, advocate alternative interpretations of such data in order to comply with a perceived constraint, from theoretical considerations, that the nebula existed only for a much shorter time, usually stated as less than or equal to 10(exp 6) yr. In this paper, we review evidence relevant to solar nebula duration which is available through three different disciplines: theoretical modeling of star formation, isotopic data from meteorites, and astronomical observations of T Tauri stars. Theoretical models based on observations of present star-forming regions indicate that stars like the Sun form by dynamical gravitational collapse of dense cores of cold molcular clouds in the interstellar clouds in the interstellar medium. The collapse to a star and disk occurs rapidly, on a time scale of the order 10(exp 5) yr. Disks evolve by dissipating energy while redistributing angular momentum, but it is difficult to predict the rate of evolution, particularly for low mass (compared to the star) disks which nonetheless still contain enough material to account for the observed planetary system. There is no compelling evidence, from available theories of disk structure and evolution, that the solar nebula must have evolved rapidly and could not have persisted for more than 1 Ma. In considering chronoloically relevant isotopic data for meteorites, we focus on three methodologies: absolute ages by U-Pb/Pb-Pb, and relative ages by short-lived radionuclides (especially Al-26) and by evolution of Sr-87/Sr-86. Two kinds of meteoritic materials-refractory inclusions such as CAIs and differential meteorites (eucrites and augrites) -- appear to have experienced potentially dateable nebular events. In both cases, the most straightforward interpretations of the available data indicate nebular events spanning several Ma. We also consider alternative interpretations, particularly the hypothesis of radically heterogeneous distribution of Al-26, which would avoid these chronological interpretations. The principal impetus for such alternative interpretations seems to be precisely the obviation of the chronological interpretation (i.e., the presumption rather than the inference of a short (less than or equal to 1 Ma) lifetime of the nebula). Astronomical observations of T Tauri stars indicate that the presence of dusty disks is a common if not universal feature, that the disks are massive enough to accomodate a planetary system such as ours, and that at least some persist for 110(exp 7) yr or more. The results are consistent with the time scales inferred from the meteoritic isotopic data. They cannot be considered conclusive with regard to solar nebula time scales, however, in part because it is difficult to relate disk observations to processes that affect meteorites, and in part because the ages assigned for these stars could be wrong by a factor of several in either direction. We conclude that the balance of available evidence favors the view that the nebula existed and was active for at least several Ma. However, because the evidence is not definitive, it is important that the issue be perceived to be an open question, whose answer should be sought rather than presumed.

Podosek, Frank A.↗