Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak form”

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 865 records · Page 48

Mesospheric winds observed by the Kyoto meteor radar

The Kyoto meteor radar has been operated almost continuously from May 1983. By using a large amount of radar data, seasonal variations of mean winds, planetary waves and atmospheric tides can be delinated. Seasonal and latitudinal variation of tidally corrected mean winds including Kyoto observations have been examined in comparison with the CIRA 1972 model. A power spectral density of the zonal wind field for 60 days of data is formed from an autocorrelation function with a maximum time lag of 30 days. A quasi-2-day oscillation was evident in the period from June to September, and seems to be repeated in 1983 and 1983. It had large amplitudes in June and July and became weak in August; then it was again amplified in September. In summer months, the period of the wave was indicated as 2.2 days, and it became shorter in September. At Kyoto (35 deg N), the diurnal tide is the main component of atmospheric tides. The monthly mean amplitude of the eastward component of diurnal tide is plotted for the period from May 1983 to September 1984. It was clearly amplified in summer months at around 92 to 96 km altitudes. This feature was repeated in two successive years except for the maximum amplitude which exceeded 40 m/sec and 30 m/sec in 1983 and 1984, respectively.

Tsuda, T.↗

Progress on a Taylor weak statement finite element algorithm for high-speed aerodynamic flows

A new finite element numerical Computational Fluid Dynamics (CFD) algorithm has matured to the point of efficiently solving two-dimensional high speed real-gas compressible flow problems in generalized coordinates on modern vector computer systems. The algorithm employs a Taylor Weak Statement classical Galerkin formulation, a variably implicit Newton iteration, and a tensor matrix product factorization of the linear algebra Jacobian under a generalized coordinate transformation. Allowing for a general two-dimensional conservation law system, the algorithm has been exercised on the Euler and laminar forms of the Navier-Stokes equations. Real-gas fluid properties are admitted, and numerical results verify solution accuracy, efficiency, and stability over a range of test problem parameters.

Baker, A. J.↗

Fully-Explicit and Self-Consistent Algebraic Reynolds Stress Models

A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the 'weak equilibrium' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the 'weak equilibrium' assumption treat the production-to-dissipation (P/epsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

Girimaji, Sharath S.↗

Nanostructured Composites: Effective Mechanical Property Determination of Nanotube Bundles

Carbon nanotubes naturally tend to form crystals in the form of hexagonally packed bundles or ropes that should exhibit a transversely isotropic constitutive behavior. Although the intratube axial stiffness is on the order of 1 TPa due to a strong network of delocalized bonds, the intertube cohesive strength is orders of magnitude less controlled by weak, nonbonding van der Waals interactions. An accurate determination of the effective mechanical properties of nanotube bundles is important to assess potential structural applications such as reinforcement in future composite material systems. A direct method for calculating effective material constants is developed in the present study. The Lennard-Jones potential is used to model the nonbonding cohesive forces. A complete set of transverse moduli are obtained and compared with existing data.

Saether, E.↗

Creation and Distribution of CAIs in the Protoplanetary Nebula

CaAl rich refractory mineral inclusions (CAIs) found at 1 - 10% mass fraction in primitive chondrites appear to be several million years older than the dominant (chondrule) components in the same parent bodies. A prevalent concern is that it is difficult to retain CAIs for this long against gas-drag-induced radial drift into the sun. We assess a hot inner (turbulent) nebula context for CAI formation, using analytical models of nebula evolution and particle diffusion. We show that outward radial diffusion in a weakly turbulent nebula can prevent significant numbers of CAI-size particles from being lost into the sun for times of 1 - 3 x 10(exp 6) years. To match the CAI abundances quantitatively, we advocate an enhancement of the inner hot nebula in silicate-forming material, due to rapid inward migration of very primitive, silicate and carbon rich, meter-sized objects. 'Combustion' of the carbon into CO would make the CAI formation environment more reduced than solar, as certain observations imply. Abundant CO might also play a role in mass-independent chemical fractionation of oxygen isotopes as seen in CAIs and associated primitive, high-temperature condensates.

Cuzzi, J. N.↗

Is Water a Universal Solvent for Life?

There are strong reasons to believe that the laws, principles and constraints of physics and chemistry are universal. It is much less clear how this universality translates into our understanding of the origins of life. Conventionally, discussions of this topic focus on chemistry that must be sufficiently rich to seed life. Although this is clearly a prerequisite for the emergence of living systems, I propose to focus instead on self-organization of matter into functional structures capable of reproduction, evolution and responding to environmental changes. In biology, most essential functions are largely mediated by noncovalent interactions (interactions that do not involve making or breaking chemical bonds). Forming chemical bonds is only a small part of what living systems do. There are specific implications of this point of view for universality. I will concentrate on one of these implications. Strength of non-covalent interactions must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. This balance, however, is not a natural property of complex chemical systems. Instead, it has to be achieved with the aid of an appropriate solvent for life. In particular, potential solvents for life must be characterized by a high dielectric constant to ensure solubility of polar species and sufficient flexibility of biological structures stabilized by electrostatic interactions. Among these solvents, water exhibits a remarkable trait that it also promotes solvophobic (hydrophobic) interactions between non-polar species, typically manifested by a tendency of these species to aggregate and minimize their contacts with the aqueous solvent. Hydrophobic interactions are responsible, at least in part, for many self-organization phenomena in biological systems, such as the formation of cellular boundary structures or protein folding. Strengths of electrostatic and hydrophobic interactions are similar and can be balanced over a wide range of temperatures, which considerably increases the repertoire of interactions that can be used to modulate biological functions. Some properties of water, e.g. its chemical activity against polymerization reactions, are considered as unfavorable to life. In actuality, this might be a favorable trait because life requires a balance between constructive and destructive processes. For example, molecules synthesized in response to specific conditions must be degraded once these conditions change. Otherwise regulation of biological processes would be virtually impossible. Water might not be the only liquid with favorable properties for supporting life. It has been proposed that formamide, which might be present elsewhere in the universe in sufficient quantities to warrant interest, could be a potential alternative to water for the origin of life. However, this will remain highly hypothetical until it is demonstrated in further studies on its physical, chemical and biological properties it is capable of mediating self-organization of matter and providing proper balance between different types of non-covalent interactions.

Pohorill, Andrew↗

Structural and Dynamical Analysis of the Quiescent Molecular Ridge in the Large Magellanic Cloud

We present a comparison of low-J13CO and CS observations of four different regions in the LMC—the quiescent Molecular Ridge, 30 Doradus, N159, and N113, all at a resolution of ∼3 pc. The regions 30 Dor, N159, and N113 are actively forming massive stars, while the Molecular Ridge is forming almost no massive stars, despite its large reservoir of molecular gas and proximity to N159 and 30 Dor. We segment the emission from each region into hierarchical structures using dendrograms and analyze the sizes, masses, and line widths of these structures. We find that the Ridge has significantly lower kinetic energy at a given size scale and also lower surface densities than the other regions, resulting in higher virial parameters. This suggests that the Ridge is not forming massive stars as actively as the other regions because it has less dense gas and not because collapse is suppressed by excess kinetic energy. We also find that these physical conditions and energy balance vary significantly within the Ridge and that this variation appears only weakly correlated with distance from sites of massive-star formation such as R136 in 30 Dor, which is ∼1 kpc away. These variations also show only a weak correlation with local star formation activity within the clouds.

Molly K Finn↗

Large Impact Features on Ganymede and Callisto as Revealed By Geological Mapping and Morphometry

The icy Galilean satellites are host to a broad range of impact feature morphologies. Hypotheses seeking to explain the diversity of these impact features consider the effects of impact melt, the physical state of the subsurface at the time of impact, and the impactor characteristics. As part of a larger effort to assess the role of these factors in the formation and evolution of these impact features, we have performed topographic and geological mapping of 19 large impact features on Ganymede and Callisto. These are divided into two main morphological groups: craters (subdivided into pit, dome, and anomalous dome craters), and penepalimpsests/palimpsests. The transitions from pit, dome, to anomalous dome craters appear to be size-dependent up to diameters of ∼170 km. The morphologies of pit and dome craters appear to be independent of their age or geologic context. The impacts that formed them only affected a cold, rigid ice layer, with the development of pits and raised annuli on their floors possibly stemming from the evolution of a pocket of impact melt. The subdued rims and floors of anomalous dome craters indicate the increasing effect of a weak, warm ice layer on impact feature morphology with increasing size, but their prominent annuli and pits indicate that mobilization of impact melt is also a factor. The very low topographic relief of older penepalimpsests and palimpsests indicates that their impacts penetrated the ice shell to mobilize very large volumes of pre-existing liquid from a subsurface layer, with little contribution to the final feature morphology from impact melt. Penepalimpsests are distinguished from palimpsests by the higher frequency of concentric ridges within their interiors, indicating a generally more robust state of the subsurface that could better support the rotation and uplift of solid material during impact, even if a crater-like depression could not be supported. A few impact features seem to be transitional between anomalous dome craters and penepalimpsests, and the overlap of anomalous dome craters, penepalimpsests, and palimpsests in terms of diameter as well as age indicates that impactor size and subsurface properties over time are major factors in determining which of these morphologies emerges.

Oliver L White↗

Long-range three-dimensional magnetic structures of the spin S =1 hexamer cluster fedotovite-like A 2 Cu 3 O ( SO 4 ) 3 ( A 2 = K 2 , NaK, Na 2 ): A neutron diffraction study

The crystal and magnetic structures of the spin S = 1 hexamer cluster fedotovite-like A 2 Cu 3 O(SO 4 ) 3 (A 2 = K 2 , NaK, Na 2 ) were studied by neutron powder diffraction at temperatures 1.6–290 K. The crystal structures in all compounds are well refined in the monoclinic space group C2/c. The basic magnetic units of the compounds are copper hexamers, which are coupled by weak superexchange interactions giving rise to three-dimensional long-range magnetic order below 3.0 < T N < 4.7 K. We have found that for A 2 = K 2 and NaK the propagation vector of the magnetic structure is k = [0,0,0], and the coupling of the Cu hexamers is ferromagnetic (FM) along the ab diagonal and antiferromagnetic along the bc diagonal. In contrast, for A = Na the propagation vector is k = [0,1,0], and the Cu hexamers are coupled antiferromagnetically (AFM) along the ab diagonal. Further, the hexamers are formed by three Cu pairs arranged along the b axis. The calculated spin expectation values $\langle s \rangle$ for the simplest symmetric spin Hamiltonian (obtained from inelastic neutron spectroscopy) of the isolated hexamers in the mean field amounted to $\langle s \rangle$ = 3/8 for the side spins Cu1 and Cu2 and $\langle s \rangle$ = 1/4 for Cu3 in the middle. The Cu spins are FM coupled in pairs and AFM between neighboring pairs. The experimental magnetic moments of the Cu 2+ ions turn out to be not completely collinear due to spin frustrations within the weak interhexamer interactions. The sizes of magnetic moments of Cu in the hexamers determined from the diffraction data are in fair agreement with the calculated values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Facet-Dependent Water Inhibition of Alkanol Dehydration on TiO2 via Distinct Water–Alkanol Complexes

Water is ubiquitous in biomass-derived feeds, yet its molecular impact on oxygen-elimination reactions remains poorly understood, particularly for catalysts exposing different facets. Here, we utilize well-defined TiO2 nanocrystals with dominant (101) and (001) facets to reveal a pronounced facet-dependent effect of water, where inhibition for dehydration of isopropanol (IPA) on the TiO2(001) surface is about four times more severe than TiO2(101). Through a combination of in situ solid state NMR, in situ infrared spectroscopy, kinetics studies, and theoretical calculations, we demonstrate that this disparity arises from the formation of distinct alkanol-water complex intermediates. On TiO2(001), IPA undergoes dissociative adsorption to form an isopropoxide-H2O complex that readily drives the surface into a complex-dominated regime. This pathway increases the activation barrier for C–H cleavage by 40 kJ mol-1 by inducing a disordered transition state. In contrast, TiO2(101) favors molecular IPA adsorption with weak hydrogen bonding to water, resulting in a smaller complex formation constant and a much smaller activation barrier increase (25 kJ mol-1). By quantitatively linking facet-dependent complex coverage to transition-state destabilization, this work moves beyond simple site-blocking models and provides a conceptual framework for designing catalysts that remain active in water-containing environments

Hu, Wenda↗

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗

Electron Acceleration during Macroscale Magnetic Reconnection

Here, the first self-consistent simulations of electron acceleration during magnetic reconnection in a macroscale system are presented. Consistent with solar flare observations, the spectra of energetic electrons take the form of power laws that extend more than two decades in energy. The drive mechanism for these nonthermal electrons is Fermi reflection in growing and merging magnetic flux ropes. A strong guide field suppresses the production of nonthermal electrons by weakening the Fermi drive mechanism. For a weak guide field the total energy content of nonthermal electrons dominates that of the hot thermal electrons even though their number density remains small. Our results are benchmarked with the hard x-ray, radio, and extreme ultraviolet observations of the X8.2-class solar flare on September 10, 2017.

79 ASTRONOMY AND ASTROPHYSICS↗

lanl/nubhlight

nubhlight is a code for simulating disks of materials around black holes formed after the in-spiral and merger of two-black holes. It is based on the open-source code bhlight. nubhlight models all relevant physics required to solve the post-merger disk problem. This includes: general relativity in a static spacetime, neutrino and weak physics, relativistic magnetohydrodynamics, and the ability to use open-sourced tabulated equations of state provided by stellarcollapse.org.

Miller, Jonah↗

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗

Mass loss, levitation, accretion, and the sharp-lined features in hot white dwarfs

A study has been conducted of eight white dwarfs, including seven DA and one He-rich types. The study is based on high-resolution observations conducted with the aid of the International Ultraviolet Explorer. Four of the dwarfs show features related to heavy elements which are not interstellar in origin. It is tentatively suggested that, at least in the hottest low-gravity DA white dwarfs, the observed narrow-lined features are formed in expanding halos or winds associated with the white dwarfs. Theoretically, stable white dwarf halos should actually be coronae with temperatures in excess of 1,000,000 K. However, the observed narrow-lined features do not suggest such high temperatures. The observed radial velocities suggest weak stellar winds in two hot white dwarfs, namely, G191-B2B and 2111+49. It is tentatively proposed that radiative levitation can explain the appearance of the observed metallic lines in the hot DA white dwarfs.

Bruhweiler, F. C.↗

A frontal cyclogenesis case study from Seasat scatterometer data

In June 1978, the Seasat satellite was launched carrying, among other instruments, the Seasat-A scatterometer system (SASS), which produced ambiguous wind speed and direction data at the ocean surface. A fifteen day subset of dealiased wind vector data with the inherent ambiguities removed was produced for the period of September 6-20, 1978. On September 8, SASS began to observe a development of frontal cyclogenesis in the South Pacific off the east coast of New Zealand, in an area of few surface observations. A large mature cyclone contained weak warm and cold fronts and an occlusion with a strong horizontal wind shear. Satellite imagery shows that a strong upper-level jet streak was moving rapidly over the area of the surface frontal occlusion and as the jet passed over this area a new vortex formed. This cyclogenesis event was studied using 50-km resolution scatterometer surface wind data. High-resolution fields of wind vectors, divergence and vorticity are computed and plotted from the scatterometer data to study the structure and development of the newly formed cyclonic vortex, not otherwise possible using conventional observations.

Cunningham, Glenn F.↗

YBa2Cu3O(7-x)/Au/Nb device structures

Fabrication and testing of planar and edge geometry YBaCuO/Au/Nb superconductor/normal-metal/superconductor (SNS) device structures is described. Weak-link devices of this type serve as sensitive probes of the electrical quality of the YBaCuO/Au interface. The devices are fabricated using laser-ablated, in situ, c-axis-oriented YBaCuO thin films, with both annealed and unannealed YBaCuO/Au interfaces. The planar SNS structures are formed by sequential, in situ deposition of YBaCuO, Au, and Nb, followed by etching, planarization, and wiring electrode definition to produce junctions ranging from 5 to 20 micron on a side. Resulting RnA products are 1-10 x 10 to the -8th ohm-sq cm with critical current densities up to 5 kA/sq cm. For the edge geometry devices, the YBaCuO film edges are patterned using Ar ion milling, followed by a low energy ion cleaning step and in situ deposition of Au and Nb. Devices with areas in the 10 to the -7th to 10 to the -8th sq cm range have been fabricated with RnA products lower than 10 to the -8th ohm-sq cm and critical current densities up to 3kA/sq cm. Both types of devices show ac Josephson steps under microwave irradiation. The best results have been obtained with annealed YBaCuO/Au interfaces.

Hunt, B. D.↗