Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “different dimensions”

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 433 records · Page 24

Space-tiled colloidal crystals from DNA-forced shape-complementary polyhedra pairing

Generating space-filling arrangements of most discrete polyhedra nanostructures of the same shape is not possible. However, if the appropriate individual building blocks are selected (e.g., cubes), or multiple shapes of the appropriate dimensions are matched (e.g., octahedra and tetrahedra) and their pairing interactions are subsequently forced, space-filled architectures may be possible. With flexible molecular ligands (polyethylene glycol–modified DNA), the shape of a polyhedral nanoparticle can be deliberately altered and used to realize geometries that favor space tessellation. In this work, 10 new colloidal crystals were synthesized from DNA-modified nanocrystal building blocks that differed in shapes and sizes, designed to form space-filling architectures with micron-scale dimensions. The insights and capabilities provided by this new strategy substantially expand the scope of colloidal crystals possible and provide an expanded tool kit for researchers interested in designing metamaterials.

Science & Technology - Other Topics↗

Correlation of plastic events with local structure in jammed packings across spatial dimensions

Significance Mean-field theories, exact in the limit of infinite spatial dimensions, succeed in describing many features of glasses and amorphous solids in low dimensions, leading to considerable effort to understand how behavior evolves with dimension. Until now, all evidence has supported a picture in which “localized physics,” responsible for deviations from mean-field behavior in low dimensions, fades away with rising dimension. Our work shows that rearrangements, in which particles change relative positions leading to fluid-like response, reveal a different picture of dimensional cross-over. We find that rearrangements, which are localized in two- and three-dimensional systems and correlated with local structure, remain just as correlated with local structure up to five dimensions, suggesting that local structure is important even in high dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct E&M Field Solver for a Plane Wave without Sources

A three dimensional direct field solver for the electric and magnetic fields using the Finite Difference Time Domain (FDTD) method is described. Example simulations using periodic and absorbing boundary conditions in a vacuum medium for a plane wave propagating in three dimensions are explored and are compared to theory.

97 MATHEMATICS AND COMPUTING↗

High-precision quantum algorithms for partial differential equations

Quantum computers can produce a quantum encoding of the solution of a system of differential equations exponentially faster than a classical algorithm can produce an explicit description. However, while high-precision quantum algorithms for linear ordinary differential equations are well established, the best previous quantum algorithms for linear partial differential equations (PDEs) have complexity poly(1/ϵ), where ϵ is the error tolerance. By developing quantum algorithms based on adaptive-order finite difference methods and spectral methods, we improve the complexity of quantum algorithms for linear PDEs to be poly(d,log(1/ϵ)), where d is the spatial dimension. Our algorithms apply high-precision quantum linear system algorithms to systems whose condition numbers and approximation errors we bound. We develop a finite difference algorithm for the Poisson equation and a spectral algorithm for more general second-order elliptic equations.

97 MATHEMATICS AND COMPUTING↗

Energy-enhanced dimension eight SMEFT effects in VBF Higgs production

We study Higgs boson production via vector boson fusion at the LHC, focusing on the process pp → H + jj and capturing the leading energy-enhanced contributions within the Standard Model Effective Field Theory (SMEFT) up to order 1/Λ$^{4}$. Employing energy-scaling arguments, we predict the magnitude of each higher-dimensional operator’s contribution. Utilizing the geometric formulation of SMEFT, our analysis incorporates dimension-eight operators not previously considered. We find that the kinematics of vector boson fusion — characterized by two highly forward jets — tend to suppress contributions from higher-dimensional operators, requiring a lower scale Λ for SMEFT effects to become observable. This suggests that the SMEFT remains valid for lower Λ than expected. Combined with the fact that LEP constrains the dimension-six operators with the most considerable impact on vector boson fusion, a regime exists where dimension-eight operators can have significant effects. In many cases, these dimension-eight operators also influence associated production processes like pp → HV(jj), though differences in analysis cuts and kinematics mean this is not always the case. Our findings provide insights that could refine the search for SMEFT signals in collider experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effective entropy of quantum fields coupled with gravity

Entanglement entropy, or von Neumann entropy, quantifies the amount of uncertainty of a quantum state. For quantum fields in curved space, entanglement entropy of the quantum field theory degrees of freedom is well-defined for a fixed background geometry. In this paper, we propose a generalization of the quantum field theory entanglement entropy by including dynamical gravity. The generalized quantity named effective entropy, and its Renyi entropy generalizations, are defined by analytic continuation of a replica calculation. The replicated theory is defined as a gravitational path integral with multiple copies of the original boundary conditions, with a co-dimension-2 brane at the boundary of region we are studying. We discuss different approaches to define the region in a gauge invariant way, and show that the effective entropy satisfies the quantum extremal surface formula. When the quantum fields carry a significant amount of entanglement, the quantum extremal surface can have a topology transition, after which an entanglement island region appears. Our result generalizes the Hubeny-Rangamani-Takayanagi formula of holographic entropy (with quantum corrections) to general geometries without asymptotic AdS boundary, and provides a more solid framework for addressing problems such as the Page curve of evaporating black holes in asymptotic flat spacetime. We apply the formula to two example systems, a closed two-dimensional universe and a four-dimensional maximally extended Schwarzchild black hole. We discuss the analog of the effective entropy in random tensor network models, which provides more concrete understanding of quantum information properties in general dynamical geometries. We show that, in absence of a large boundary like in AdS space case, it is essential to introduce ancilla that couples to the original system, in order for correctly characterizing quantum states and correlation functions in the random tensor network. Using the superdensity operator formalism, we study the system with ancilla and show how quantum information in the entanglement island can be reconstructed in a state-dependent and observer-dependent map. We study the closed universe (without spatial boundary) case and discuss how it is related to open universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗

The Development of a Generalized Riser Flow Regime Map Based Upon Higher Moment and Chaotic Statistics Using Electrical Capacitance Volume Tomography (ECVT)

Dynamic analyses have been applied to the temporal signals from an Electro Capacitance Volume Tomography instrument located near mid-height on the riser of an industrial-scale cold-flow circulating fluidized bed to characterize gas-solids flow behavior in the riser. Twelve capacitance electrodes surround the cylindrical riser over a height of 1.3 m. The instrument used a neural network deconvolution algorithm to determine the spatially resolved solids fraction recorded at 52 Hz. Experiments were carried out over a range of gas and solids flows in the transport regime using a Geldart Group B bed material, high density polyethylene with mean particle size of 880 μm. The radial solids distribution was found to vary from one-time step to the next between profiles typical of laminar and turbulent flow. The duration of time spent in each of these flow profiles depended upon the operating regime – dilute, core-annular, or fast fluidized bed. The chaotic structure of the temporal data was characterized using the three conventional approaches: the first 4 moments from the distribution of signal in time, system memory parameters from the autocorrelation function and the Hurst exponent, and analysis of the correlationentropy and correlation dimension of the attractor. These signal analysis techniques were used to clearly distinguish differences between different transport operating regimes. Specifically, it was experimentally observed that a riser transitions from core annular flow profile to dilute and dense regimes via increasing the frequency of short term transients to either dilute or dense flow profiles, respectively. A regime map was generated based upon these dynamics using solids flux and gas velocity axes. Fast fluidized, core annular, and dilute each exhibited different degree of dynamic characteristics typical of fluid dominated or particle compromising behavior. It should be noted that the magnitude for the different statistics was in the same range regardless of the regime, it was the radial profile for the statistic that changed and subsequently identified that there was a change in the regime. Finally, a reduced regime map was developed consisting of plotting the gas velocity normalized by the upper transport velocity versus the solids flux normalized by the saturation carrying capacity. The use of this reduced plot allowed the data from widely different conditions to be plotted and compared on the same<p>graph. Note that in many instances, some of the statistics identified the operating point as being in one regime while others indicated that it was in another indicating a transition region between dilute or core annular regimes and between the core annular and fast fluidization regimes. This now provides a tool that can be used to optimize process performance, identify changes in operating states, or replicate process dynamics during process scaling or changing operating parameters. </p>

Breault, Ronald↗

Evanescent and inertial-like waves in rigidly rotating odd viscous liquids

Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves ( J. Fluid Mech ., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity 𝛺, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection ( J. Fluid Mech ., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers 𝜅 which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies 𝜔 < 2⁢𝛺. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.

Mechanics↗

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗

A comparison of artificial solar granules with real solar granules

The properties of computer-generated images of solar granules were compared with data from the literature and with observations of granules from 1975 and 1979. The lifetimes, shapes, and dimensions of the granules were estimated, and the results are discussed. No significant differences were found between the artificial images and the observed granules. The ratios of width to length among the artificial granules are given in a table.

Woehl, H.↗

A computational fluid dynamics algorithm on a massively parallel computer

The implementation and performance of a finite-difference algorithm for the compressible Navier-Stokes equations in two or three dimensions on the Connection Machine are described. This machine is a single-instruction multiple-data machine with up to 65536 physical processors. The implicit portion of the algorithm is of particular interest. Running times and megadrop rates are given for two- and three-dimensional problems. Included are comparisons with the standard codes on a Cray X-MP/48.

Jespersen, Dennis C.↗

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

Neurovestibular Session Summary

Oman - The early mission operational problems caused by space motion sickness have been largely resolved in recent years. This has been achieved by appropriate timeline adjustments, voluntary head movement restriction, and judicious use of promethazine. Crew members now simply accept that some symptoms "come with the job," and usually last only a few days. But as more people have flown longer flights, we've seen cases of space sickness and inversion illusion that take several weeks to resolve. Visual reorientation illusions continue throughout long flights, and occasionally cause difficulties. EVA astronauts sometimes suddenly fear they will fall out of the payload bay or off of the RMS or Strella arms. Orientation and navigation in three dimensions in the MIR station reportedly does not come naturally, because modules have different visual verticals. It is clear that the neurovestibular problems of spaceflight have not disappeared. After return to Earth, many crew members are disoriented and ataxic in the first hour after return, and require assistance leaving the vehicle, Flight surgeons say that the longer the mission, the stronger the aftereffects, certain of which last for weeks. We do not yet know how to predict who will be afflicted. Looking ahead to 3-4 month long voyages to Mars, it seems obvious that if cruise is in O-G, the crew may encounter neurovestibular problems on arrival. Artificial G may be broadly effective as a countermeasure for many of the physiological changes of spaceflight, but from the neurovestibular perspective, it is a double-edged sword. We know that the Coriolis stimulus resulting from rotation is potentially disorienting and nauseogenic. But we don't yet know how much artificial G will be enough, nor how successfully people can adapt to a specific angular velocity and hypo G level. Development of countermeasures remains a big challenge for our neurovestibular community. Maintaining an interdisciplinary perspective is important. Three examples were presented at this meeting: 1) Transgenic animal experiments suggest that in addition to the light illumination cycle, vestibular inputs may also serve as an important input to the circadian system. 2) Radiation can cause important CNS effects in animals, including loss of spatial memory. 3) As described in our session, otolith inputs may contribute to cardiovascular regulation of orthostatic tolerance. Over the past three days, we've all enjoyed catching up with old friends, and making many new ones. On behalf of my colleagues, I want to thank Al Coats and the USRA DSLS staff for the great job they did in running this meeting. And keeping the emphasis on fun. And also my Co- Chair, Mal Cohen, who had more stamina than many of us, despite major surgery only three weeks ago. Mal and I have written a few lines describing each of the seventeen papers in our session, to give you a quick over-view, and as a guide to the full abstracts, We have grouped them under five themes: preflight and inflight countermeasurements, postlanding posture and locomotion deficits: assessment and prediction, adaptive processes, relationships among physical simuli, perceptions, and eye movements, vestibular contribution to human autonomic responses, and implications and recommendations.

Oman, Charles↗

Steady Nuclear Combustion in Rockets

The astrophysical theory of stationary nuclear reactions in stars is applied to the conditions that would be met in the practical engineering cases that would differ from the former, particularly with respect to the much lower combustion pressures, dimensions of the reacting volume, and burnup times. This application yields maximum rates of hear production per unit volume of reacting gas occurring at about 10(exp 8) K in the cases of reactions between the hydrogen isotopes, but yields higher rates for heavier atoms. For the former, with chamber pressures of the order of 100 atmospheres, the energy production for nuclear combustion reaches values of about 10(exp 4) kilocalories per cubic meter per second, which approaches the magnitude for the familiar chemical fuels. The values are substantially lower for heavier atoms, and increase with the square of the combustion pressure. The half-life of the burnup in the fastest reactions may drop to values as low as those for chemical fuels so that, despite the high temperature, the radiated energy can remain smaller than the energy produced, particularly if an inefficiently radiating (i.e., easily completely ionized reacting material like hydrogen), is used. On the other hand, the fraction of completely ionized particles in the gases undergoing nuclear combustion must not exceed a certain upper limit because the densities (approximately 10(exp -10) grams per cubic centimeter)) lie in the range of high vacua and only for the previously mentioned fraction of nonionized particles can mean free paths be retained small enough so that the chamber diameters of several dozen meters will suffice. Under these conditions it appears that continuously maintained stable nuclear reactions at practical pressures and dimensions are fundamentally possible and their application can be visualized as energy sources for power plants and propulsion units.

Saenger, E.↗

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗

How gravity and joint scaling affect dynamics response

If a scale model of a space structure is ground tested to establish dynamic characteristics of full-scale space structure, two factors need to be taken into consideration: the joint tolerances may not be scaled, and joint preloads (from gravity) may be present. In some situations, these differences may significantly affect dynamic response so that the model is not representative of the structure in space. In this study, joint effects are investigated by considering the response of a single-degree-of-freedom system that is restrained by a joint linkage. The joint is characterized by its initial stiffness, the load at which it slips, its hardening rate, and the distance it can slip (joint 'slop' or gap). The steady-state response of the system is calculated using a harmonic balance approach. It is found that the system response is multivalued for certain ranges of gap dimensions. When this is the case, small changes in system properties can lead to large differences in dynamic response and make it difficult to obtain meaningful results directly from model testing.

Hsu, S. T.↗

Relative Navigation In Elliptical Orbits Using An Iterative Nonlinear Filter

The two step filter is applied to process intersatellite radar measurements to determine the motion of one satellite relative to another in close elliptical orbits. This filter breaks a nonlinear estimation problem into two state vectors. The 'first step' state is chosen so as to have a linear measurement equation. This is nonlinearly related to the 'second step' state which describes the dynamics. Two different forms are used. In one, the first step state is the second step state vector augmented by the measurement equation. In the other, the first step and second step state vectors are of equal dimension. The two step filter is compared against an iterated extended Kalman filter and a Kalman filter using a change of variables. Analytical differences between the two step estimator and these conventional filters are highlighted. Special concerns for initializing the first step state covariance matrix and handling the possibility of numerically rank deficient covariance matrices are addressed. Numerical simulations are performed which show that the Two Step estimator produces a lower estimation bias under two circumstances; large apriori initial error; and small dimension observation vectors which require a longer arc of measurements to generate observability of the state.

Garrison, James L.↗