Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 721 records · Page 40

Rotational symmetry protected edge and corner states in Abelian topological phases

Spatial symmetries can enrich the topological classification of interacting quantum matter and endow systems with non-trivial strong topological invariants (protected by internal symmetries) with additional "weak" topological indices. In this paper, we study the edge physics of systems with a non-trivial shift invariant, which is protected by either a continuous $\text{U}(1)_r$ or discrete $\text{C}_n$ rotation symmetry, along with internal $\text{U}(1)_c$ charge conservation. Specifically, we construct an interface between two systems which have the same Chern number but are distinguished by their Wen-Zee shift and, through analytic arguments supported by numerics, show that the interface hosts counter-propagating gapless edge modes which cannot be gapped by arbitrary local symmetry-preserving perturbations. Using the Chern-Simons field theory description of two-dimensional Abelian topological orders, we then prove sufficient conditions for continuous rotation symmetry protected gapless edge states using two complementary approaches. One relies on the algebraic Lagrangian sub-algebra framework for gapped boundaries while the other uses a more physical flux insertion argument. For the case of discrete rotation symmetries, we extend the field theory approach to show the presence of fractional corner charges for Abelian topological orders with gappable edges, and compute them in the case where the Abelian topological order is placed on the two-dimensional surface of a Platonic solid. Finally, our work paves the way for studying the edge physics associated with spatial symmetries in symmetry enriched topological phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗

Automation of measurement of spatial damping parameters

A traveling wave unit was built to process measurement data on spatial damping parameters and to simplify statistical treatment of the results. Mean values of the loss coefficient, lengths of flexual waves, and wave propagation velocities for a given set of discrete functions were determined. Results are given in a table.

Akolzin, A. N.↗

A space-time discretization procedure for wave propagation problems

Higher order compact algorithms are developed for the numerical simulation of wave propagation by using the concept of a discrete dispersion relation. The dispersion relation is the imprint of any linear operator in space-time. The discrete dispersion relation is derived from the continuous dispersion relation by examining the process by which locally plane waves propagate through a chosen grid. The exponential structure of the discrete dispersion relation suggests an efficient splitting of convective and diffusive terms for dissipative waves. Fourth- and eighth-order convection schemes are examined that involve only three or five spatial grid points. These algorithms are subject to the same restrictions that govern the use of dispersion relations in the constructions of asymptotic expansions to nonlinear evolution equations. A new eighth-order scheme is developed that is exact for Courant numbers of 1, 2, 3, and 4. Examples are given of a pulse and step wave with a small amount of physical diffusion.

Davis, Sanford↗

Implementing system simulation of C3 systems using autonomous objects

The basis of all conflict recognition in simulation is a common frame of reference. Synchronous discrete-event simulation relies on the fixed points in time as the basic frame of reference. Asynchronous discrete-event simulation relies on fixed-points in the model space as the basic frame of reference. Neither approach provides sufficient support for autonomous objects. The use of a spatial template as a frame of reference is proposed to address these insufficiencies. The concept of a spatial template is defined and an implementation approach offered. Discussed are the uses of this approach to analyze the integration of sensor data associated with Command, Control, and Communication systems.

Rogers, Ralph V.↗

An optimized finite-difference scheme for wave propagation problems

Two fully-discrete finite-difference schemes for wave propagation problems are presented, a maximum-order scheme and an optimized (or spectral-like) scheme. Both combine a seven-point spatial operator and an explicit six-stage time-march method. The maximum-order operator is fifth-order in space and is sixth-order in time for a linear problem with periodic boundary conditions. The phase and amplitude errors of the schemes obtained using Fourier analysis are given and compared with a second-order and a fourth-order method. Numerical experiments are presented which demonstrate the usefulness of the schemes for a range of problems. For some problems, the optimized scheme leads to a reduction in global error compared to the maximum-order scheme with no additional computational expense.

Zingg, D. W.↗

Semi-discrete Galerkin solution of the compressible boundary-layer equations with viscous-inviscid interaction

A semi-discrete Galerkin (SDG) method is under development to model attached, turbulent, and compressible boundary layers for transonic airfoil analysis problems. For the boundary-layer formulation the method models the spatial variable normal to the surface with linear finite elements and the time-like variable with finite differences. A Dorodnitsyn transformed system of equations is used to bound the infinite spatial domain thereby providing high resolution near the wall and permitting the use of a uniform finite element grid which automatically follows boundary-layer growth. The second-order accurate Crank-Nicholson scheme is applied along with a linearization method to take advantage of the parabolic nature of the boundary-layer equations and generate a non-iterative marching routine. The SDG code can be applied to any smoothly-connected airfoil shape without modification and can be coupled to any inviscid flow solver. In this analysis, a direct viscous-inviscid interaction is accomplished between the Euler and boundary-layer codes through the application of a transpiration velocity boundary condition. Results are presented for compressible turbulent flow past RAE 2822 and NACA 0012 airfoils at various freestream Mach numbers, Reynolds numbers, and angles of attack.

Day, Brad A.↗

Aeroacoustic computation of a gust-plate interaction via MacCormack schemes

We apply second, fourth, and sixth order spatially accurate variations of the MacCormack scheme to calculate the noise radiated by the interaction of a flat plate with an oncoming gust. For the given gust wavelength and numerical discretization, the fourth and sixth order schemes each are effective in capturing the occurring acoustic waves. However, for the case of the sixth order scheme, the high order of extrapolation applied at the far field boundaries necessitates an extension of the computational boundaries from their prescribed location.

Martin, James E.↗

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

Clouds on Neptune: Motions, Evolution, and Structure

The aims of our original proposal were these: (1) improving measurements of Neptune's circulation, (2) understanding the spatial distribution of cloud features, (3) discovery of new cloud features and understanding their evolutionary process, (4) understanding the vertical structure of zonal cloud patterns, (5) defining the structure of discrete cloud features, and (6) defining the near IR albedo and light curve of Triton. Towards these aims we proposed analysis of existing 1996 groundbased NSFCAM/IRTF observations and nearly simultaneous WFPC2 observations from the Hubble Space Telescope. We also proposed to acquire new observations from both HST and the IRTF.

Sromovsky, Larry A.↗

Computational Methods for Jet Noise Simulation

The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Goodrich, John W.↗

Low-Dissipation Advection Schemes Designed for Large Eddy Simulations of Hypersonic Propulsion Systems

The 2nd-order upwind inviscid flux scheme implemented in the multi-block, structured grid, cell centered, finite volume, high-speed reacting flow code VULCAN has been modified to reduce numerical dissipation. This modification was motivated by the desire to improve the codes ability to perform large eddy simulations. The reduction in dissipation was accomplished through a hybridization of non-dissipative and dissipative discontinuity-capturing advection schemes that reduces numerical dissipation while maintaining the ability to capture shocks. A methodology for constructing hybrid-advection schemes that blends nondissipative fluxes consisting of linear combinations of divergence and product rule forms discretized using 4th-order symmetric operators, with dissipative, 3rd or 4th-order reconstruction based upwind flux schemes was developed and implemented. A series of benchmark problems with increasing spatial and fluid dynamical complexity were utilized to examine the ability of the candidate schemes to resolve and propagate structures typical of turbulent flow, their discontinuity capturing capability and their robustness. A realistic geometry typical of a high-speed propulsion system flowpath was computed using the most promising of the examined schemes and was compared with available experimental data to demonstrate simulation fidelity.

White, Jeffrey A.↗

Beam combining of high-power terahertz lasers with semiconductor metasurface gratings

Spectral beam combining (SBC) of an array of single-mode lasers has the potential for realizing a chip-scale electrically tunable multispectral laser source without any moving parts and with emission at closely spaced discrete frequencies. Here, we demonstrate SBC of an array of single-mode THz quantum-cascade lasers (QCLs) with small frequency spacing. Semiconductor-based metasurface THz blazed gratings are developed for integration with QCLs in a compact cryocooler setup. Single-lobed beams of four QCLs emitting around 3.2 THz and spectrally separated by ~14 GHz are combined to achieve collinear propagation of spatially overlapped and collimated beams in the far field with ≲ 0.1° of pointing error. Peak optical power in the range 50–100 mW is detected from each of the QCLs after combining. The dense frequency spacing in our method provides a pathway to incorporate tens of QCLs on a chip to realize a high-power discretely tunable THz laser source for applications in THz spectroscopy.

Jia, Fei [Lehigh Univ., Bethlehem, PA (United Stat↗

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.

Lemieux, Jean-François (ORCID:0000000320845759)↗

Spatial modeling algorithms for reactions and transport in biological cells

Biological cells rely on precise spatiotemporal coordination of biochemical reactions to control their functions. Such cell signaling networks have been a common focus for mathematical models, but they remain challenging to simulate, particularly in realistic cell geometries. Here we present Spatial Modeling Algorithms for Reactions and Transport (SMART), a software package that takes in high-level user specifications about cell signaling networks and then assembles and solves the associated mathematical systems. SMART uses state-of-the-art finite element analysis, via the FEniCS Project software, to efficiently and accurately resolve cell signaling events over discretized cellular and subcellular geometries. We demonstrate its application to several different biological systems, including yes-associated protein (YAP)/PDZ-binding motif (TAZ) mechanotransduction, calcium signaling in neurons and cardiomyocytes, and ATP generation in mitochondria. Throughout, we utilize experimentally derived realistic cellular geometries represented by well-conditioned tetrahedral meshes. These scenarios demonstrate the applicability, flexibility, accuracy and efficiency of SMART across a range of temporal and spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗

The microstrip proportional counter

Microstrip detectors in which the usual discrete anode and cathode wires are replaced by conducting strips on an insulating or partially insulating substrate are fabricated using integrated circuit-type photolithographic techniques and hence offer very high spatial accuracy and uniformity, together with the capability of producing extremely fine electrode structures. Microstrip proportional counters have now been variously reported having an energy resolution of better than 11 percent FWHM at 5.9 keV. They have been fabricated with anode bars down to 2 microns and on a variety of substrate materials including thin films which can be molded to different shapes. This review will examine the development of the microstrip detector with emphasis on the qualities which make this detector particularly interesting for use in astronomy.

Ramsey, B. D.↗

Development of coronal mass ejections and association with interplanetary events

Results are presented on the development of two coronal mass ejections (CMEs) obtained by comparing the observations of the large angle spectroscopic coronagraph (LASCO) and the extreme ultraviolet imaging telescope (EIT) instrument onboard the SOHO with those of the Nancay radioheliograph. The radioheliograph provides images at five levels in the corona. An excellent spatial association is found between the position and extent of the type 4 radio sources and the CMEs seen by LASCO. One result is the existence for these two events of discrete successive phases in their development. For these events, Ulysses and SOHO missions measured interplanetary particles of coronal origin. The coronal acceleration site was attempted to be identified, as well as the path of these particles from the corona to the interplanetary medium.

Pick, M.↗