Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “N CODES”

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 145 records · Page 8

A unified N-body and statistical treatment of stellar dynamics. III - Early postcollapse evolution of globular clusters

The period immediately following the core collapse phase in the evolution of a globular cluster is studied using a hybrid N-body/Fokker-Planck stellar dynamical code. Several core oscillations of the type predicted in earlier work are seen. The oscillations are driven by the formation, hardening, and ejection of binaries by three-body processes, and appear to decay on a timescale of about 10 to the 7th yr, for the choice of 'typical' cluster parameters made here. There is no evidence that they are gravothermal in nature. The mechanisms responsible for the decay are discussed in some detail. The distribution of hard binaries produced by the oscillations is compared with theoretical expectations and the longer term evolution of the system is considered.

Mcmillan, S. L. W.↗

Dynamical experiments on models of colliding disk galaxies

Collisions between galaxies can induce large morphological changes in the participants and, in the case of colliding disk galaxies, bridges and tails are often formed. Observations of such systems indicate a wide variation in color (see Larson and Tinsley, 1978) and that some of the particpants are experiencing enhanced rates of star formation, especially in their central regions (Bushouse 1986, 1987; Kennicutt et al., 1987, Bushouse, Lamb, and Werner, 1988). Here the authors describe progress made in understanding some of the dynamics of interacting galaxies using N-body stellar dynamical computer experiments, with the goal of extending these models to include a hydrodynamical treatment of the gas so that a better understanding of globally enhanced star formation will eventually be forthcoming. It was concluded that close interactions between galaxies can produce large perturbations in both density and velocity fields. The authors measured, via computational experiments that represent a galaxy's stars, average radial velocity flows as large as 100 km/sec and 400 percent density increases. These can occur in rings that move outwards through the disk of a galaxy, in roughly homologous inflows toward the nucleus, and in off center, non-axisymmetric regions. Here the authors illustrate where the gas is likely to flow during the early stages of interaction and in future work they plan to investigate the fate of the gas more realistically by using an N-body/Smoothed Particle Hydrodynamics code to model both the stellar and gaseous components of a disk galaxy during a collision. Specifically, they will determine the locations of enhanced gas density and the strength and location of shock fronts that form during the interaction.

Gerber, Richard A.↗

A mean field neural network for hierarchical module placement

This paper proposes a mean field neural network for the two-dimensional module placement problem. An efficient coding scheme with only O(N log N) neurons is employed where N is the number of modules. The neurons are evolved in groups of N in log N iteration steps such that the circuit is recursively partitioned in alternating vertical and horizontal directions. In our simulations, the network was able to find optimal solutions to all test problems with up to 128 modules.

Unaltuna, M. Kemal↗

Monte Carlo Analysis of Pion Contribution to Absorbed Dose from Galactic Cosmic Rays

Accurate knowledge of the physics of interaction, particle production and transport is necessary to estimate the radiation damage to equipment used on spacecraft and the biological effects of space radiation. For long duration astronaut missions, both on the International Space Station and the planned manned missions to Moon and Mars, the shielding strategy must include a comprehensive knowledge of the secondary radiation environment. The distribution of absorbed dose and dose equivalent is a function of the type, energy and population of these secondary products. Galactic cosmic rays (GCR) comprised of protons and heavier nuclei have energies from a few MeV per nucleon to the ZeV region, with the spectra reaching flux maxima in the hundreds of MeV range. Therefore, the MeV - GeV region is most important for space radiation. Coincidentally, the pion production energy threshold is about 280 MeV. The question naturally arises as to how important these particles are with respect to space radiation problems. The space radiation transport code, HZETRN (High charge (Z) and Energy TRaNsport), currently used by NASA, performs neutron, proton and heavy ion transport explicitly, but it does not take into account the production and transport of mesons, photons and leptons. In this paper, we present results from the Monte Carlo code MCNPX (Monte Carlo N-Particle eXtended), showing the effect of leptons and mesons when they are produced and transported in a GCR environment.

Aghara, S.K.↗

Using Neutron Spectroscopy to Constrain the Composition and Provenance of Phobos and Deimos

The origin of the Martian moons Phobos and Deimos is obscure and enigmatic. Hypotheses include the capture of asteroids originally from the outer main belt or beyond, residual material left over from Mars' formation, and accreted ejecta from a large impact on Mars, among others. Measurements of reflectance spectra indicate a similarity to dark, red D-type asteroids, but could indicate a highly space-weathered veneer. Here we suggest a way of constraining the near-surface composition of the two moons, for comparison to known meteoritic compositions. Neutron spectroscopy, particularly the thermal and epithermal neutron flux, distinguishes clearly between various classes of meteorites and varying hydrogen (water) abundances. Perhaps most surprising of all, a rendezvous with Phobos or Deimos is not necessary to achieve this. A low-cost mission based on the LADEE spacecraft design in an eccentric orbit around Mars can encounter Phobos every 2 weeks. As few as five flyby encounters at speeds of 2.3 kilometers per second and closest-approach distance of 3 kilometers provide sufficient data to distinguish between ordinary chondrite, water-bearing carbonaceous chondrite, ureilite, Mars surface, and aubrite compositions. A one-Earth year mission design includes many more flybys at lower speeds and closer approach distances, as well as similar multiple flybys at Deimos in the second mission phase, as described in the Phobos And Deimos Mars Environment (PADME) mission concept. This presentation will describe the expected thermal and epithermal neutron fluxes based on MCNP6 (Monte Carlo N (i.e. Neutron)-Particle transport code (version 6) simulations of different meteorite compositions and their uncertainties.

thtermal and epithermal neutron fluxes↗

Modeling Radiation Sources for Experimental Validation

To build on previously completed experimental work, I have been working to implement a model of our lab’s Americium-Beryllium neutron source and Ludlum detection instruments computationally to provide secondary verification of our experimental results. This has been done using the radiation transport code MCNP6 (Monte Carlo N-Particle). I have developed files to calculate the mass ratio of various composite materials that we have exposed to the radiation source, as well as the components of the detection instruments and converted all of these into geometric coordinates that define the experimental setup. Identifying the proper energy distribution and detector efficiencies as well as the correct interpretation of the detector physics have proven to be the most difficult and time-consuming challenges of the work. It is currently too early for results based on the most updated and accurate models. With this experience, the next step is to develop another source model of a proton beam to simulate exposures on other composite targets. This will be developed based on existing experimental datasets, but not used to validate them. Each of these projects will further define the radiation shielding needs for humans in space radiation environments that astronauts will someday be exposed to beyond Low Earth Orbit.

Juliana Simon↗

A finite element conjugate gradient FFT method for scattering

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

Collins, Jeffery D.↗

Fault-Tolerant Coding for State Machines

Two reliable fault-tolerant coding schemes have been proposed for state machines that are used in field-programmable gate arrays and application-specific integrated circuits to implement sequential logic functions. The schemes apply to strings of bits in state registers, which are typically implemented in practice as assemblies of flip-flop circuits. If a single-event upset (SEU, a radiation-induced change in the bit in one flip-flop) occurs in a state register, the state machine that contains the register could go into an erroneous state or could hang, by which is meant that the machine could remain in undefined states indefinitely. The proposed fault-tolerant coding schemes are intended to prevent the state machine from going into an erroneous or hang state when an SEU occurs. To ensure reliability of the state machine, the coding scheme for bits in the state register must satisfy the following criteria: 1. All possible states are defined. 2. An SEU brings the state machine to a known state. 3. There is no possibility of a hang state. 4. No false state is entered. 5. An SEU exerts no effect on the state machine. Fault-tolerant coding schemes that have been commonly used include binary encoding and "one-hot" encoding. Binary encoding is the simplest state machine encoding and satisfies criteria 1 through 3 if all possible states are defined. Binary encoding is a binary count of the state machine number in sequence; the table represents an eight-state example. In one-hot encoding, N bits are used to represent N states: All except one of the bits in a string are 0, and the position of the 1 in the string represents the state. With proper circuit design, one-hot encoding can satisfy criteria 1 through 4. Unfortunately, the requirement to use N bits to represent N states makes one-hot coding inefficient.

Naegle, Stephanie Taft↗

Radiation Transport Tools for Space Applications: A Review

This slide presentation contains a brief discussion of nuclear transport codes widely used in the space radiation community for shielding and scientific analyses. Seven radiation transport codes that are addressed. The two general methods (i.e., Monte Carlo Method, and the Deterministic Method) are briefly reviewed.

Cosmic Ray Effects on Micro-Electronics (CREME96)↗

Computer simulation of plasma and N-body problems

The following FORTRAN language computer codes are presented: (1) efficient two- and three-dimensional central force potential solvers; (2) a three-dimensional simulator of an isolated galaxy which incorporates the potential solver; (3) a two-dimensional particle-in-cell simulator of the Jeans instability in an infinite self-gravitating compressible gas; and (4) a two-dimensional particle-in-cell simulator of a rotating self-gravitating compressible gaseous system of which rectangular coordinate and superior polar coordinate versions were written.

Harries, W. L.↗

N-Factor Computations for the X-33 Vehicle

N-Factor computations for the X-33 vehicle on the lower side of the body were per-formed. The transition prediction code e-Malik was used to compute the N-Factors and the meanflow field was computed using the CFD code LAURA. The computations are done at two angles of attack 30 and 40 degrees and the freestream Mach number and Reynolds number are M= 6 and Re = 7.7* 10(exp 6). The N-Factors obtained are very low in the range of 2-3 near the experimentally observed transition onset regions. The stability computations show that at 30 degrees angle of attack, the transition occurs due to crossflow instability and at higher angles of attack the transition occurs due to first mode inviscid instability.

Balakumar, P.↗

Codes with Parity Conditions on Subsets of Coordinates

Binary codes with the constraint that the codes restricted to certain subsets of columns must be contained in particular codes of the shorter lengths are considered. In particular, codes of even length 2k, and of minimum distance approximately greater than d, where in the code obtained by restricting to the first k positions has even weight and at the same time the code obtained by restricting to the last k positions also has even weight are considered. If k = 2n, n odd, and d = 2n, it is proved that the code has at most 8n - 4 codewords, and 8n - 4 is attainable for n = 3. This permits a file-transfer protocol control function assignment for personal computers to be chosen for 20 control functions using essentially just pairs of upper-case alphabetic ASCII characters where the Hamming distance between the binary forms of every two different control functions is at least six.

Posner, E. C.↗

Hydrodynamic solutions at a sonic-boom focus

Study of the focusing of sonic-boom N waves using a numerical shock-following method. The computer code used in the calculations employs a moving mesh so that weak N wave shocks can be followed for large distances without numerical dissipation. In several sample calculations, overpressures at the focus ranging from 2 to 20 times the nominal overpressure have been computed. It is also found that weak shocks tend to cusp at the axis in accord with the geometric-acoustics description, whereas strong shocks straighten out without cusping.

Parker, L. W.↗

Large CMOS imager using hadamard transform based multiplexing

We have developed a concept design for a large (~10k x 10k) CMOS imaging array whose elements are grouped in small subarrays with N pixels in each. The subarrays are code-division multiplexed using the Hadamard Transform (HT) based encoding. The Hadamard code improves the signal-to-noise (SNR) ratio to the reference of the read-out amplifier by a factor of N^1/2. This way of grouping pixels reduces the number of hybridization bumps by N. A single chip layout has been designed and the architecture of the imager has been developed to accommodate the HT base multiplexing into the existing CMOS technology. The imager architecture allows for a trade-off between the speed and the sensitivity. The envisioned imager would operate at a speed >100 fps with the pixel noise < 20 e-. The power dissipation would be ~100 pW/pixe1. The combination of the large format, high speed, high sensitivity and low power dissipation can be very attractive for space reconnaissance applications.

CMOS imager↗

Accelerated panel methods using the fast multipole method

Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Leathrum, James F., Jr.↗

Interpolation in numerical optimization

The present work discusses the generation of the cubic-spline interpolator in numerical optimization methods which use a variable-step integrator with step size control based on local relative truncation error. An algorithm for generating the cubic spline with successive over-relaxation is presented which represents an improvement over that given by Ralston and Wilf (1967). Rewriting the code reduces the number of N-vectors from eight to one. The algorithm is formulated in such a way that the solution of the linear system set up yields the first derivatives at the nodal points. This method is as accurate as other schemes but requires the minimum amount of storage.

Hall, K. R.↗