Sample records for one-dimensional compass model

In order to investigate the quantum phase transition in the one-dimensional quantum compassmodel, we numerically calculate non-local string correlations, entanglement entropy and fidelity per lattice site by using the infinite matrix product state representation with the infinite time evolving block decimation method. In the whole range of the interaction parameters, we find that four distinct string orders characterize the four different Haldane phases and the topological quantum phase transition occurs between the Haldane phases. The critical exponents of the string order parameters β = 1/8 and the cental charges c = 1/2 at the critical points show that the topological phase transitions between the phases belong to an Ising type of universality classes. In addition to the string order parameters, the singularities of the second derivative of the ground state energies per site, the continuous and singular behaviors of the Von Neumann entropy and the pinch points of the fidelity per lattice site manifest that the phase transitions between the phases are of the second-order, in contrast to the first-order transition suggested in previous studies.

The onedimensionalmodel of neutron dynamic in reactor core was developed. The core was divided in several axial nodes. The one group neutron diffusion equation for each node is solved. Feedback affects of fuel and water temperatures is calculated. The influence of xenon, boron and control rods is included in cross section calculations for each node. The system of equations is solved implicitly. The model is used in basic principle Training Simulator of NPP Krsko. (author)

The general expression for the dispersion relation for a polyatomic onedimensional crystal obtained by the Laplace Transform Method is applied to materials with the fcc and hcp structures, both consisting of close-packed planes of atoms with the stacking sequence of plane ABC/ABC... and AB/AB... respectively. The expression is also applied to polytypes, that is materials caracterized by a stacking sequence with longer repeat unit. The effective mass is cast in a condensed form useful for further calculations. The results from this simple model are only qualitative. (Author) [pt

A one-dimensionalmodel is introduced in this paper for problems of internal ballistics involving solid propellant combustion. First, the work presents the physical approach and equations adopted. Closure relationships accounting for the physical phenomena taking place during combustion (interfacial friction, interfacial heat transfer, combustion) are deeply discussed. Secondly, the numerical method proposed is presented. Finally, numerical results provided by this code (UXGun) are compared with results of experimental tests and with the outcome from a well-known zero-dimensional code. The model provides successful results in firing tests of artillery guns, predicting with good accuracy the maximum pressure in the chamber and muzzle velocity what highlights its capabilities as prediction/design tool for internal ballistics.

The subject of transport near the separatrix in FRC devices is important for determining the performance to be expected from an FRC reactor or from FRC experiments. A computer code was constructed for studying the micro-stability properties of FRCs near the separatrix as a first step in obtaining quasilinear transport coefficients that can be used in a transport code. We consider collisionless ions and electrons, without an expansion in powers of a parameter, like the electron or ion gyroradius, and we approximate the equilibrium with an infinitely long axially and translationally symmetric equilibrium. Thus, in our equilibria, there are only an axial magnetic field and a radial electric field. Our equilibria are collisionless, two-species, diffuse-profile, one-dimensional, theta-pinch equilibria. We allow the possibility that there be a magnetic field null in order to be able to model FRC devices more realistically

A one-dimensionalmodel of inertial pumping is introduced and solved. The pump is driven by a high-pressure vapor bubble generated by a microheater positioned asymmetrically in a microchannel. The bubble is approximated as a short-term impulse delivered to the two fluidic columns inside the channel. Fluid dynamics is described by a Newton-like equation with a variable mass, but without the mass derivative term. Because of smaller inertia, the short column refills the channel faster and accumulates a larger mechanical momentum. After bubble collapse the total fluid momentum is nonzero, resulting in a net flow. Two different versions of the model are analyzed in detail, analytically and numerically. In the symmetrical model, the pressure at the channel-reservoir connection plane is assumed constant, whereas in the asymmetrical model it is reduced by a Bernoulli term. For low and intermediate vapor bubble pressures, both models predict the existence of an optimal microheater location. The predicted net flow in the asymmetrical model is smaller by a factor of about 2. For unphysically large vapor pressures, the asymmetrical model predicts saturation of the effect, while in the symmetrical model net flow increases indefinitely. Pumping is reduced by nonzero viscosity, but to a different degree depending on the microheater location.

RETRAN-02 is a modular code system that has been designed for one-dimensional, transient thermal-hydraulics analysis. In RETRAN-02, core power behavior may be treated using a one-dimensional reactor kinetics model. This model allows the user to investigate the interaction of time- and space-dependent effects in the reactor core on overall system behavior for specific LWR operational transients. The purpose of this paper is to review the recent analysis and development activities related to the onedimensional kinetics model in RETRAN-02

This paper describes a one-dimensional spatial neutron kinetics model that was developed for the RETRAN code. The RETRAN -01 code has a point kinetics model to describe the reactor core behavior during thermal-hydraulic transients. A one-dimensional neutronics model has been developed for RETRAN-02. The ability to account for flux shape changes will permit an improved representation of the thermal and hydraulic feedback effects for many operational transients. 19 refs

Excitons in carbon nanotubes may be modeled by two oppositely charged particles living on the surface of a cylinder. We derive three onedimensional effective Hamiltonians which become exact as the radius of the cylinder vanishes. Two of them are solvable.......Excitons in carbon nanotubes may be modeled by two oppositely charged particles living on the surface of a cylinder. We derive three onedimensional effective Hamiltonians which become exact as the radius of the cylinder vanishes. Two of them are solvable....

Excitons in carbon nanotubes may be modeled by two oppositely charged particles living on the surface of a cylinder. We derive three one-dimensional effective Hamiltonians which become exact as the radius of the cylinder vanishes. Two of them are solvable.......Excitons in carbon nanotubes may be modeled by two oppositely charged particles living on the surface of a cylinder. We derive three one-dimensional effective Hamiltonians which become exact as the radius of the cylinder vanishes. Two of them are solvable....

Full Text Available Stirling engines are regarded as an efficient and promising power system for underwater devices. Currently, many researches on one-dimensionalmodel is used to evaluate thermodynamic performance of Stirling engine, but in which there are still some aspects which cannot be modeled with proper mathematical models such as mechanical loss or auxiliary power. In this paper, a four-cylinder double-acting Stirling engine for Unmanned Underwater Vehicles (UUVs is discussed. And a one-dimensionalmodel incorporated with empirical equations of mechanical loss and auxiliary power obtained from experiments is derived while referring to the Stirling engine computer model of National Aeronautics and Space Administration (NASA. The P-40 Stirling engine with sufficient testing results from NASA is utilized to validate the accuracy of this one-dimensionalmodel. It shows that the maximum error of output power of theoretical analysis results is less than 18% over testing results, and the maximum error of input power is no more than 9%. Finally, a Stirling engine for UUVs is designed with Schmidt analysis method and the modified one-dimensionalmodel, and the results indicate this designed engine is capable of showing desired output power.

We study quasi-one-dimensional scattering of one and two particles with short-range interactions on a discrete lattice model in two dimensions. One of the directions is tightly confined by an arbitrary trapping potential. We obtain the collisional properties of these systems both at finite and zero...

This paper presents an approximate inelastic analysis based on geometric simplification with emphasis on its applicability, modeling, and the method of defining the loading conditions. Two problems are investigated: a one-dimensional axisymmetric model of generalized plane strain thick-walled cylinder is applied to the primary sodium inlet nozzle of the Clinch River Breeder Reactor Intermediate Heat Exchanger (CRBRP-IHX), and a finite cylindrical shell is used to simulate the branch shell forging (Y) junction. The results are then compared with the available detailed inelastic analyses under cyclic loading conditions in terms of creep and fatigue damages and inelastic ratchetting strains per the ASME Code Case N-47 requirements. In both problems, the one-dimensional simulation is able to trace the detailed stress-strain response. The quantitative comparison is good for the nozzle, but less satisfactory for the Y junction. Refinements are suggested to further improve the simulation

The SAIL model for one-dimensional homogeneous vegetation canopies has been modified to include the specular reflectance and hot spot effects. This modified model and the Nilson-Kuusk model are evaluated by comparing the reflectances given by them against those given by a radiosity-based computer model, Diana, for a set of canopies, characterized by different leaf area index (LAI) and leaf angle distribution (LAD). It is shown that for homogeneous canopies, the analytical models are generally quite accurate in the visible region, but not in the infrared region. For architecturally realistic heterogeneous canopies of the type found in nature, these models fall short. These shortcomings are quantified.

A Green's function formalism is used to calculate the energy of impurity modes associated with one and/or two magnetic impurities in the one-dimensional Heisenberg XXZ magnetic chain. The system can be tuned from the Heisenberg to the Ising model varying a parameter λ. A numerical study is performed showing two types of localized modes (s and p). The modes depend on λ and the degeneracy of the acoustic modes is broken.

Turbulent cascade processes are studied in terms of a one-dimensional forest-fire model. A hier- archy of steady-state equations for the forests and the holes between them is constructed and solved within a mean-field closure scheme. The exact hole distribution function is found to be N H (s)=4N/[s(s+1)(s+2)], where N is the number of forests

This paper presents a one-dimensional energy flow model to investigate the energy behavior for poroelastic media coupled with acoustical media. The proposed energy flow model is expressed by an independent energy governing equation that is classified into each wave component propagating in poroelastic media. The energy governing equation is derived using the General Energetic Method (GEM). To facilitate a comparison with the classical solution based on the conventional displacement-base formulation, approximate solutions of energy density and intensity are obtained. Furthermore, the limitations and usability of the proposed energy flow model for poroelastic media are described.

In this article, we present a computational modeling, which gives us a dynamic view of some applications of Nuclear Engineering, specifically in the power distribution and the effective multiplication factor (keff) calculations. We work with one-dimensional problems of deterministic neutron transport theory, with the linearized Boltzmann equation in the discrete ordinates (SN) formulation, independent of time, with isotropic scattering and then built a software (Simulator) for modeling computational problems used in a typical calculations. The program used in the implementation of the simulator was Matlab, version 7.0. (author)

We study the one-dimensional behavior of a cellular automaton aimed at the description of the formation and evolution of cultural domains. The model exhibits a non-equilibrium transition between a phase with all the system sharing the same culture and a disordered phase of coexisting regions with different cultural features. Depending on the initial distribution of the disorder the transition occurs at different values of the model parameters. This phenomenology is qualitatively captured by a mean-field approach, which maps the dynamics into a multi-species reaction-diffusion problem.

Field II is a program for simulating ultrasound transducer fields. It is capable of calculating the emitted and pulse-echoed fields for both pulsed and continuous wave transducers. To make it fully calibrated a model of the transducer’s electro-mechanical impulse response must be included. We...... examine an adapted onedimensional transducer model originally proposed by Willatzen [9] to calibrate Field II. This model is modified to calculate the required impulse responses needed by Field II for a calibrated field pressure and external circuit current calculation. The testing has been performed...... to the calibrated Field II program for 1, 4, and 10 cycle excitations. Two parameter sets were applied for modeling, one real valued Pz27 parameter set, manufacturer supplied, and one complex valued parameter set found in literature, Alguer´o et al. [11]. The latter implicitly accounts for attenuation. Results show...

Previous versions of RETRAN have had only a point kinetics model to describe the reactor core behavior during thermal-hydraulic transients. The principal assumption in deriving the point kinetics model is that the neutron flux may be separated into a time-dependent amplitude funtion and a time-independent shape function. Certain types of transients cannot be correctly analyzed under this assumption, since proper definitions for core average quantities such as reactivity or lifetime include the inner product of the adjoint flux with the perturbed flux. A one-dimensional neutronics model has been included in a preliminary version of RETRAN-02. The ability to account for flux shape changes will permit an improved representation of the thermal and hydraulic feedback effects. This paper describes the neutronics model and discusses some of the analyses

In this work, we use a model developed to deduce a one-dimensionalmodel for the description of the poling of ferroelectric ceramics. This is built within the scheme of the thermodynamical theory of internal variables. The model produces both plastic and electric hysteresis effects in the form of ''plasticity'', i.e., rate-independent evolution equations for the plastic strain, and the residual electric polarization and both mechanical and electric hardenings. The influence of stresses on ferroelectric hysteresis loops through piezoelectricity and electrostriction is a natural outcome of this model. Some simple experimental methods for the determination of the material coefficients of the considered ceramics are suggested. (author). 21 refs, 3 figs

We study the spin-1/2 Ising chain with multispin interactions K involving the product of m successive spins, for general values of m. Using a change of spin variables the zero-field partition function of a finite chain is obtained for free and periodic boundary conditions and we calculate the two-spin correlation function. When placed in an external field H the system is shown to be self-dual. Using another change of spin variables the one-dimensional Ising model with multispin interactions in a field is mapped onto a zero-field rectangular Ising model with first-neighbour interactions K and H. The 2D system, with size m × N/m, has the topology of a cylinder with helical BC. In the thermodynamic limit N/m\\to ∞ , m\\to ∞ , a 2D critical singularity develops on the self-duality line, \\sinh 2K\\sinh 2H=1.

We propose feasible scenarios for revealing the modified exchange statistics in one-dimensional anyon models in optical lattices based on an extension of the multicolor lattice-depth modulation scheme introduced in [Phys. Rev. A 94, 023615 (2016), 10.1103/PhysRevA.94.023615]. We show that the fast modulation of a two-component fermionic lattice gas in the presence a magnetic field gradient, in combination with additional resonant microwave fields, allows for the quantum simulation of hardcore anyon models with periodic boundary conditions. Such a semisynthetic ring setup allows for realizing an interferometric arrangement sensitive to the anyonic statistics. Moreover, we show as well that simple expansion experiments may reveal the formation of anomalously bound pairs resulting from the anyonic exchange.

As a joint is loaded, the tangent stiffness of the joint reduces due to slip at interfaces. This stiffness reduction continues until the direction of the applied load is reversed or the total interface slips. Total interface slippage in joints is called macro-slip. For joints not undergoing macro-slip, when load reversal occurs the tangent stiffness immediately rebounds to its maximum value. This occurs due to stiction effects at the interface. Thus, for periodic loads, a softening and rebound hardening cycle is produced which defines a hysteretic, energy absorbing trajectory. For many jointed sub-structures, this hysteretic trajectory can be approximated using simple polynomial representations. This allows for complex joint substructures to be represented using simple non-linear models. In this paper a simple onedimensionalmodel is discussed.

Full Text Available Studies of the power-law relations of seismicity and earthquake source parameters based on the one-dimensional (1-D Burridge-Knopoff¡¦s (BK dynamical lattice model, especially those studies conducted by Taiwan¡¦s scientists, are reviewed in this article. In general, velocity- and/or state-dependent friction is considered to control faulting. A uniform distribution of breaking strengths (i.e., the static friction strength is taken into account in some studies, and inhomogeneous distributions in others. The scaling relations in these studies include: Omori¡¦s law, the magnitude-frequency or energy-frequency relation, the relation between source duration time and seismic moment, the relation between rupture length and seismic moment, the frequency-length relation, and the source power spectra. The main parameters of the one-dimensional (1-D Burridge-Knopoff¡¦s (BK dynamical lattice model include: the decreasing rate (r of dynamic friction strength with sliding velocity; the type and degree of heterogeneous distribution of the breaking strengths, the stiffness ratio (i.e., the ratio between the stiffness of the coil spring connecting two mass elements and that of the leaf spring linking a mass element and the moving plate; the frictional drop ratio of the minimum dynamic friction strength to the breaking strength; and the maximum breaking strength. For some authors, the distribution of the breaking strengths was considered to be a fractal function. Hence, the fractal dimension of such a distribution is also a significant parameter. Comparison between observed scaling laws and simulation results shows that the 1-D BK dynamical lattice model acceptably approaches fault dynamics.

In 1982, one of the authors (Okazaki), of Toshiba Corporation, wrote a one-dimensional, rigid-disk model computer program to serve as a reliable design tool for the 150 MW klystron development project. This is an introductory note for the users of this program. While reviewing the so-called disk programs presently available, hypotheses such as gridded interaction gaps, a linear relation between phase and position, and so on, were found. These hypotheses bring serious limitations and uncertainties into the computational results. JPNDISK was developed to eliminate these defects, to follow the equations of motion as rigorously as possible, and to obtain self-consistent solutions for the gap voltages and the electron motion. Although some inaccuracy may be present in the relativistic region, JPNDISK, in its present form, seems a most suitable tool for klystron design; it is both easy and inexpensive to use

The quantum mechanical description of the optical properties of crystalline materials typically requires extensive numerical computation. Including excitonic and non-perturbative field effects adds to the complexity. In one dimension, however, the analysis simplifies and optical spectra can be computed exactly. In this paper, we apply the Wannier exciton formalism to derive analytical expressions for the optical response in four cases of increasing complexity. Thus, we start from free carriers and, in turn, switch on electrostatic fields and electron–hole attraction and, finally, analyze the combined influence of these effects. In addition, the optical response of impurity-localized excitons is discussed. - Highlights: • Optical response of one-dimensional semiconductors including excitons. • Analytical model of excitonic Franz–Keldysh effect. • Computation of optical response of impurity-localized excitons

method to obtain the zero-temperature phase diagram of the one-dimensional, extended ... Progress in this field has been driven by an interplay between ... superconductor-insulator transition in thin films of superconducting materials like bis-.

Full Text Available In modern information technology, as integration density increases rapidly and the dimension of materials reduces to nanoscale, interfacial thermal transport (ITT has attracted widespread attention of scientists. This review introduces the latest theoretical development in ITT through one-dimensional (1D atomic junction model to address the thermal transport across an interface. With full consideration of the atomic structures in interfaces, people can apply the 1D atomic junction model to investigate many properties of ITT, such as interfacial (Kapitza resistance, nonlinear interface, interfacial rectification, and phonon interference, and so on. For the ballistic ITT, both the scattering boundary method (SBM and the non-equilibrium Green’s function (NEGF method can be applied, which are exact since atomic details of actual interfaces are considered. For interfacial coupling case, explicit analytical expression of transmission coefficient can be obtained and it is found that the thermal conductance maximizes at certain interfacial coupling (harmonic mean of the spring constants of the two leads and the transmission coefficient is not a monotonic decreasing function of phonon frequency. With nonlinear interaction—phonon–phonon interaction or electron–phonon interaction at interface, the NEGF method provides an efficient way to study the ITT. It is found that at weak linear interfacial coupling, the nonlinearity can improve the ITT, but it depresses the ITT in the case of strong-linear coupling. In addition, the nonlinear interfacial coupling can induce thermal rectification effect. For interfacial materials case which can be simulated by a two-junction atomic chain, phonons show interference effect, and an optimized thermal coupler can be obtained by tuning its spring constant and atomic mass.

The incentives for burning alternative fuels in lime kilns are growing. An increasing demand on thorough investigations of alternative fuel impact on lime kiln performance have been recognized, and the purpose of this project has been to develop a lime kiln CFD model with the possibility to fire fuel oil and lignin. The second part of the project consists of three technical studies. Simulated data from a one-dimensional steady state program has been used to support theories on the impact of biofuels and lime mud dryness. The CFD simulations was carried out in the commercial code FLUENT. Due to difficulties with the convergence of the model the calcination reaction is not included. The model shows essential differences between the two fuels. Lignin gives a different flame shape and a longer flame length compared to fuel oil. Mainly this depends on how the fuel is fed into the combustion chamber and how much combustion air that is added as primary and secondary air. In the case of lignin combustion the required amount of air is more than in the fuel oil case. This generates more combustion gas and a different flow pattern is created. Based on the values from turbulent reaction rate for the different fuels an estimated flame length can be obtained. For fuel oil the combustion is very intense with a sharp peak in the beginning and a rapid decrease. For lignin the combustion starts not as intense as for the fuel oil case and has a smoother shape. The flame length appears to be approximately 2-3 meter longer for lignin than for fuel oil based on turbulent reaction rate in the computational simulations. The first technical study showed that there are many benefits of increasing dry solids content in the lime mud going into a kiln such as increased energy efficiency, reduced TRS, and reduced sodium in the kiln. However, data from operating kilns indicates that these benefits can be offset by increasing exit gas temperature that can limit kiln production capacity. Simulated

Full Text Available We investigate the initial-value problem for one-dimensional compressible fluids endowed with internal capillarity. We focus on the isothermal inviscid case with variable capillarity. The resulting equations for the density and the velocity, consisting of the mass conservation law and the momentum conservation with Korteweg stress, are a system of third order nonlinear dispersive partial differential equations. Additionally, this system is Hamiltonian and admits travelling solutions, representing propagating phase boundaries with internal structure. By change of unknown, it roughly reduces to a quasilinear Schrodinger equation. This new formulation enables us to prove local well-posedness for smooth perturbations of travelling profiles and almost-global existence for small enough perturbations. A blow-up criterion is also derived.

The evolution of a particle which moves on a discrete one-dimensional lattice, according to a random walk low, approximates better the diffusion process smaller the steps of the spatial lattice and time are. For a sufficiently large assembly of particles one can assume that their relative frequency at lattice knots approximates the distribution function of the diffusion process. This assumption has been tested by simulating on computer two analytical solutions of the diffusion equation: the Brownian motion and the steady state linear distribution. To evaluate quantitatively the similarity between the numerical and analytical solutions we have used a norm given by the absolute value of the difference of the two solutions. Also, a diffusion coefficient at any lattice knots and moment of time has been calculated, by using the numerical solution both from the diffusion equation and the particle flux given by Fick's low. The difference between diffusion coefficient of analytical solution and the spatial lattice mean coefficient of numerical solution constitutes another quantitative indication of the similarity of the two solutions. The results obtained show that the approximation depends first on the number of particles at each knot of the spatial lattice. In conclusion, the random walk is a microscopic process of the molecular dynamics type which permits simulations precision of the diffusion processes with given precision. The numerical method presented in this work may be useful both in the analysis of real experiments and for theoretical studies

Full Text Available - dimensional models presented here may illuminate the study of more realistic models. For the model in which as many dislocations are poised for backward jumps as for forward jumps, the experimental activation volume Vye(C27a) under applied stresses close to C...27a is different from the true activation volume V(C27) evaluated at C27 ?C27a. The relations between the two are developed. A model is then discussed in which fewer dislocations are available for backward than for forward jumps. Finally...

Full Text Available If a turbine loses its electrical load, it will rotate freely and increase speed, eventually achieving that rotational speed which produces zero net torque. This is known as a runaway situation. Unlike many other types of turbine, a marine current turbine will typically overshoot the final runaway speed before slowing down and settling at the runaway speed. Since the hydrodynamic forces acting on the turbine are dependent on rotational speed and acceleration, turbine behaviour during runaway becomes important for load analyses during turbine design. In this article, we consider analytical and numerical models of marine current turbine runaway behaviour in one dimension. The analytical model is found not to capture the overshoot phenomenon, while still providing useful estimates of acceleration at the onset of runaway. The numerical model incorporates turbine wake build-up and predicts a rotational speed overshoot. The predictions of the models are compared against measurements of runaway of a marine current turbine. The models are also used to recreate previously-published results for a tidal turbine and applied to a wind turbine. It is found that both models provide reasonable estimates of maximum accelerations. The numerical model is found to capture the speed overshoot well.

We propose a stochastic particle model in (1+1) dimensions, with one dimension corresponding to rapidity and the other one to the transverse size of a dipole in QCD, which mimics high-energy evolution and scattering in QCD in the presence of both saturation and particle-number fluctuations, and hence of pomeron loops. The model evolves via non-linear particle splitting, with a non-local splitting rate which is constrained by boost-invariance and multiple scattering. The splitting rate saturates at high density, so like the gluon emission rate in the JIMWLK evolution. In the mean field approximation obtained by ignoring fluctuations, the model exhibits the hallmarks of the BK equation, namely a BFKL-like evolution at low density, the formation of a traveling wave, and geometric scaling. In the full evolution including fluctuations, the geometric scaling is washed out at high energy and replaced by diffusive scaling. It is likely that the model belongs to the universality class of the reaction-diffusion process. The analysis of the model sheds new light on the pomeron loops equations in QCD and their possible improvements

In this report, some classical and new simplifications in mathematical and numerical models for river morphology are compared for conditions representing rivers in mountainous areas (high values of Froude numbers and relatively large values of sediment transport rates). Options for simplification

We use the vector recursion method of Haydock to obtain the transmittance of a class of generalized Aubry models in one-dimension. We also study the phase change of the wavefunctions as they travel through the chain and also the behaviour of the conductance with changes in size. (author). 10 refs, 9 figs

This paper is an overview of the notions of mathematical and simulation model applied to flows met in pipes and in several components of power plants (valves, pump, turbine). Finally, the results of a computer code based on the equations previously presented are given [fr

The asymmetric exclusion model describes a system of particles hopping in a preferred direction with hard core repulsion. These particles can be thought of as charged particles in a field, as steps of an interface, as cars in a queue. Several exact results concerning the steady state of this system have been obtained recently. The solution consists of representing the weights of the configurations in the steady state as products of non-commuting matrices. (author)

The 1D isotropic s = ½XY-model ( N sites), with random exchange interaction in a transverse random field is considered. The random variables satisfy bimodal quenched distributions. The solution is obtained by using the Jordan-Wigner fermionization and a canonical transformation, reducing the problem to diagonalizing an N × N matrix, corresponding to a system of N noninteracting fermions. The calculations are performed numerically for N = 1000, and the field-induced magnetization at T = 0 is obtained by averaging the results for the different samples. For the dilute case, in the uniform field limit, the magnetization exhibits various discontinuities, which are the consequence of the existence of disconnected finite clusters distributed along the chain. Also in this limit, for finite exchange constants J A and J B, as the probability of J A varies from one to zero, the saturation field is seen to vary from Γ A to Γ B, where Γ A(Γ B) is the value of the saturation field for the pure case with exchange constant equal to J A(J B) .

We study two particles in an infinite chain and coupled to phonons by interactions that modulate their hopping as described by the Peierls/Su-Schrieffer-Heeger (SSH) model. In the case of hard-core bare particles, we show that exchange of phonons generates effective nearest-neighbor repulsion between particles and also gives rise to interactions that move the pair as a whole. The two-polaron phase diagram exhibits two sharp transitions, leading to light dimers at strong coupling and the flattening of the dimer dispersion at some critical values of the parameters. This dimer (quasi)self-trapping occurs at coupling strengths where single polarons are mobile. On the other hand, in the case of soft-core particles/ spinfull fermions, we show that phonon-mediated interactions are attractive and result in strongly bound and mobile bipolarons in a wide region of parameter space. This illustrates that, depending on the strength of the phonon-mediated interactions and statistics of bare particles, the coupling to phonons may completely suppress or strongly enhance quantum transport of correlated particles. This work was supported by NSERC of Canada and the Stewart Blusson Quantum Matter Institute.

By the first approximation analyzing stability conditions of unperturbed solution of one-dimensional dynamic model with magnetic interaction between two superconducting rings obtained. The stability region in the frozen magnetic flux parameters space was constructed.

Three one-dimensional (1D) numerical wave models are evaluated for wave transformation over reefs and estimates of wave setup, runup, and ponding levels in an island setting where the beach is fronted by fringing reef and lagoons...

A radially dependent, angular leakage correction was applied to a one-dimensional, multigroup neutron diffusion theory computer code to accurately model hemispherical geometry. This method allows the analyst to model hemispherical geometry, important in nuclear criticality safety analyses, with one-dimensional computer codes, which execute very quickly. Rapid turnaround times for scoping studies thus may be realized. This method uses an approach analogous to an axial leakage correction in a one-dimensional cylinder calculation. The two-dimensional Laplace operator was preserved in spherical geometry using a leakage correction proportional to 1/r 2 , which was folded into the one-dimensional spherical calculation on a mesh-by-mesh basis. Hemispherical geometry is of interest to criticality safety because of its similarity to piles of spilled fissile material and accumulations of fissile material in process containers. A hemisphere also provides a more realistic calculational model for spilled fissile material than does a sphere

In this thesis a one-dimensional (1-D) plasma and impurity transport model is developed to address issues related to impurity behavior in Reversed Field Pinch (RFP) fusion plasmas. A coronal non-equilibrium model is used for impurities. The impurity model is incorporated into an existing onedimensional plasma transport model creating a multi-species plasma transport model which treats the plasma and impurity evolution self-consistently. Neutral deuterium particles are treated using a one-dimensional (slab) model of neutral transport. The resulting mode, RFPBI, is then applied to existing RFP devices such as ZT-40M and MST, and also to examine steady state behavior of ZTH based on the design parameters. A parallel algorithm for the impurity transport equations is implemented and tested to determine speedup and efficiency

The one-dimensional Holstein model of spinless fermions interacting with dispersionless phonons is solved by using a Bethe ansatz in analogue to that for the one-dimensional spinless Fermi-Hubbard model. Excitation energies and the corresponding wavefunctions of the model are determined by a set of partial differential equations. It is shown that the model is, at least, quasi-exactly solvable for the two-site case, when the phonon frequency, the electron-phonon coupling strength and the hopping integral satisfy certain relations. As examples, some quasi-exact solutions of the model for the two-site case are derived. (letter to the editor)

We show a simple but computationally exact solution of the one-dimensional plasma model, so-called 'Dawson's model'. Using this solution, we can describe the evolution of the plasma and find the relative stabilization of a big hole after the instability of two streams. (author)

A one-dimensional (1D) laminar oscillating flow heat transfer model is derived and applied to parallel-plate thermoacoustic heat exchangers. The model can be used to estimate the heat transfer from the solid wall to the acoustic medium, which is required for the heat input/output of thermoacoustic

The authors derive a regularized exactly solvable version of the one-dimensional quantum sine-Gordon model proceeding from the exact solution of the U(1)-symmetric Thirring model. The ground state and the excitation spectrum are obtained in the region ν 2 < 8π. (Auth.)

To understand the connection between the dynamics of microscopic turbulence and the macroscale power scaling in the L-I-H transition in magnetically confined plasmas, a new time-dependent, one-dimensional (in radius) model has been developed. The model investigates the radial force balance equati...

We present a simple one-dimensional quantum-mechanical model for a particle attached to a surface. It leads to the Schrodinger equation for a harmonic oscillator bounded on one side that we solve in terms of Weber functions and discuss the behaviour of the eigenvalues and eigenfunctions. We derive the virial theorem and other exact relationships…

By solving numerically a one-dimensionalmodel, the range of validity of some non-perturbative treatments proposed for the problem of atomic ionization by strong laser fields is examined. Some scalling properties of the ionization probability are stablished and a new approximation, which converges to the exact results in the limit of very strong fields is proposed. (Author) [pt

A one-dimensional map is proposed for modeling some of the neuronal activities, including different spiking and bursting behaviors. The model is obtained by applying some modifications on the well-known Logistic map and is named the Modified and Confined Logistic (MCL) model. Map-based neuron models are known as phenomenological models and recently, they are widely applied in modeling tasks due to their computational efficacy. Most of discrete map-based models involve two variables representing the slow-fast prototype. There are also some one-dimensional maps, which can replicate some of the neuronal activities. However, the existence of four bifurcation parameters in the MCL model gives rise to reproduction of spiking behavior with control over the frequency of the spikes, and imitation of chaotic and regular bursting responses concurrently. It is also shown that the proposed model has the potential to reproduce more realistic bursting activity by adding a second variable. Moreover the MCL model is able to replicate considerable number of experimentally observed neuronal responses introduced in Izhikevich (2004) [23]. Some analytical and numerical analyses of the MCL model dynamics are presented to explain the emersion of complex dynamics from this one-dimensional map

We consider onedimensional coupled classical-quantum models for quantum semiconductor device simulations. The coupling occurs in the space variable : the domain of the device is divided into a region with strong quantum effects (quantum zone) and a region where quantum effects are negligible (classical zone). In the classical zone, transport in diffusive approximation is modeled through diffusive limits of the Boltzmann transport equation. This leads to a hierarchy of classical model. The qu...

The numerical prediction of fog requires a very high vertical resolution of the atmosphere. Owing to a prohibitive computational effort of high resolution three dimensional models, operational fog forecast is usually done by means of onedimensional fog models. An important condition for a successful fog forecast with onedimensionalmodels consists of the proper integration of observational data into the numerical simulations. The goal of the present study is to introduce new methods for the consideration of these data in the onedimensional radiation fog model PAFOG. First, it will be shown how PAFOG may be initialized with observed visibilities. Second, a nudging scheme will be presented for the inclusion of measured temperature and humidity profiles in the PAFOG simulations. The new features of PAFOG have been tested by comparing the model results with observations of the German Meteorological Service. A case study will be presented that reveals the importance of including local observations in the model calculations. Numerical results obtained with the modified PAFOG model show a distinct improvement of fog forecasts regarding the times of fog formation, dissipation as well as the vertical extent of the investigated fog events. However, model results also reveal that a further improvement of PAFOG might be possible if several empirical model parameters are optimized. This tuning can only be realized by comprehensive comparisons of model simulations with corresponding fog observations.

A computer-oriented analytical method for predicting the rewetting rate of a hot dry wall is proposed. The wall, which is modeled as a thin flat plate with internal heat generation, receives a variable heat flux from one side while it is cooled from the other side. The model accounts for the large variations of the heat transfer coefficient near the wet front and for the temperature dependence of the thermal and physical properties of the wall. The one-dimensional heat-conduction equation is solved by dividing the quenching zone into small segments of arbitrary temperature increment and constant properties and heat transfer coefficient. A trial-and-error method is developed to predict the velocity of the wet front, the length of the quenching zone and the temperature profile. The one-dimensionalmodels of other authors can be obtained as particular cases of the present model. (Auth.)

Full Text Available A one-dimensionalmodel is developed to represent the ash-melting phenomenon, which was not considered in the previous one-dimensional (1-D entrained-flow gasifier model. We include sensible heat of slag and the fusion heat of ash in the heat balance equation. To consider the melting of ash, we propose an algorithm that calculates the energy balance for three scenarios based on temperature. We also use the composition and the thermal properties of anorthite mineral to express ash. gPROMS for differential equations is used to solve this algorithm in a simulation; the results include coal conversion, gas composition, and temperature profile. Based on the Texaco pilot plant gasifier, we validate our model. Our results show good agreement with previous experimental data. We conclude that the sensible heat of slag and the fusion heat of ash must be included in the entrained flow gasifier model.

We calculate the density of states and various characteristic lengths of the one-dimensional Anderson model in the limit of weak disorder. All these quantities show anomalous fluctuations near the band centre. This has already been observed for the density of states in a different model by Gorkov and Dorokhov, and is in close agreement with a Monte-Carlo calculation for the localization length by Czycholl, Kramer and Mac-Kinnon.

Single-fluid transport in the plasma scrape-off layer is modeled for poloidal divertor and mechanically limited discharges. This numerical model is one-dimensional along a field line and time-independent. Conductive and convective transport, as well as impurity and neutral source (sink) terms are included. A simple shooting method technique is used for obtaining solutions. Results are shown for the case of the proposed Alcator DCT tokamak

An analytical model is presented for the electronic conductance in a onedimensional atomic chain across an isolated defect. The model system consists of two semi infinite lead atomic chains with the defect atom making the junction between the two leads. The calculation is based on a linear combination of atomic orbitals in the tight-binding approximation, with a single atomic one s-like orbital chosen in the present case. The matching method is used to derive analytical expressions for the scattering cross sections for the reflection and transmission processes across the defect, in the Landauer-Buttiker representation. These analytical results verify the known limits for an infinite atomic chain with no defects. The model can be applied numerically for onedimensional atomic systems supported by appropriate templates. It is also of interest since it would help establish efficient procedures for ensemble averages over a field of impurity configurations in real physical systems.

Presently there is quite few analytical studies in traffic flows due to the non-linearity of the governing equations. In the present paper we introduce travelling wave solutions for the homogeneous one-dimensional FREFLO model, which are expressed in the form of series and describe the procedure that vehicles/pedestrians move with a negative velocity and decelerate until rest, then accelerate inversely to positive velocities. This method is expect to be extended to more complex situations in the future.

The one-dimensional Hubbard model with arbitrary boundary magnetic fields is solved exactly via the Bethe ansatz methods. With the coordinate Bethe ansatz in the charge sector, the second eigenvalue problem associated with the spin sector is constructed. It is shown that the second eigenvalue problem can be transformed into that of the inhomogeneous XXX spin chain with arbitrary boundary fields which can be solved via the off-diagonal Bethe ansatz method.

In the present paper, a methodology to calculate the heat fluxes inside a turbocharger from diesel passenger car is presented. The heat transfer phenomenon is solved by using a onedimensional lumped model that takes into account both the heat fluxes between the different turbocharger elements, as well as the heat fluxes between the working fluids and the turbocharger elements. This heat transfer study is supported by the high temperature differences between the working fluids passing thr...

The ground-state configurations of the one-dimensional Falicov-Kimball model are studied exactly with numerical calculations revealing unexpected effects for small interaction strength. In neutral systems we observe molecular formation, phase separation, and changes in the conducting properties; while in nonneutral systems the phase diagram exhibits Farey tree order (Aubry sequence) and a devil's staircase structure. Conjectures are presented for the boundary of the segregated domain and the general structure of the ground states

In this paper, we apply the two-time Green's function method, and provide a simple way to study the magnetic properties of one-dimensional spin-(S,s) Heisenberg ferromagnets. The magnetic susceptibility and correlation functions are obtained by using the Tyablikov decoupling approximation. Our results show that the magnetic susceptibility and correlation length are a monotonically decreasing function of temperature regardless of the mixed spins. It is found that in the case of S=s, our results of one-dimensional mixed-spin model is reduced to be those of the isotropic ferromagnetic Heisenberg chain in the whole temperature region. Our results for the susceptibility are in agreement with those obtained by other theoretical approaches. (condensed matter: electronic structure, electrical, magnetic, and optical properties)

We present an analytical and numerical study of the competition between double and super-exchange interactions in a one-dimensionalmodel. For low super-exchange interaction energy we find phase separation between ferromagnetic and anti-ferromagnetic phases. When the super-exchange interaction energy gets larger, the conduction electrons are self-trapped within separate small magnetic polarons. These magnetic polarons contain a single electron inside two or three sites depending on the conduction electron density and form a Wigner crystallization. A new phase separation is found between these small polarons and the anti-ferromagnetic phase. Spin-glass behavior is obtained consistent with experimental results of the nickelate one-dimensional compound Y 2-x Ca x BaNiO 5 .

Transient quantum dynamics in an interacting fermion-phonon system are investigated with a focus on a charge order (CO) melting after a short optical-pulse irradiation and the roles of the quantum phonons in the transient dynamics. A spinless-fermion model in a one-dimensional chain coupled with local phonons is analyzed numerically. The infinite time-evolving block decimation algorithm is adopted as a reliable numerical method for one-dimensional quantum many-body systems. Numerical results for the photoinduced CO melting dynamics without phonons are well interpreted by the soliton picture for the CO domains. This interpretation is confirmed by numerical simulation of an artificial local excitation and the classical soliton model. In the case of large phonon frequencies corresponding to the antiadiabatic condition, CO melting is induced by propagations of the polaronic solitons with the renormalized soliton velocity. On the other hand, in the case of small phonon frequencies corresponding to the adiabatic condition, the first stage of the CO melting dynamics occurs due to the energy transfer from the fermionic to phononic systems, and the second stage is brought about by the soliton motions around the bottom of the soliton band. The analyses provide a standard reference for photoinduced CO melting dynamics in one-dimensional many-body quantum systems.

The dynamic modeling of a sliding joint on a one-dimensional medium, such as a cable or a beam, is studied in this paper. The sliding joint is implemented by positioning it at a moving node on the one-dimensional medium, which is realized by variable-length elements at either side of the joint. The variable-length element is established with an absolute nodal coordinate formulation (ANCF) in the framework of the Arbitrary Lagrange–Euler (ALE) description. The sliding of the joint is described by the increasing of the length on one side of the one-dimensional medium and a corresponding decreasing of the other side. In order to capture the discontinuity of the slopes at the position of the sliding joint, the moving node has two slopes as generalized coordinates which are equal to each other in the case of a beam but not in the case of a cable, and in order to avoid the addition–deletion constraint, the node adjacent to the moving node is added or deleted if the element is too long or too short. The governing equations for the coupled system are derived in terms of D’Alembert’s principle and the resulting equations of motion are formulated in the standard form of differential algebraic equations of multibody systems. Numerical examples are presented to validate the method proposed by comparing with analytical results which are available or are made possible by simplifying the model.

A one-dimensional (1D) laminar oscillating flow heat transfer model is derived and applied to parallel-plate thermoacoustic heat exchangers. The model can be used to estimate the heat transfer from the solid wall to the acoustic medium, which is required for the heat input/output of thermoacoustic systems. The model is implementable in existing (quasi-)1D thermoacoustic codes, such as DeltaEC. Examples of generated results show good agreement with literature results. The model allows for arbitrary wave phasing; however, it is shown that the wave phasing does not significantly influence the heat transfer.

In our work we justify the applicability of a dielectric mirror model to the description of a real photonic crystal. We demonstrate that a simple one-dimensionalmodel of a multilayer mirror can be employed for modeling of a slab waveguide with periodically changing width. It is shown that this width change can be recalculated to the effective refraction index modulation. The applicability of transfer matrix method of reflection properties calculation was demonstrated. Finally, our 1-D model was employed to analyze reflection properties of a 2-D structure - a slab photonic crystal with a number of elliptic holes.

We extend the Bethe Ansatz solution of a one-dimensional integrable fermionic model with correlated hopping to the parameter regime Δt > 1. It is found that the model is equivalent to one with interaction 2 - Δt, but with twisted boundary conditions. Apart from the ground state energy we investigate the low-lying excitations and the asymptotic behaviour of the correlation functions. As in the case of Δt < 1 we find dominating superconducting correlations for small doping. The behaviour in this regime therefore differs from that of the non-integrable model with symmetric bond-charge interaction (Hirsch model). (orig.)

We have shown that by using a model of the gas of partition points on one-dimensional lattice, we can find some results about the enzyme kinetics or the average domain-size, which we have obtained before by using a correlated Walks' theory or a probabilistic (combinatoric) way. We have discussed also the problem related with the spread of an infection of disease and the stochastic model of partition points. We think that this model, as a very simple model and mathematically transparent, can be advantageous for other theoretical investigations in chemistry or modern biology. (author). 14 refs, 6 figs, 1 tab

We study the ground-state phase diagram of the strongly interacting Harper–Hofstadter–Mott model at quarter flux on a quasi-one-dimensional lattice consisting of a single magnetic flux quantum in y-direction. In addition to superfluid phases with various density patterns, the ground-state phase diagram features quasi-one-dimensional analogs of fractional quantum Hall phases at fillings ν = 1/2 and 3/2, where the latter is only found thanks to the hopping anisotropy and the quasi-one-dimensional geometry. At integer fillings—where in the full two-dimensional system the ground-state is expected to be gapless—we observe gapped non-degenerate ground-states: at ν = 1 it shows an odd ‘fermionic’ Hall conductance, while the Hall response at ν = 2 consists of the transverse transport of a single particle–hole pair, resulting in a net zero Hall conductance. The results are obtained by exact diagonalization and in the reciprocal mean-field approximation.

A point-kinetics model with one-dimensional (radial) heat conduction formalism has been developed. The heat conduction formalism is based on corner-mesh finite difference method. To get average temperatures in various conducting regions, a novel weighting scheme has been devised. The heat conduction model has been incorporated in the point-kinetics code MRTF-FUEL. The point-kinetics equations are solved using the method of real integrating factors. It has been shown by analysing the simulation of hypothetical loss of regulation accident in NAPP reactor that the model is superior to the conventional one in accuracy and speed of computation. (author). 3 refs., 3 tabs

We present the method for determining the velocity model of the Earth's crust and the parameters of earthquakes in the Middle Kura Depression from the data of network telemetry in Azerbaijan. Application of this method allowed us to recalculate the main parameters of the hypocenters of the earthquake, to compute the corrections to the arrival times of P and S waves at the observation station, and to significantly improve the accuracy in determining the coordinates of the earthquakes. The model was constructed using the VELEST program, which calculates one-dimensional minimal velocity models from the travel times of seismic waves.

We have shown that by using a model of the partition points gas on a one-dimensional lattice, we can study, besides the saturation curves obtained before for the enzyme kinetics, also the denaturation process, i.e. the breaking of the hydrogen bonds connecting the two strands, under treatment by heat of DNA. We think that this model, as a very simple model and mathematically transparent, can be advantageous for pedagogic goals or other theoretical investigations in chemistry or modern biology. (author). 29 refs, 4 figs

The Glauber dynamics of interacting Ising spins (the single spin-flip dynamics) is generalized to p spin-flip dynamics with a simultaneous flip of up to p spins in a single configuration move. The p spin-flip dynamics is studied of the one-dimensional Ising model with uniform nearest-neighbour interaction. For this case, an exact relation is given for the time dependence of magnetization. It was found that the critical slowing down in this model could be avoided when p spin-flip dynamics with p>2 was considered. (author). 17 refs

We investigate the one-dimensional p-wave superconducting model with periodically modulated hopping and show that under time-reversal symmetry, the number of the Majorana zero modes (MZMs) strongly depends on the modulation period. If the modulation period is odd, there can be at most one MZM. However if the period is even, the number of the MZMs can be zero, one and two. In addition, the MZMs will disappear as the chemical potential varies. We derive the condition for the existence of the MZMs and show that the topological properties in this model are dramatically different from the one with periodically modulated potential.

Full Text Available A stretching behavior of knitted and woven textiles is modeled. In our work, the yarns are modeled as one-dimensional hyperelastic strings with frictional contact. Capstan law known for Coulomb’s friction of yarns is extended to an additional adhesion due to gluing of filaments on the yarn surface or some chemical reaction. Two-step Newton’s method is applied for the solution of the large stretching with sliding evolution in the contact nodes. The approach is illustrated on a hysteresis of knitted textile and on the force-strain curve for a woven pattern and both compared with experimental effective curves.

In this paper, the relation between synchronization and pattern formation in one-dimensional discrete and continuous open flow models is investigated in detail. Firstly a sufficient condition for globally asymptotical stability of lag/anticipating synchronization among lattices of these models is proved by analytic method. Then, by analyzing and simulating lag/anticipating synchronization in discrete case, three kinds of pattern of wave (it is called wave pattern) travelling in the lattices are discovered. Finally, a proper definition for these kinds of pattern is proposed

The novel method proposed by one of the authors to calculate exactly the response functions of the one-dimensional Tomonaga-model is described in more detail. The method is generalized for the case of a system of coupled chains where both the interchain and interchain interactions have forward scattering components only. The model does not show real phase transition at any finite temperature indicating that the interchain backward scattering or hopping is needed to have an ordering of the chains at finite temperature. (author)

In this Letter, we present a one-dimensionalmodel that includes a hard core at the origin, a Dirac delta barrier at a point in the positive semiaxis and a mass jump at the same point. We study the effect of this mass jump in the behavior of the resonances of the model. We obtain an infinite number of resonances for this situation, showing that for the case of a mass jump the imaginary part of the resonance poles tend to a fixed value depending on the quotient of masses, and demonstrate that none of these resonances is degenerated.

A derivation of turbulent flow parameters, combined with data from erosive burning test motors and blowing wall tests results in erosive burning model candidates useful in one-dimensional internal ballistics analysis capable of scaling across wide ranges of motor size. The real-time burn rate data comes from three test campaigns of subscale segmented solid rocket motors tested at two facilities. The flow theory admits the important effect of the blowing wall on the turbulent friction coefficient by using blowing wall data to determine the blowing wall friction coefficient. The erosive burning behavior of full-scale motors is now predicted more closely than with other recent models.

Type I X-ray bursts are caused by thermonuclear explosions occurring on the surface of an accreting neutron star in a binary star system. Observations and simulations of these phenomena are of great importance for understanding the fundamental properties of neutron stars and dense matter because the equation of state for cold dense matter can be constrained by the mass-radius relationship of neutron stars. During the bursts, turbulence plays a key role in mixing the fuels and driving the unstable nuclear burning process. This dissertation presents one-dimensionalmodels of photospheric radius expansion bursts with a new approach to simulate turbulent advection. Compared with the traditional mixing length theory, the one-dimensional turbulence (ODT) model represents turbulent motions by a sequence of maps that are generated according to a stochastic process. The light curves I obtained with the ODT models are in good agreement with those of the KEPLER model in which the mixing length theory and various diffusive processes are applied. The abundance comparison, however, indicates that the differences in turbulent regions and turbulent diffusivities result in more 12C survival during the bursts in the ODT models, which can make a difference in the superbursts phenomena triggered by unstable carbon burning.

The Q-machine is a nontrivial bounded plasma system which is excellently suited not only for fundamental plasma physics investigations but also for the development and testing of new theoretical methods for modeling such systems. However, although Q-machines have now been around for over thirty years, it appears that there exist no comprehensive theoretical models taking into account their considerable geometrical and physical complexity with a reasonable degree of self-consistency. In the present context we are concerned with the low-density, single-emitter Q-machine, for which the most widely used model is probably the (one-dimensional) ''collisionless plane-diode model'', which has originally been developed for thermionic diodes. Although the validity of this model is restricted to certain ''axial'' phenomena, we consider it a suitable starting point for extensions of various kinds. While a generalization to two-dimensional geometry (with still collisionless plasma) is being reported elsewhere, the present work represents a first extension to collisional plasma (with still one-dimensional geometry). (author) 12 refs., 2 figs

Type I X-ray bursts are caused by thermonuclear explosions occurring on the surface of an accreting neutron star in a binary star system. Observations and simulations of these phenomena are of great importance for understanding the fundamental properties of neutron stars and dense matter because the equation of state for cold dense matter can be constrained by the mass-radius relationship of neutron stars. During the bursts, turbulence plays a key role in mixing the fuels and driving the unstable nuclear burning process. This dissertation presents one-dimensionalmodels of photospheric radius expansion bursts with a new approach to simulate turbulent advection. Compared with the traditional mixing length theory, the one-dimensional turbulence (ODT) model represents turbulent motions by a sequence of maps that are generated according to a stochastic process. The light curves I obtained with the ODT models are in good agreement with those of the KEPLER model in which the mixing length theory and various diffusive processes are applied. The abundance comparison, however, indicates that the differences in turbulent regions and turbulent diffusivities result in more 12 C survival during the bursts in the ODT models, which can make a difference in the superbursts phenomena triggered by unstable carbon burning.

Since the LLNL one-dimensional coupled transport and chemical kinetics model of the troposphere and stratosphere was originally developed in 1972 (Chang et al., 1974), there have been many changes to the model's representation of atmospheric physical and chemical processes. A brief description is given of the current LLNL one-dimensional coupled transport and chemical kinetics model of the troposphere and stratosphere

This paper presents the simulation of a planar, one-dimensional expanding Ge x-ray laser plasma using a new code which combines hydrodynamics, laser absorption, and detailed level population calculations within the same simulation. Previously, these simulations were performed in separate steps. We will present the effect of line transfer on gains and excited level populations and compare the line transfer result with simulations using escape probabilities. We will also discuss the impact of different atomic models on the accuracy of our simulation

In order to simplify condensation heat transfer calculations for feed water heaters, onedimensional (1D) analyses were compared with three dimensional (3D) analyses. The results showed that average condensation heat transfer coefficients by 1D analyses with 1/2 rows of heat transfer tubes agreed with those by 3D analyses within 7%. Using the 1D analysis model, effects of the pitch of heat transfer tubes were evaluated. The results showed that the pitch did not affect much on heat transfer rates and that the size of heat transfer tube bundle could be decreased by a small pitch. (author)

The results reported here, using a standard numerical algorithm and a simple low temperature extrapolation, appear consistent with numerical results of Sorella et al. for the one-dimensional Hubbard model in the half-filled and quarter-filled band cases. However, it is argued that the discontinuity at the Fermi level found in the quarter-filled case is likely to come from the zero-temperature finite-size dependence of the quasiparticle weight Z, which is also discussed here. (orig.).

We analyze the so-called Kovacs effect in the one-dimensional Ising model with Glauber dynamics. We consider small enough temperature jumps, for which a linear response theory has been recently derived. Within this theory, the Kovacs hump is directly related to the monotonic relaxation function of the energy. The analytical results are compared with extensive Monte Carlo simulations, and an excellent agreement is found. Remarkably, the position of the maximum in the Kovacs hump depends on the fact that the true asymptotic behavior of the relaxation function is different from the stretched exponential describing the relevant part of the relaxation at low temperatures.

Mathematical modeling for composition distal arterial pulse wave in the blood vessels of the upper limbs was considered. Formation of distal arterial pulse wave is represented as a composition of forward and reflected pulse waves propagating along the arterial vessels. The formal analogy between pulse waves propagation along the human arterial system and the propagation of electrical oscillations in electrical transmission lines with distributed parameters was proposed. Dependencies of pulse wave propagation along the human arterial system were obtained by solving the one-dimensional Navier-Stokes equations for a few special cases.

The effect of the chain and the dimer anisotropies on the ground state energy and the energy gap of the spin-1/2 quasi-one-dimensional antiferromagnetic Heisenberg model is investigated using a mean field theory. The dependence of the magnetization and the effective hopping parameters on the anisotropy α xy (=J xy perpendicular /J xy parallel ) are presented for several values of the chain anisotropy. However, such a system exhibits a transition from antiferromagnetic ordered to disordered phases for arbitrary chain anisotropy and dimer anisotropy. (author). 22 refs, 11 figs

The phase transition of the classical two-dimensional Potts model, in particular the order of the transition as the number of components q increases, is studied by constructing renormalization group transformations on the equivalent one-dimensional quatum problem. It is shown that the block transformation with two sites per cell indicates the existence of a critical qsub(c) separating the small q and large q regions with different critical behaviours. The physically accessible fixed point for q>qsub(c) is a discontinuity fixed point where the specific heat exponent α=1 and therefore the transition is of first order. (author)

Plasma-immersion ion implantation (also known as plasma-source ion implantation) is a process in which a target is immersed in a plasma and a series of large negative-voltage pulses are applied to it to extract ions from the plasma and implant them into the target. A general one-dimensionalmodel is developed to study this process in different coordinate systems for the case in which the pressure of the neutral gas is large enough that the ion motion in the sheath can be assumed to be highly collisional

We propose a thermodynamic version of the Axelrod model of social influence. In one-dimensional (1D) lattices, the thermodynamic model becomes a coupled Potts model with a bonding interaction that increases with the site matching traits. We analytically calculate thermodynamic and critical properties for a 1D system and show that an order-disorder phase transition only occurs at T=0 independent of the number of cultural traits q and features F. The 1D thermodynamic Axelrod model belongs to the same universality class of the Ising and Potts models, notwithstanding the increase of the internal dimension of the local degree of freedom and the state-dependent bonding interaction. We suggest a unifying proposal to compare exponents across different discrete 1D models. The comparison with our Hamiltonian description reveals that in the thermodynamic limit the original out-of-equilibrium 1D Axelrod model with noise behaves like an ordinary thermodynamic 1D interacting particle system.

We introduce and analyze an exactly soluble one-dimensional Ising model with long range interactions that exhibits a mixed-order transition, namely a phase transition in which the order parameter is discontinuous as in first order transitions while the correlation length diverges as in second order transitions. Such transitions are known to appear in a diverse classes of models that are seemingly unrelated. The model we present serves as a link between two classes of models that exhibit a mixed-order transition in one dimension, namely, spin models with a coupling constant that decays as the inverse distance squared and models of depinning transitions, thus making a step towards a unifying framework.

Over the last few decades Greenland Ice sheet mass balance has become increasingly negative, caused by enhanced surface melting and speedup of the marine-terminating outlet glaciers at the ice sheet margins. Glaciers speedup has been related, among other factors, to enhanced submarine melting, which in turn is caused by warming of the surrounding ocean and less obviously, by increased subglacial discharge. While ice-ocean processes potentially play an important role in recent and future mass balance changes of the Greenland Ice Sheet, their physical understanding remains poorly understood. In this work we performed numerical experiments with a one-dimensional plume model coupled to a one-dimensional iceflow model. First we investigated the sensitivity of submarine melt rate to changes in ocean properties (ocean temperature and salinity), to the amount of subglacial discharge and to the glacier's tongue geometry itself. A second set of experiments investigates the response of the coupled model, i.e. the dynamical response of the outlet glacier to altered submarine melt, which results in new glacier geometry and updated melt rates.

A three species one-dimensional kinetic model is presented for a spatially homogeneous weakly ionized plasma subjected to the action of a time varying electric field. Planar geometry is assumed, which means that the plasma evolves in the privileged direction of the field. The energy transmitted to the electric charges is channelized to the neutrals thanks to collisions, a mechanism that influences the plasma dynamics. Charge-charge interactions have been designed as a one-dimensional collision term equivalent to the Landau operator used for fully ionized plasmas. Charge-neutral collisions are modelled by a conservative drift-diffusion operator in the Dougherty's form. The resulting set of coupled integro-differential equations is solved with the stable and robust propagator integral method. This semi–analytical method feasibility accounts for non–linear effects without appealing to linearisation or simplifications, providing conservative physically meaningful solutions even for initial or emerging sharp velocity distribution function profiles. It is found that charge-neutral collisions exert a significant effect since a quite different plasma evolution arises if compared to the collisionless limit. In addition, substantial differences in the system motion are found for constant and temperature dependent collision frequencies cases.

We show that single-file water in nanopores can be viewed as a one-dimensional (1D) Ising model, and we investigate, on the basis of this, the static dielectric response of a chain of hydrogen-bonded water molecules to an external field. To achieve this, we use a recently developed dipole lattice model that accurately captures the free energetics of nanopore water. In this model, the total energy of the system can be expressed as the sum of the effective interactions of chain ends and orientational defects. Neglecting these interactions, we essentially obtain the 1D Ising model, which allows us to derive analytical expressions for the free energy as a function of the total dipole moment and for the dielectric susceptibility. Our expressions, which agree very well with simulation results, provide the basis for the interpretation of future dielectric spectroscopy experiments on water-filled nanopore membranes.

A user oriented computer model (TRANS1D) was developed for application to the analysis and design of vertical soil-bentonite barriers. TRANS1D is a collection of analytical and numerical solutions to the onedimensional advective-dispersive-reactive (ADR) equation. The primary objective in developing TRANS1D was to enable the designer of a barrier system to evaluate the potential system performance with respect to contaminant transport, without performing difficult and time consuming field or laboratory experiments. Several issues related to model application are discussed, including identification of governing transport processes, specification of boundary conditions, and parameter estimation. Model predictions are compared with the results of laboratory column experiments conducted with soil bentonite barrier material under diffusion-dominated conditions. Good agreement between model calibrations and experimental results was noted, with calibrated diffusion coefficients for organic contaminants consistent with literature values

Recently Seevinck and Uffink argued that genuine multipartite entanglement (GME) had not been established in the experiments designed to confirm GME. In this paper, we use the Bell-type inequalities introduced by Seevinck and Svetlichny [M. Seevinck, G. Svetlichny, Phys. Rev. Lett. 89 (2002) 060401] to investigate the GME problem in the one-dimensional transverse-field Ising model. We show explicitly that the ground states of this model violate the inequality when the external transverse magnetic field is weak, which indicate that the ground states in this model with weak magnetic field are fully entangled. Since this model can be simulated with nuclear magnetic resonance, our results provide a fresh approach to experimental test of GME.

The Axelrod model for the dissemination of culture exhibits a rich spatial distribution of cultural domains, which depends on the values of the two model parameters: F, the number of cultural features and q, the common number of states each feature can assume. In the one-dimensionalmodel with F = q = 2, which is closely related to the constrained voter model, Monte Carlo simulations indicate the existence of multicultural absorbing configurations in which at least one macroscopic domain coexist with a multitude of microscopic ones in the thermodynamic limit. However, rigorous analytical results for the infinite system starting from the configuration where all cultures are equally likely show convergence to only monocultural or consensus configurations. Here we show that this disagreement is due simply to the order that the time-asymptotic limit and the thermodynamic limit are taken in the simulations. In addition, we show how the consensus-only result can be derived using Monte Carlo simulations of finite chains.

“Pure” homogeneous linear supermultiplets (minimal and non-minimal) of the N=4-extended one-dimensional supersymmetry algebra are classified. “Pure” means that they admit at least one graphical presentation (the corresponding graph/graphs are known as “Adinkras”). We further prove the existence of “entangled” linear supermultiplets which do not admit a graphical presentation, by constructing an explicit example of an entangled N=4 supermultiplet with field content (3, 8, 5). It interpolates between two inequivalent pure N=4 supermultiplets with the same field content. The one-dimensional N=4 sigma-model with a three-dimensional target based on the entangled supermultiplet is presented. The distinction between the notion of equivalence for pure supermultiplets and the notion of equivalence for their associated graphs (Adinkras) is discussed. Discrete properties such as “chirality” and “coloring” can discriminate different supermultiplets. The tools used in our classification include, among others, the notion of field content, connectivity symbol, commuting group, node choice group, and so on.

We investigate the effect of dimensional crossover in the ground state of the antiferromagnetic spin-1 Heisenberg model on the anisotropic triangular lattice that interpolates between the regime of weakly coupled Haldane chains (J'≪J ) and the isotropic triangular lattice (J'=J ). We use the density-matrix renormalization group (DMRG) and Schwinger boson theory performed at the Gaussian correction level above the saddle-point solution. Our DMRG results show an abrupt transition between decoupled spin chains and the spirally ordered regime at (J'/J) c˜0.42 , signaled by the sudden closing of the spin gap. Coming from the magnetically ordered side, the computation of the spin stiffness within Schwinger boson theory predicts the instability of the spiral magnetic order toward a magnetically disordered phase with one-dimensional features at (J'/J) c˜0.43 . The agreement of these complementary methods, along with the strong difference found between the intra- and the interchain DMRG short spin-spin correlations for sufficiently large values of the interchain coupling, suggests that the interplay between the quantum fluctuations and the dimensional crossover effects gives rise to the one-dimensionalization phenomenon in this frustrated spin-1 Hamiltonian.

In the present work, we make use of simplified nonlinear models based on the compressibility factor (Peter et al., Phys. Plasmas, vol. 20 (12), 2013, 123104) to predict the gain of one-dimensional (1-D) free-electron lasers (FELs), considering space-charge and thermal effects. These models proved to be reasonable to estimate some aspects of 1-D FEL theory, such as the position of the onset of mixing, in the case of a initially cold electron beam, and the position of the breakdown of the laminar regime, in the case of an initially warm beam (Peter et al., Phys. Plasmas, vol. 21 (11), 2014, 113104). The results given by the models are compared to wave-particle simulations showing a reasonable agreement.

A simple two-particle model of excitons in conjugated polymers is proposed as an alternative to usual highly computationally demanding quantum chemical methods. In the two-particle model, the exciton is described as an electron-hole pair interacting via Coulomb forces and confined to the polymer backbone by rigid walls. Furthermore, by integrating out the transverse part, the two-particle equation is reduced to one-dimensional form. It is demonstrated how essentially exact solutions are obtained in the cases of short and long conjugation length, respectively. From a linear combination of these cases an approximate solution for the general case is obtained. As an application of the model the influence of a static electric field on the electron-hole overlap integral and exciton energy is considered.

HECTR (Hydrogen Event Containment Transient Response) is a lumped-parameter computer code developed for calculating the pressure-temperature response to combustion in a nuclear power plant containment building. The code uses a control-volume approach and subscale models to simulate the mass, momentum, and energy transfer occurring in the containment during a loss-of-collant-accident (LOCA). This document describes one-dimensional subscale models for mass and momentum transfer, and the modifications to the code required to implement them. Two problems were analyzed: the first corresponding to a standard problem studied with previous HECTR versions, the second to experiments. The performance of the revised code relative to previous HECTR version is discussed as is the ability of the code to model the experiments. 8 refs., 5 figs., 3 tabs

We propose a mean to obtain computationally useful resource states also known as cluster states, for measurement-based quantum computation, via transitionless quantum driving algorithm. The idea is to cool the system to its unique ground state and tune some control parameters to arrive at computationally useful resource state, which is in one of the degenerate ground states. Even though there is set of conserved quantities already present in the model Hamiltonian, which prevents the instantaneous state to go to any other eigenstate subspaces, one cannot quench the control parameters to get the desired state. In that case, the state will not evolve. With involvement of the shortcut Hamiltonian, we obtain cluster states in fast-forward manner. We elaborate our proposal in the one-dimensional Kitaev honeycomb model, and show that the auxiliary Hamiltonian needed for the counterdiabatic driving is of M-body interaction.

There are some particular one-dimensionalmodels, such as the Ising-Heisenberg spin models with a variety of chain structures, which exhibit unexpected behaviors quite similar to the first and second order phase transition, which could be confused naively with an authentic phase transition. Through the analysis of the first derivative of free energy, such as entropy, magnetization, and internal energy, a "sudden" jump that closely resembles a first-order phase transition at finite temperature occurs. However, by analyzing the second derivative of free energy, such as specific heat and magnetic susceptibility at finite temperature, it behaves quite similarly to a second-order phase transition exhibiting an astonishingly sharp and fine peak. The correlation length also confirms the evidence of this pseudo-transition temperature, where a sharp peak occurs at the pseudo-critical temperature. We also present the necessary conditions for the emergence of these quasi-phases and pseudo-transitions.

Over the past several years, a picture has emerged of transport in field reversed configuration (FRC) which explains many, though not all, of the loss phenomena observed in that device. That picture is complicated by the geometry, which includes both magnetically connected and magnetically isolated regions, and by the transport process, which includes a substantial contribution from short wavelength, fast time scale processes. This paper extends our previous work on this topic by carrying a one-dimensionalmodel as far as it can be carried, in terms of goemetrical and physical consistency, and isolates the difference between the model and experiment as coming from phenomena beyond the scope of 1-D anomalous transport

Phonons in a one-dimensional chain of charged microspheres suspended in a plasma were studied in an experiment. The phonons correspond to random particle motion in the chain; no external manipulation was applied to excite the phonons. Two modes were observed, longitudinal and transverse. The velocity fluctuations in the experiment are analyzed using current autocorrelation functions and a phonon spectrum. The phonon energy was found to be unequally partitioned among phonon modes in the dusty plasma experiment. The experimental phonon spectrum was characterized by a dispersion relation that was found to differ from the dispersion relation for externally excited phonons. This difference is attributed to the presence of frictional damping due to gas, which affects the propagation of externally excited phonons differently from phonons that correspond to random particle motion. A model is developed and fit to the experiment to explain the features of the autocorrelation function, phonon spectrum, and the dispersion relation

A one-dimensionalmodel of steady, compressible channel flow with mass, momentum and energy addition is discussed. An exact solution to the governing equations was found and from it a similarity parameter relating dimensionless mass, momentum and energy addition identified. This similarity parameter is used to make two flows having different dimensionless mass, momentum and energy additions equivalent. Application of the similarity parameter to the LASL Intense Neutron Source experiment and the Sandia simulation of that experiment results in an expression relating the dimensionless mass addition of combustible gas required in the Sandia experiment to dimensionless energy addition in the LASL experiment. Results of the analysis indicate that the Sandia experiment can realistically simulate the energy addition in the LASL Intense Neutron Source experiment

The Yang-Lee, Fisher and Potts zeros of the one-dimensional Q-state Potts model are studied using the theory of dynamical systems. An exact recurrence relation for the partition function is derived. It is shown that zeros of the partition function may be associated with neutral fixed points of the recurrence relation. Further, a general equation for zeros of the partition function is found and a classification of the Yang-Lee, Fisher and Potts zeros is given. It is shown that the Fisher zeros in a nonzero magnetic field are located on several lines in the complex temperature plane and that the number of these lines depends on the value of the magnetic field. Analytical expressions for the densities of the Yang-Lee, Fisher and Potts zeros are derived. It is shown that densities of all types of zeros of the partition function are singular at the edge singularity points with the same critical exponent

Finite-size scaling functions are investigated both for the mean-square magnetization fluctuations and for the probability distribution of the magnetization in the one-dimensional Ising model. The scaling functions are evaluated in the limit of the temperature going to zero (T {yields} 0), the size of the system going to infinity (N {yields} {infinity}) while N[1 - tanh(J/k{sub B}T)] is kept finite (J being the nearest neighbour coupling). Exact calculations using various boundary conditions (periodic, antiperiodic, free, block) demonstrate explicitly how the scaling functions depend on the boundary conditions. We also show that the block (small part of a large system) magnetization distribution results are identical to those obtained for free boundary conditions.

Through the analysis of many experimental post-dryout data, it is shown that the most probable flow regime near dryout or quench front is not annular flow but churn-turbulent flow when the mass flux is low. A correlation describing the initial droplet size just after the CHF position at low mass flux is low. A correlation describing the initial droplet size just after the CHF position at low mass flux is suggested through regression analysis. In the post-dryout region at low pressure and low flow, it is found that the suggested one-dimensional mechanistic model is not applicable when the vapor superficial velocity is very low, i. e., when the flow is bubbly or slug flow regime. This is explained by the change of main entrainment mechanism with the change of flow regime. Therefore, the suggested correlation is valid only in the churn-turbulent flow regime (j * g = 0.5 ∼ 4.5)

The critical exponent of the momentum distribution near k F , 3k F and 5k F are studied numerically for one-dimensional U → ∞ Hubbard model, using finite size systems and extrapolating them to the thermodynamic limit. Results at k F agree with earlier calculations, while at 3k F exponents less than 1 are obtained for finite size systems with extrapolation to 1 (regular behaviour) in the thermodynamic limit, in contrast to earlier analytic prediction 9/8. The distribution is regular at 5k F even for finite systems. The singularity near 3k F is interpreted as due to low energy excitations near 3k F in finite systems. (author). 18 refs, 4 figs, 1 tab

Post-test calculations for twelve selected SSRTF, SCTF and CCTF tests were performed to assess the predictive capability of the one-dimensional reflood model in the REFLA/TRAC code for core thermal behavior during the reflood in a PWR LOCA. Both core void fraction profile and clad temperature transients were predicted excellently by the REFLA/TRAC code including parameter effect of core inlet subcooling, core flooding rate, core configuration, core power, system pressure, initial clad temperature and so on. The peak clad temperature was predicted within an error of 50 K. Based on these assessment results, it is verified that the core thermal hydraulic behaviors during the reflood can be predicted excellently with the REFLA/TRAC code under various conditions where the reflood may occur in a PWR LOCA. (author)

We are concerned with the global existence and large time behavior of entropy solutions to the one-dimensional unipolar hydrodynamic model for semiconductors in the form of Euler-Poisson equations in a bounded interval. In this paper, we first prove the global existence of entropy solution by vanishing viscosity and compensated compactness framework. In particular, the solutions are uniformly bounded with respect to space and time variables by introducing modified Riemann invariants and the theory of invariant region. Based on the uniform estimates of density, we further show that the entropy solution converges to the corresponding unique stationary solution exponentially in time. No any smallness condition is assumed on the initial data and doping profile. Moreover, the novelty in this paper is about the unform bound with respect to time for the weak solutions of the isentropic Euler-Poisson system.

One-dimensional secondary settling tank models are widely used in current engineering practice for design and optimization, and usually can be expressed as a nonlinear hyperbolic or nonlinear strongly degenerate parabolic partial differential equation (PDE). Reliable numerical methods are needed to produce approximate solutions that converge to the exact analytical solutions. In this study, we introduced a reliable numerical technique, the Yee-Roe-Davis (YRD) method as the governing PDE solver, and compared its reliability with the prevalent Stenstrom-Vitasovic-Takács (SVT) method by assessing their simulation results at various operating conditions. The YRD method also produced a similar solution to the previously developed Method G and Enquist-Osher method. The YRD and SVT methods were also used for a time-to-failure evaluation, and the results show that the choice of numerical method can greatly impact the solution. Reliable numerical methods, such as the YRD method, are strongly recommended.

We apply a variational method devised for the nuclear many-body problem to the one-dimensional Hubbard model with nearest neighbor hopping and periodic boundary conditions. The test wave function consist for each state out of a single Hartree-Fock determinant mixing all the sites (or momenta) as well as the spin projections of the electrons. Total spin and linear momentum are restored by projection methods before the variation. It is demonstrated that this approach reproduces the results of exact diagonalizations for half-filled N=12 and N=14 lattices not only for the energies and occupation numbers of the ground but also of the lowest excited states rather well. Furthermore, a system of ten electrons in an N=12 lattice is investigated and, finally, an N=30 lattice is studied. In addition to energies and occupation numbers we present the spectral functions computed with the help of the symmetry-projected wave functions as well

The specific heat of an infinite one-dimensional polymer chain bearing periodically arranged side radicals connected to the even sites is studied by means of quantum transfer-matrix method based on a Ising-Heisenberg model. In the absence of the exchange interactions between side radicals and the main chain, the curves of specific heat show a round peak due to the antiferromagnetic excitations for the all antiferromagnetic interactions along the polymer chain. Considering the exchange interactions between the side radicals and the main chain, the curves of the specific heat show double-peak structure for ferromagnetic interactions between the radicals and main chain, indicating that a competition between ferromagnetic and antiferromagnetic interactions and the possibility of the occurrence of the stable ferrimagnetic state along the polymer chain

In this study, the bidirectional reflectance distribution function (BRDF) of a one-dimensional conducting rough surface and a dielectric rough surface are calculated with different frequencies and roughness values in the microwave band by using the method of moments, and the relationship between the bistatic scattering coefficient and the BRDF of a rough surface is expressed. From the theory of the parameters of the rough surface BRDF, the parameters of the BRDF are obtained using a genetic algorithm. The BRDF of a rough surface is calculated using the obtained parameter values. Further, the fitting values and theoretical calculations of the BRDF are compared, and the optimization results are in agreement with the theoretical calculation results. Finally, a reference for BRDF modeling of a Gaussian rough surface in the microwave band is provided by the proposed method. (electromagnetism, optics, acoustics, heat transfer, classical mechanics, and fluid dynamics)

Seed dispersal is a primary factor influencing plant community development, and thus plays a critical role in maintaining wetland ecosystem functioning. However, compared with fluvial seed dispersal of riparian plants, dispersal of saltmarsh plant seeds in tidal channels is much less studied due to its complex behavior, and relevant mathematical modelling is particularly lacking. In this study, we developed a one-dimensional advection-dispersion model to explore the patterns of tidal seed dispersal. Oscillatory tidal current and water depth were assumed to represent the tidal effects. An exponential decay coefficient λ was introduced to account for seed deposition and retention. Analytical solution in integral form was derived using Green's function and further evaluated using numerical integration. The developed model was applied to simulate Spartina densiflora seed dispersal in a tidal channel located at the Mad River Slough in North Humboldt Bay, California, USA, to demonstrate its practical applicability. Model predictions agree satisfactorily with field observation and simulation results from Delft3D numerical model. Sensitivity analyses were also conducted to evaluate the effects of varying calibrated parameters on model predictions. The range of the seed dispersion as well as the distribution of the seed concentration were further analyzed through statistical parameters such as centroid displacement and variance of the seed cloud together with seed concentration contours. Implications of the modelling results on tidal marsh restoration and protection, e.g., revegetation through seed addition, were also discussed through scenario analysis. The developed analytical model provides a useful tool for ecological management of tidal marshes.

Program WH was developed to calculate transient pressure and velocities in hydraulic networks. It is based on one-dimensional approximation of conservation laws of mass and momentum. the energy equation is ignored which means that heat transfer effects are no included. When calculating the velocity of pressure wave, compressibility of liquid, elasticity of pipe and possible minimal presence of gas in bubble or dissolved form are included. (author)

A generalization of the one-dimensional frustrated quantum XY model is considered in which the inter and intra-chain coupling constants of the two infinite XY (planar rotor) chains have different strengths. The model can describe the superconductor-insulator transition due to charging effects in a ladder of Josephson junctions in a magnetic field with half a flux quantum per plaquette. From a fluctuation-effective action, this transition is expected to be in the universality class of the two-dimensional classical XY-Ising model. The critical behaviour is studied using a Monte Carlo transfer matrix applied to the path-integral representation of the model and a finite-size-scaling analysis of data on small system sizes. It is found that, unlike the previous studied case of equal inter and intra-chain coupling constants, the XY and Ising-like excitations of the quantum model decouple for large interchain coupling, giving rise to pure Ising model critical behaviour for the chirality order parameter in good agreement with the results for the XY-Ising model. (author). 18 refs, 4 figs

The mechanical behaviour of solid biological tissues has long been described using models based on classical continuum mechanics. However, the classical continuum theories of elasticity and viscoelasticity cannot easily capture the continual remodelling and associated structural changes in biological tissues. Furthermore, models drawn from plasticity theory are difficult to apply and interpret in this context, where there is no equivalent of a yield stress or flow rule. In this work, we describe a novel one-dimensional mathematical model of tissue remodelling based on the multiplicative decomposition of the deformation gradient. We express the mechanical effects of remodelling as an evolution equation for the effective strain, a measure of the difference between the current state and a hypothetical mechanically relaxed state of the tissue. This morphoelastic model combines the simplicity and interpretability of classical viscoelastic models with the versatility of plasticity theory. A novel feature of our model is that while most models describe growth as a continuous quantity, here we begin with discrete cells and develop a continuum representation of lattice remodelling based on an appropriate limit of the behaviour of discrete cells. To demonstrate the utility of our approach, we use this framework to capture qualitative aspects of the continual remodelling observed in fibroblast-populated collagen lattices, in particular its contraction and its subsequent sudden re-expansion when remodelling is interrupted.

This paper presents a revised one-dimensional (1D) pulse flow modeling of twin-scroll turbocharger turbine under pulse flow operating conditions. The proposed methodology in this paper provides further consideration for the turbine partial admission performance during model characterization. This gives rise to significant improvement on the model pulse flow prediction quality compared to the previous model. The results show that a twin-scroll turbine is not operating at full admission throughout the in-phase pulse flow conditions. Instead, they are operating at unequal admission state due to disparity in the magnitude of turbine inlet flow. On the other hand, during out-of-phase pulse flow, a twin-scroll turbine is working at partial admission state for majority of the pulse cycle. An amended mathematical correlation in calculating the twin-scroll turbine partial admission characteristics is also presented in the paper. The impact of its accuracy on the pulse flow model prediction is explored. - Highlights: • Paper presents a 1D modeling for twin-scroll turbine under pulsating flow. • Predicted pulse pressure propagation is in good agreement with experimental data. • A methodology is proposed to consider the turbine partial admission performance. • Prediction shows twin-scroll turbine operates at unequal admission during in-phase flow. • During out-of-phase flow a twin-scroll turbine mainly operates at partial admission.

Procedures have been established for modeling the thermal response of TRU container walls (TRUPACT) exposed to a fire environment. The effort included simulation testing and thermal modeling of the wall material. In this study, both testing and modeling were directed at determining a one-dimensional thermal model for undamaged polyurethane foam. The foam was assumed to exist in a nonoxidizing environment and was exposed to an almost step change in surface temperature. Results indicate that if the TRU waste container wall includes a polyurethane foam (64 kg/m 3 density) of thickness greater than 20 cm and the wall is otherwise undamaged, there will be no change in the waste content temperature where the container is subjected to a surface temperature as high as 1333 K for times less than 3600 s. Further improvements are needed in the thermal model to include transpiration, better estimates of the temperature-dependent thermal conductivity, effects of damaged wall structure and radiation absorption effects for the charged foam. 10 figures

Sea surface current has a significant influence on electromagnetic (EM) backscattering signals and may constitute a dominant synthetic aperture radar (SAR) imaging mechanism. An effective EM backscattering model for a one-dimensional drifting fractal sea surface is presented in this paper. This model is used to simulate EM backscattering signals from the drifting sea surface. Numerical results show that ocean currents have a significant influence on EM backscattering signals from the sea surface. The normalized radar cross section (NRCS) discrepancies between the model for a coupled wave-current fractal sea surface and the model for an uncoupled fractal sea surface increase with the increase of incidence angle, as well as with increasing ocean currents. Ocean currents that are parallel to the direction of the wave can weaken the EM backscattering signal intensity, while the EM backscattering signal is intensified by ocean currents propagating oppositely to the wave direction. The model presented in this paper can be used to study the SAR imaging mechanism for a drifting sea surface. (paper)

The Calogero model can be regularised through the introduction of a cutoff parameter which removes the divergence in the interaction term. In this work we show that the one-dimensional two-particle regularised Calogero model is quasi-exactly solvable and that for certain values of the Hamiltonian parameters the eigenfunctions can be written in terms of Heun’s confluent polynomials. These eigenfunctions are such that the reduced density matrix of the two-particle density operator can be obtained exactly as well as its entanglement spectrum. We found that the number of non-zero eigenvalues of the reduced density matrix is finite in these cases. The limits for the cutoff distance going to zero (Calogero) and infinity are analysed and all the previously obtained results for the Calogero model are reproduced. Once the exact eigenfunctions are obtained, the exact von Neumann and Rényi entanglement entropies are studied to characterise the physical traits of the model. The quasi-exactly solvable character of the model is assessed studying the numerically calculated Rényi entropy and entanglement spectrum for the whole parameter space.

We develop a simplified model to better explain electric current response when direct current (DC) is applied to a flame. In particular, different current responses have been observed by changing the polarity of the DC in a sub-saturated current regime that results from the presence of ions and electrons in the flame zone. A flame zone was modeled as a thin, ionized layer located in one-dimensional DC electric fields. We derived simplified model-governing equations from species equations by implementing mobility differences dependent on the type of charged particle, particularly between ions and electrons; we performed experiments to substantiate the model. Results showed that the sub-saturated current and local field intensity were significantly influenced by the polarity of the DC because of the combined effect of unequal mobility of charged particles and the position of the ionized layer in the gap relative to two electrodes. When an energized electrode is close to the ionized layer, applying a negative DC causes a more rapid increase in current than by applying a positive DC to the same electrode. Results from our experimental measurement of current using counterflow diffusion flames agreed qualitatively well with the model predictions. A sensitivity analysis using dimensional and non-dimensional parameters also supported the importance of the mobility difference and the relative location of the ionized layer on the electric current response.

The Axelrod model for the dissemination of culture exhibits a rich spatial distribution of cultural domains, which depends on the values of the two model parameters: F, the number of cultural features and q, the common number of states each feature can assume. In the one-dimensionalmodel with F = q = 2, which is closely related to the constrained voter model, Monte Carlo simulations indicate the existence of multicultural absorbing configurations in which at least one macroscopic domain coexist with a multitude of microscopic ones in the thermodynamic limit. However, rigorous analytical results for the infinite system starting from the configuration where all cultures are equally likely show convergence to only monocultural or consensus configurations. Here we show that this disagreement is due simply to the order that the time-asymptotic limit and the thermodynamic limit are taken in the simulations. In addition, we show how the consensus-only result can be derived using Monte Carlo simulations of finite chains. (paper)

A dynamic model for PWR nuclear power plants is presented. The plant is assumed to consist of one-dimensional single-channel core, a counterflow once-through steam generator (represented by two nodes according to the nonboiling and boiling region) and the necessary connection coolant lines. The model describes analytically the frequency response behaviour of important parameters of such a plant with respect to perturbations in reactivity, subcooling or mass flow (both at the entrances to the reactor core and/or the secondary steam generator side), the perturbations in steam load or system pressure (on the secondary side of the steam generator). From corresponding 'open' loop considerations it can then be concluded - by applying the Nyquist criterion - upon the degree of the stability behaviour of the underlying system. Based on this theoretical model, a computer code named ADYPMO has been established. From the knowledge of the frequency response behaviour of such a system, the corresponding transient behaviour with respect to a stepwise or any other perturbation signal can also be calculated by applying an appropriate retransformation method, e.g. by using digital code FRETI. To demonstrate this procedure, a transient experimental curve measured during the pre-operational test period at the PWR nuclear power plant KKS Stade was recalculated using the combination ADYPMO-FRETI. Good agreement between theoretical calculations and experimental results give an insight into the validity and efficiency of the underlying theoretical model and the applied retransformation method. (Auth.)

The damage spreading (DS) transitions of two one-dimensional stochastic cellular automata suggested by Grassberger (A and B) and the kinetic Ising model of Menyhárd (NEKIM) have been investigated on the level of kinks and spins. On the level of spins the parity conservation is not satisfied and therefore studying these models provides a convenient tool to understand the dependence of DS properties on symmetries. For the model B the critical point and the DS transition point is well separated and directed percolation damage spreading transition universality was found for spin damage as well as for kink damage in spite of the conservation of damage variables modulo 2 in the latter case. For the A stochastic cellular automaton, and the NEKIM model the two transition points coincide with drastic effects on the damage of spin and kink variables showing different time dependent behaviours. While the kink DS transition is continuous and shows regular PC class universality, the spin damage exhibits a discontinuous p...

We investigate the properties of low-lying photoexcited states of a one-dimensional (1D) ionic extended Hubbard model at half-filling. Numerical analysis by using the full and Lanczos diagonalization methods shows that, in the ionic phase, there exist low-lying photoexcited states below the charge transfer gap. As a result of comparison with numerical data for the 1D antiferromagnetic (AF) Heisenberg model, it was found that, for a small alternating potential Δ, these low-lying photoexcited states are spin excitations, which is consistent with a previous analytical study [Katsura et al., link ext-link-type="uri" xlink:href="https://doi.org/10.1103/PhysRevLett.103.177402" xlink:type="simple">Phys. Rev. Lett. 103, 177402 (2009)link>]. As Δ increases, the spectral intensity of the 1D ionic extended Hubbard model rapidly deviates from that of the 1D AF Heisenberg model and it is clarified that this deviation is due to the neutral-ionic domain wall, an elementary excitation near the neutral-ionic transition point.

A one-dimensional homogenized model for dynamic fluid-structure interaction in a pressurized water reactor core is used to study the influence of the virtual density and spacer's stiffness. The model consists of a linear system of partial differential equations for fluid velocity, rod velocity and pressure. For these equations analytical solutions are deduced for boundary conditions prescribing either periodic wall oscillations or linearly growing wall accelerations from rest. The theoretical model for the virtual density is verified by comparison to an experiment. For zero spacer stiffness, purely acoustic oscillations appear. For positive spacer stiffness, additional oscillations arise with relative rod motions. The wavelengths of the latter oscillations are small for weak spacers. Large numerical effort would be required in a more complete three-dimensional core-model to resolve such short wave lengths. In fact in a typical core the spacer's stiffness csub(S) is small in comparison to the fluid bulk modulus K. For csub(s)/K <= 0.1 it might be appropriate to neglect the influence of the spacers. (orig.)

We study the interplay between the electron-phonon (e -ph) and on-site electron-electron (e-e) interactions in a three-orbital Hubbard-Holstein model on an extended one-dimensional lattice using determinant quantum Monte Carlo. For weak e-e and e -ph interactions, we observe a competition between an orbital-selective Mott phase (OSMP) and a (multicomponent) charge-density-wave (CDW) insulating phase, with an intermediate metallic phase located between them. For large e-e and e -ph couplings, the OSMP and CDW phases persist, while the metallic phase develops short-range orbital correlations and becomes insulating when both the e-e and e -ph interactions are large but comparable. Many of our conclusions are in line with those drawn from a prior dynamical mean-field theory study of the two-orbital Hubbard-Holstein model [Phys. Rev. B 95, 121112(R) (2017), 10.1103/PhysRevB.95.121112] in infinite dimension, suggesting that the competition between the e -ph and e-e interactions in multiorbital Hubbard-Holstein models leads to rich physics, regardless of the dimension of the system.

Ferrogels are smart soft materials, consisting of a polymeric network and embedded magnetic particles. Novel phenomena, such as the variation of the overall mechanical properties by external magnetic fields, emerge consequently. However, the dynamic behavior of ferrogels remains largely unveiled. In this paper, we consider a one-dimensional chain consisting of magnetic dipoles and elastic springs between them as a simple model for ferrogels. The model is evaluated by corresponding simulations. To probe the dynamics theoretically, we investigate a continuum limit of the energy governing the system and the corresponding equation of motion. We provide general classification scenarios for the dynamics, elucidating the touching/detachment dynamics of the magnetic particles along the chain. In particular, it is verified in certain cases that the long-time relaxation corresponds to solutions of shock-wave propagation, while formations of particle pairs underlie the initial stage of the dynamics. We expect that these results will provide insight into the understanding of the dynamics of more realistic models with randomness in parameters and time-dependent magnetic fields.

Estimating the magnitude and the intensity of rapid landslides like debris flows is fundamental to evaluate quantitatively the hazard in a specific location. Intensity varies through the travelled course of the flow and can be described by physical features such as deposited volume, velocities, height of the flow, impact forces and pressures. Dynamic run-out models are able to characterize the distribution of the material, its intensity and define the zone where the elements will experience an impact. These models can provide valuable inputs for vulnerability and risk calculations. However, most dynamic run-out models assume a constant volume during the motion of the flow, ignoring the important role of material entrained along its path. Consequently, they neglect that the increase of volume enhances the mobility of the flow and can significantly influence the size of the potential impact area. An appropriate erosion mechanism needs to be established in the analyses of debris flows that will improve the results of dynamic modeling and consequently the quantitative evaluation of risk. The objective is to present and test a simple 1D debris flow model with a material entrainment concept based on limit equilibrium considerations and the generation of excess pore water pressure through undrained loading of the in situ bed material. The debris flow propagation model is based on a onedimensional finite difference solution of a depth-averaged form of the Navier-Stokes equations of fluid motions. The flow is treated as a laminar one phase material, which behavior is controlled by a visco-plastic Coulomb-Bingham rheology. The model parameters are evaluated and the model performance is tested on a debris flow event that occurred in 2003 in the Faucon torrent (Southern French Alps).

Generally, one observes an anomaly of temperature before a big earthquake. In this paper, we established the expression of thermal energy produced by friction forces between the walls of a seismic fault while considering the dynamic of a one-dimensional spring-block model. It is noted that, before the rupture of a seismic fault, displacements are caused by microseisms. The curves of variation of this thermal energy with time show that, for oscillatory and aperiodic displacement, the thermal energy is accumulated in the same way. The study reveals that thermal energy as well as temperature increases abruptly after a certain amount of time. We suggest that the corresponding time is the start of the anomaly of temperature observed which can be considered as precursory effect of a big seism. We suggest that the thermal energy can heat gases and dilate rocks until they crack. The warm gases can then pass through the cracks towards the surface. The cracks created by thermal energy can also contribute to the rupture of the seismic fault. We also suggest that the theoretical model of thermal energy, produced in seismic fault, associated with a large quantity of experimental data may help in the prediction of earthquakes.

The principles of the δ perturbation theory were first proposed in the context of self-interacting scalar quantum field theory. There it was shown how to expand a (φ 2 ) 1+δ theory as a series in powers of δ and how to recover nonperturbative information about a φ 4 field theory from the δ expansion at δ=1. The purpose of this series of papers is to extend the notions of δ perturbation theory from boson theories to theories having a local gauge symmetry. In the case of quantum electrodynamics one introduces the parameter δ by generalizing the minimal coupling terms to bar ψ(∂-ieA) δ ψ and expanding in powers of δ. This interaction preserves local gauge invariance for all δ. While there are enormous benefits in using the δ expansion (obtaining nonperturbative results), gauge theories present new technical difficulties not encountered in self-interacting boson theories because the expression (∂-ieA) δ contains a derivative operator. In the first paper of this series a one-dimensionalmodel whose interaction term has the form bar ψ[d/dt-igφ(t)] δ ψ is considered. The virtue of this model is that it provides a laboratory in which to study fractional powers of derivative operators without the added complexity of γ matrices. In the next paper of this series we consider two-dimensional electrodynamics and show how to calculate the anomaly in the δ expansion

In this paper we determine the nonequilibrium magnetic work performed on a Ising model and relate it to the fluctuation theorem derived some years ago by Jarzynski. The basic idea behind this theorem is the relationship connecting the free energy difference between two thermodynamic states of a system and the average work performed by an external agent, in a finite time, through nonequilibrium paths between the same thermodynamic states. We test the validity of this theorem by considering the one-dimensional Ising model where the free energy is exactly determined as a function of temperature and magnetic field. We have found that the Jarzynski theorem remains valid for all the values of the rate of variation of the magnetic field applied to the system. We have also determined the probability distribution function for the work performed on the system for the forward and reverse processes and verified that predictions based on the Crooks relation are equally correct. We also propose a method to calculate the lag between the current state of the system and that of the equilibrium based on macroscopic variables. We have shown that the lag increases with the sweeping rate of the field at its final value for the reverse process, while it decreases in the case of the forward process. The lag increases linearly with the size of the chain and with a slope decreasing with the inverse of the rate of variation of the field.

A one-dimensional (1D) model of oxygen transport in the diffusion media of proton exchange membrane fuel cells (PEMFC) is presented, which considers convection perpendicular to the electrode in addition to diffusion. The resulting analytical expression of the convecto-diffusive impedance is obtained using a convection-diffusion equation instead of a diffusion equation in the case of classical Warburg impedance. The main hypothesis of the model is that the convective flux is generated by the evacuation of water produced at the cathode which flows through the porous media in vapor phase. This allows the expression of the convective flux velocity as a function of the current density and of the water transport coefficient {alpha} (the fraction of water being evacuated at the cathode outlet). The resulting 1D oxygen transport impedance neglects processes occurring in the direction parallel to the electrode that could have a significant impact on the cell impedance, like gas consumption or concentration oscillations induced by the measuring signal. However, it enables us to estimate the impact of convection perpendicular to the electrode on PEMFC impedance spectra and to determine in which conditions the approximation of a purely diffusive oxygen transport is valid. Experimental observations confirm the numerical results. (Copyright copyright 2012 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim)

A field investigation characterizing contamination at the Rocky Flats Plant (Rocky Flats Environmental Technology Site) near Golden, Colorado revealed unexpectedly high moisture contents in the unsaturated soil column (vadose zone) beneath several of the Plant's Waste Water Treatment Plant (WWTP) sludge drying beds. Because these beds were seldom in use, researchers had hypothesized that the water required to maintain the saturated conditions observed beneath several of the sludge drying beds was coming from sources other than the beds themselves. In an effort to substantiate this hypothesis, a one-dimensional physically-based unsaturated flow model was utilized to simulate the vertical movement of moisture from the sludge drying beds into the unsaturated soil column below. The model was run to simulate vertical flow over a two-year period and results indicated that no significant changes from initial conditions were apparent. This evidence supports the hypothesis that the high moisture contents found beneath the sludge drying beds are being fed by sources other than infiltration of sludge applied to the beds themselves. This paper presents the details of the simulation and provides further evidence of the hypothesized flow regime

This thesis presents the one-dimensional equations, numerical method and simulations of a model to characterize the dynamical operation of an electrochemical cell. This model extends the current state-of-the art in that it accounts, in a primitive way, for the physics of the electrolyte/electrode interface and incorporates diffuse-charge dynamics, temperature coupling, surface coverage, and polarization phenomena. The one-dimensional equations account for a system with one or two mobile ions of opposite charge, and the electrode reaction we consider (when one is needed) is a one-electron electrodeposition reaction. Though the modeled system is far from representing a realistic electrochemical device, our results show a range of dynamics and behaviors which have not been observed previously, and explore the numerical challenges required when adding more complexity to a model. Furthermore, the basic transport equations (which are developed in three spatial dimensions) can in future accomodate the inclusion of additional physics, and coupling to more complex boundary conditions that incorporate two-dimensional surface phenomena and multi-rate reactions. In the model, the Poisson-Nernst-Planck equations are used to model diffusion and electromigration in an electrolyte, and the generalized Frumkin-Butler-Volmer equation is used to model reaction kinetics at electrodes. An energy balance equation is derived and coupled to the diffusion-migration equation. The model also includes dielectric polarization effects by introducing different values of the dielectric permittivity in different regions of the bulk, as well as accounting for surface coverage effects due to adsorption, and finite size "crowding", or steric effects. Advection effects are not modeled but could in future be incorporated. In order to solve the coupled PDE's, we use a variable step size second order scheme in time and finite differencing in space. Numerical tests are performed on a simplified system and

A method for allocating waste package performance requirements among waste package components with regard to radionuclide isolation has been developed. Modification or change in this approach can be expected as the understanding of radionuclide behavior in the waste package improves. Thus, the performance requirements derived in this document are preliminary and subject to change. However, this kind of analysis is a useful starting point. It has also proved useful for identifying a small group of radionuclides which should be emphasized in a laboratory experimental program designed to characterize the behavior of specific radionuclides in the waste package environment. A simple one-dimensional, two media transport model has been derived and used to calculate radionuclide transport from the waste form-packing material interface of the waste package into the host rock. Cumulative release over 10,000 years, maximum yearly releases and release rates at the packing material-host rock interface were evaluated on a radionuclide-by radionuclide basis. The major parameters controlling radionuclide release were found to be: radionuclide solubility, porosity of the rock, isotopic ratio of the radionuclide and surface area of the waste form-packing material interface. 15 refs., 2 figs., 16 tabs

We discuss the nonequilibrium dynamics in two interaction quantum quenches in the one-dimensional Fermi-Hubbard model. First, we study the decay of the Néel state as a function of interaction strength. We observe a fast charge dynamics over which double occupancies are built up, while the long-time decay of the staggered moment is controlled by spin excitations, corroborated by the analysis of the entanglement dynamics. Second, we investigate the formation of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) correlations in a spin-imbalanced system in quenches from the noninteracting case to attractive interactions. Even though the quench puts the system at a finite energy density, peaks at the characteristic FFLO quasimomenta are visible in the quasi-momentum distribution function, albeit with an exponential decay of s-wave pairing correlations. We also discuss the imprinting of FFLO correlations onto repulsively bound pairs and their rapid decay in ramps. Supported by the DFG (Deutsche Forschungsgemeinschaft) via FOR 1807.

Full Text Available The one-dimensional two-fluid model approach has been traditionally used in thermal-hydraulics codes for the analysis of transients and accidents in water–cooled nuclear power plants. This paper investigates the performance of RELAP5/MOD3 predicting vertical upward bubbly flow at low velocity conditions. For bubbly flow and vertical pipes, this code applies the drift-velocity approach, showing important discrepancies with the experiments compared. Then, we use a classical formulation of the drag coefficient approach to evaluate the performance of both approaches. This is based on the critical Weber criteria and includes several assumptions for the calculation of the interfacial area and bubble size that are evaluated in this work. A more accurate drag coefficient approach is proposed and implemented in RELAP5/MOD3. Instead of using the Weber criteria, the bubble size distribution is directly considered. This allows the calculation of the interfacial area directly from the definition of Sauter mean diameter of a distribution. The results show that only the proposed approach was able to predict all the flow characteristics, in particular the bubble size and interfacial area concentration. Finally, the computational results are analyzed and validated with cross-section area average measurements of void fraction, dispersed phase velocity, bubble size, and interfacial area concentration.

Some of the previous studies have shown that the use of compressed natural gas (CNG) in diesel engines potentially produce engine performance improvement and exhaust gas emission reduction, especially nitrogen oxides, unburned hydrocarbons, and carbon dioxide. On the other hand, there are other researchers who claimed that the use of CNG increases exhaust gas emissions, particularly nitrogen oxides. In this study, a one-dimensionalmodel of a diesel-CNG dual fuel engine was made based on a 4-cylinder 2.5L common rail direct injection diesel engine. The software used is GT-Power, and it was used to analyze the engine performance and exhaust gas emissions of several diesel-CNG dual fuel blend ratios, i.e. 100:0, 90:10, 80:20, 70:30, 60:40 and 50:50. The effect of 100%, 75%, 50% engine loads on the exhaust gas emissions were also studied. The result shows that all diesel-CNG fuel blends produces higher brake torque and brake power at engine speed of 2000-3000 rpm compared with 100% diesel. The 50:50 diesel-CNG blend produces the highest brake torque and brake power, but also has the highest brake specific fuel consumption. As a higher percentage of CNG added to the dual fuel blend, unburned hydrocarbons and carbon monoxide emission increased while carbon dioxide emission decreased. The nitrogen oxides emission concentration is generally unaffected by any change of the dual fuel ratio.

A one-dimensionalmodel that solves the time-dependent equations for the zonal mean wind and a wave of specified zonal wavenumber has been used to illustrate the ability of gravity waves forced by time-dependent tropospheric heating to produce a semiannual oscillation (SAO) in the middle atmosphere. When the heating has a strong diurnal cycle, as observed over tropical landmasses, gravity waves with zonal wavelengths of a few thousand kilometers and phase velocities in the range +/- 40-50 m/sec are excited efficiently by the maximum vertical projection criterion (vertical wavelength approximately equals 2 x forcing depth). Calculations show that these waves can account for large zonal mean wind accelerations in the middle atmosphere, resulting in realistic stratopause and mesopause oscillations. Calculations of the temporal evolution of a quasi-conserved tracer indicate strong down-welling in the upper stratosphere near the equinoxes, which is associated with the descent of the SAO westerlies. In the upper mesosphere, there is a semiannual oscillation in tracer mixing ratio driven by seasonal variability in eddy mixing, which increases at the solstices and decreases at the equinoxes.

We study the real-time evolution of a charge-excited state in a one-dimensional e{sub g}-orbital degenerate Hubbard model, by a time-dependent density-matrix renormalization group method. Considering a chain along the z direction, electrons hop between adjacent 3z{sup 2}-r{sup 2} orbitals, while x{sup 2}-y{sup 2} orbitals are localized. For the charge-excited state, a holon-doublon pair is introduced into the ground state at quarter filling. At initial time, there is no electron in a holon site, while a pair of electrons occupies 3z{sup 2}-r{sup 2} orbital in a doublon site. As the time evolves, the holon motion is governed by the nearest-neighbor hopping, but the electron pair can transfer between 3z{sup 2}-r{sup 2} orbital and x{sup 2}-y{sup 2} orbital through the pair hopping in addition to the nearest-neighbor hopping. Thus holon and doublon propagate at different speed due to the pair hopping that is characteristic of multi-orbital systems.

In the stratosphere the ozone layer is very sensitive to the NOx abundance. The ionisation of N2 and O2 molecules by TLE's (Transient Luminous Events) is a source of NOx which is currently not well quantified and could act as a loss of ozone. In this study a onedimensional explicit parameterization of a Blue-Jet propagation based on that proposed by Raizer et al. (2006 and 2007) has been developed. This parameterization considers Blue-Jet as a streamer initiated by a bidirectional leader discharge, emerging from the anvil and sustained by moderate cloud charge. The streamer growth varies with the electrical field induced by initial cloud charge and the initial altitude. This electrical parameterization and the chemical mechanisms associated with the discharge have been implemented into a detailed chemical model of stratospheric ozone including evolution of nitrogen, chlorine and bromine species. We will present several tests performed to validate the electrical code and evaluate the propagation velocity and the maximum altitude attains by the blue jet as a function of electrical parameters. The results obtained giving the spatiotemporal evolution of the electron density are then used to initiate the specific chemistry associated with the Blue Jet. Preliminary results on the impact of such discharge on the ozone content and the whole stratospheric system will be presented.

It is ideal to use exact equilibrium solutions of the steady state Vlasov-Maxwell system to intialise collsionless simulations. However, exact equilibrium distribution functions (DFs) for a given macroscopic configuration are typically unknown, and it is common to resort to using `flow-shifted' Maxwellian DFs in their stead. These DFs may be consistent with a macrosopic system with the target number density and current density, but could well have inaccurate higher order moments. We present recent theoretical work on the `inverse problem in Vlasov-Maxwell equilibria', namely calculating an exact solution of the Vlasov equation for a specific given magnetic field. In particular, we focus on one-dimensional geometries in Cartesian (current sheets) coordinates.1. From 1D fields to Vlasov equilibria: Theory and application of Hermite Polynomials: (O. Allanson, T. Neukirch, S. Troscheit and F. Wilson, Journal of Plasma Physics, 82, 905820306 (2016) [28 pages, Open Access] )2. An exact collisionless equilibrium for the Force-Free Harris Sheet with low plasma beta: (O. Allanson, T. Neukirch, F. Wilson and S. Troscheit, Physics of Plasmas, 22, 102116 (2015) [11 pages, Open Access])3. Neutral and non-neutral collisionless plasma equilibria for twisted flux tubes: The Gold-Hoyle model in a background field (O. Allanson, F. Wilson and T. Neukirch, (2016)) (accepted, Physics of Plasmas)

The hypothesis (Sparenberg et al. in EPJ Web Conf 58:01016, [1]. https://doi.org/10.1051/epjconf/20135801016) that the particular linear tracks appearing in the measurement of a spherically-emitting radioactive source in a cloud chamber are determined by the (random) positions of atoms or molecules inside the chamber is further explored in the framework of a recently established one-dimensionalmodel (Carlone et al. Comm Comput Phys 18:247, [2]. https://doi.org/10.4208/cicp.270814.311214a). In this model, meshes of localized spins 1/2 play the role of the cloud-chamber atoms and the spherical wave is replaced by a linear superposition of two wave packets moving from the origin to the left and to the right, evolving deterministically according to the Schrödinger equation. We first revisit these results using a time-dependent approach, where the wave packets impinge on a symmetric two-sided detector. We discuss the evolution of the wave function in the configuration space and stress the interest of a non-symmetric detector in a quantum-measurement perspective. Next we use a time-independent approach to study the scattering of a plane wave on a single-sided detector. Preliminary results are obtained, analytically for the single-spin case and numerically for up to 8 spins. They show that the spin-excitation probabilities are sometimes very sensitive to the parameters of the model, which corroborates the idea that the measurement result could be determined by the atom positions. The possible origin of decoherence and entropy increase in future models is finally discussed.

The suppression pool in a boiling water reactor (BWR) plant not only is the major heat sink within the containment system, but also provides the major emergency cooling water for the reactor core. In several accident scenarios, such as a loss-of-coolant accident and extended station blackout, thermal stratification tends to form in the pool after the initial rapid venting stage. Accurately predicting the pool stratification phenomenon is important because it affects the peak containment pressure; the pool temperature distribution also affects the NPSHa (available net positive suction head) and therefore the performance of the Emergency Core Cooling System and Reactor Core Isolation Cooling System pumps that draw cooling water back to the core. Current safety analysis codes use zero dimensional (0-D) lumped parameter models to calculate the energy and mass balance in the pool; therefore, they have large uncertainties in the prediction of scenarios in which stratification and mixing are important. While three-dimensional (3-D) computational fluid dynamics (CFD) methods can be used to analyze realistic 3-D configurations, these methods normally require very fine grid resolution to resolve thin substructures such as jets and wall boundaries, resulting in a long simulation time. For mixing in stably stratified large enclosures, the BMIX++ code (Berkeley mechanistic MIXing code in C++) has been developed to implement a highly efficient analysis method for stratification where the ambient fluid volume is represented by one-dimensional (1-D) transient partial differential equations and substructures (such as free or wall jets) are modeled with 1-D integral models. This allows very large reductions in computational effort compared to multi-dimensional CFD modeling. One heat-up experiment performed at the Finland POOLEX facility, which was designed to study phenomena relevant to Nordic design BWR suppression pool including thermal stratification and mixing, is used for

Full Text Available Frost formation is a common, often undesired phenomenon in heat exchanges such as air coolers. Thus, air coolers have to be defrosted periodically, causing significant energy consumption. For the design and optimization, prediction of defrosting by a CFD tool is desired. This paper presents a one-dimensional transient model approach suitable to be used as a zero-dimensional wall-function in CFD for modeling the defrost process at the fin and tube interfaces. In accordance to previous work a multi stage defrost model is introduced (e.g. [1, 2]. In the first instance the multi stage model is implemented and validated using MATLAB. The defrost process of a one-dimensional frost segment is investigated. Fixed boundary conditions are provided at the frost interfaces. The simulation results verify the plausibility of the designed model. The evaluation of the simulated defrost process shows the expected convergent behavior of the three-stage sequence.

Photonic crystals and their use in exciting Bloch surface waves have received immense attention over the past few decades. This interest is mainly due to their applications in bio-sensing, wave-guiding, and other optical phenomena such as surface field enhanced Raman spectroscopy. Improvement in numerical modeling techniques, state of the art computing resources, and advances in fabrication techniques have also assisted in growing interest in this field. The ability to model photonic crystals computationally has benefited both the theoretical as well as experimental communities. It helps the theoretical physicists in solving complex problems which cannot be solved analytically and helps to acquire useful insights that cannot be obtained otherwise. Experimentalists, on the other hand, can test different variants of their devices by changing device parameters to optimize performance before fabrication. In this dissertation, we develop two commonly used numerical techniques, namely transfer matrix method, and rigorous coupled wave analysis, in C++ and MATLAB, and use two additional software packages, one open-source and another commercial, to modelone-dimensional photonic crystals. Different variants of one-dimensional multilayered structures such as perfectly periodic dielectric multilayers, quasicrystals, aperiodic multilayer are modeled, along with one-dimensional photonic crystals with gratings on the top layer. Applications of Bloch surface waves, along with new and novel aperiodic dielectric multilayer structures that support Bloch surface waves are explored in this dissertation. We demonstrate a slow light configuration that makes use of Bloch Surface Waves as an intermediate excitation in a double-prism tunneling configuration. This method is simple compared to the more usual techniques for slowing light using the phenomenon of electromagnetically induced transparency in atomic gases or doped ionic crystals operated at temperatures below 4K. Using a semi

Electronic-behavior modeling of three-dimensional (3D) p + -π-p + and n + -ν-n + semiconducting diffusible devices with highly accurate resistances for the design of analog resistors, which are compatible with the CMOS (complementary-metal-oxide-semiconductor) technologies, is performed in three dimensions with the fast tube multiplexing method (TMM). The current–voltage (I–V) curve of a silicon device is usually computed with traditional device simulators of technology computer-aided design (TCAD) based on the finite-element method (FEM). However, for the design of 3D p + -π-p + and n + -ν-n + diffusible resistors, they show a high computational cost and convergence that may fail with fully non-separable 3D dopant concentration profiles as observed in many diffusible resistors resulting from laser trimming. These problems are avoided with the proposed TMM, which divides the 3D resistor into one-dimensional (1D) thin tubes with longitudinal axes following the main orientation of the average electrical field in the tubes. The I–V curve is rapidly obtained for a device with a realistic 3D dopant profile, since a system of three first-order ordinary differential equations has to be solved for each 1D multiplexed tube with the TMM instead of three second-order partial differential equations in the traditional TCADs. Simulations with the TMM are successfully compared to experimental results from silicon-based 3D resistors fabricated by laser-induced dopant diffusion in the gaps of MOSFETs (metal-oxide-semiconductor field-effect transistors) without initial gate. Using thin tubes with other shapes than parallelepipeds as ring segments with toroidal lateral surfaces, the TMM can be generalized to electronic devices with other types of 3D diffusible microstructures

This thesis contributes to the field of strongly correlated electron systems with studies in two distinct fields thereof: the specific nature of correlations between electrons in one dimension and quantum quenches in quantum impurity problems. In general, strongly correlated systems are characterized in that their physical behaviour needs to be described in terms of a many-body description, i.e. interactions correlate all particles in a complex way. The challenge is that the Hilbert space in a many-body theory is exponentially large in the number of particles. Thus, when no analytic solution is available - which is typically the case - it is necessary to find a way to somehow circumvent the problem of such huge Hilbert spaces. Therefore, the connection between the two studies comes from our numerical treatment: they are tackled by the density matrix renormalization group (DMRG) and the numerical renormalization group (NRG), respectively, both based on matrix product states. The first project presented in this thesis addresses the problem of numerically finding the dominant correlations in quantum lattice models in an unbiased way, i.e. without using prior knowledge of the model at hand. A useful concept for this task is the correlation density matrix (CDM) which contains all correlations between two clusters of lattice sites. We show how to extract from the CDM, a survey of the relative strengths of the system's correlations in different symmetry sectors as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. We demonstrate this by a DMRG study of a one-dimensional spinless extended Hubbard model, while emphasizing that the proposed analysis of the CDM is not restricted to one dimension. The second project presented in this thesis is motivated by two phenomena under ongoing experimental and theoretical investigation in the context of quantum impurity models: optical absorption

This thesis contributes to the field of strongly correlated electron systems with studies in two distinct fields thereof: the specific nature of correlations between electrons in one dimension and quantum quenches in quantum impurity problems. In general, strongly correlated systems are characterized in that their physical behaviour needs to be described in terms of a many-body description, i.e. interactions correlate all particles in a complex way. The challenge is that the Hilbert space in a many-body theory is exponentially large in the number of particles. Thus, when no analytic solution is available - which is typically the case - it is necessary to find a way to somehow circumvent the problem of such huge Hilbert spaces. Therefore, the connection between the two studies comes from our numerical treatment: they are tackled by the density matrix renormalization group (DMRG) and the numerical renormalization group (NRG), respectively, both based on matrix product states. The first project presented in this thesis addresses the problem of numerically finding the dominant correlations in quantum lattice models in an unbiased way, i.e. without using prior knowledge of the model at hand. A useful concept for this task is the correlation density matrix (CDM) which contains all correlations between two clusters of lattice sites. We show how to extract from the CDM, a survey of the relative strengths of the system's correlations in different symmetry sectors as well as detailed information on the operators carrying long-range correlations and the spatial dependence of their correlation functions. We demonstrate this by a DMRG study of a one-dimensional spinless extended Hubbard model, while emphasizing that the proposed analysis of the CDM is not restricted to one dimension. The second project presented in this thesis is motivated by two phenomena under ongoing experimental and theoretical investigation in the context of quantum impurity models: optical absorption

Highlights: • A non-equilibrium thermal mode with considering loses is adopted in Stirling engine. • Good agreements are achieved for predicting various critical system parameters. • Differences between helium and hydrogen systems are highlighted and analyzed. • Pressure drop of helium system is much larger and more sensitive to frequency. - Abstract: A third-order numerical model based on one-dimensional computational fluid dynamics is developed for kinetic Stirling engines. Various loss mechanisms in Stirling engines, including gas spring hysteresis loss, shuttle loss, appendix displacer gap loss, gas leakage loss, finite speed loss, piston friction loss, pressure drop loss, heat conduction loss, mechanical loss and imperfect heat transfer, are considered and embedded into the basic control equations. The non-equilibrium thermal model is adopted for the regenerator to capture the oscillating features of the gas and solid temperatures. To improve the numerical stability and accuracy, the implicit second-order time difference scheme and the second-order upwind scheme are adopted for discretizing the time differential terms and convective terms, respectively. Experimental validations are then conducted on a beta-type Stirling engine with the extensive experimental data for diverse working conditions. The results show that the developed model has better accuracies than the previous second-order models. Good agreements are achieved for predicting various critical system parameters, including pressure-volume diagram, indicated power, brake power, indicated efficiency, brake efficiency and mechanical efficiency. In particular, both the experiments and simulations show that the Stirling engine charged with helium tends to have much lower optimal working frequencies and poorer performances compared to the hydrogen system. Based on the analyses of the losses, it reveals that the pressure drop in the flow channels plays a critical role in shaping the different

Evaluation of pluvial flood risk is often based on computations using 1D/2D urban flood models. However, guidelines on choice of model complexity are missing, especially for one-dimensional (1D) network models. This study presents a new automatic approach for simplification of 1D hydraulic networ...

One-dimensionalmodel for vertical profiles of longitudinal velocities in open-channel flows is verified against laboratory data obtained in an open channel with artificial plants. Those plants simulate Canadian waterweed which in nature usually forms dense stands that reach all the way to the water surface. The model works particularly well for densely spaced plants.

Based on the ground states of the one-dimensional Lipkin-Meshkov-Glick model (LMGM), we show an all-versus-nothing proof of violation of local realism in this model. Moreover, the quantum entanglement swapping is also investigated in terms of the braiding transformations. (general)

An assumption of quasi-one-dimensionality is used to derive a simple set of equations describing the comet Giacobini-Zinner tail configuration. The MHD equations are expanded in terms of a parameter representing the ratio of the length scale in the direction perpendicular to the neutral sheet over the length scale in the direction parallel to the tail. It is shown that in this way it is possible to obtain much information on the structure of the tail and to fit reasonably well the observations made by the ICE spacecraft

Full Text Available One-dimensional dynamic design of a component characterized by inertia coefficient, elastic coefficient, and coefficient of energy dispersion. The component is affected by external action in the form of time-independent initial kinematic disturbances and varying ones. Mathematical model of component dynamics as well as a new form of analytical representation of transient in terms of one-dimensional problem of kinematic effect is provided. Dynamic design of a component is being carried out according to a theory of modal control.

Full Text Available The laboratory methods used for the soil water retention curve (SWRC construction and parameter estimation is time-consuming. A vertical infiltration method was proposed to estimate parameters α and n and to further construct the SWRC. In the present study, the relationships describing the cumulative infiltration and infiltration rate with the depth of the wetting front were established, and simplified expressions for estimating α and n parameters were proposed. The one-dimensional vertical infiltration experiments of four soils were conducted to verify if the proposed method would accurately estimate α and n. The fitted values of α and n, obtained from the RETC software, were consistent with the calculated values obtained from the infiltration method. The comparison between the measured SWRCs obtained from the centrifuge method and the calculated SWRCs that were based on the infiltration method displayed small values of root mean square error (RMSE, mean absolute percentage error (MAPE, and mean absolute error. SWMS_2D-based simulations of cumulative infiltration, based on the calculated α and n, remained consistent with the measured values due to small RMSE and MAPE values. The experiments verified the proposed one-dimensional vertical infiltration method, which has applications in field hydraulic parameter estimation.

Full Text Available A one-dimensional electrostatic sheet model of a coaxial geometry Virtual Cathode Oscillator (VCO) is presented. The cathode is centrally located and connects to a peripherally located plate electrode to form a resonant cavity, and is thus...

This investigation is to show that two-dimensional steady state heat transfer problems of composite walls should not be solved by the conventionally one-dimensional parallel thermal resistance circuits (PTRC) model because the interface temperatures are not unique. Thus PTRC model cannot be used like its conventional recognized analogy, parallel electrical resistance circuits (PERC) model which has the unique node electric voltage. Two typical composite wall examples, solved by CFD software, are used to demonstrate the incorrectness. The numerical results are compared with those obtained by PTRC model, and very large differences are observed between their results. This proves that the application of conventional heat transfer PTRC model to two-dimensional composite walls, introduced in most heat transfer text book, is totally incorrect. An alternative one-dimensional separately series thermal resistance circuit (SSTRC) model is proposed and applied to the two-dimensional composite walls with isothermal boundaries. Results with acceptable accuracy can be obtained by the new model

This study is to prove that two-dimensional steady state heat transfer problems of composite circular pipes cannot be appropriately solved by the conventional one-dimensional parallel thermal resistance circuits (PTRC) model because its interface temperatures are not unique. Thus, the PTRC model is definitely different from its conventional recognized analogy, parallel electrical resistance circuits (PERC) model, which has unique node electric voltages. Two typical composite circular pipe examples are solved by CFD software, and the numerical results are compared with those obtained by the PTRC model. This shows that the PTRC model generates large error. Thus, this conventional model, introduced in most heat transfer text books, cannot be applied to two-dimensional composite circular pipes. On the contrary, an alternative one-dimensional separately series thermal resistance circuit (SSTRC) model is proposed and applied to a two-dimensional composite circular pipe with isothermal boundaries, and acceptable results are returned

A GIS-based one-dimensional flood simulation model is presented and applied to the centre of the city of Nîmes (Gard, France), for mapping flow depths or velocities in the streets network. The geometry of the one-dimensional elements is derived from the Digital Elevation Model (DEM). The flow is routed from one element to the next using the kinematic wave approximation. At the crossroads, the flows in the downstream branches are computed using a conceptual scheme. This scheme was previously designed to fit Y-shaped pipes junctions, and has been modified here to fit X-shaped crossroads. The results were compared with the results of a two-dimensional hydrodynamic model based on the full shallow water equations. The comparison shows that good agreements can be found in the steepest streets of the study zone, but differences may be important in the other streets. Some reasons that can explain the differences between the two models are given and some research possibilities are proposed.

Localization properties of noninteracting quantum particles in one-dimensional incommensurate lattices are investigated with an exponential short-range hopping that is beyond the minimal nearest-neighbor tight-binding model. Energy dependent mobility edges are analytically predicted in this model and verified with numerical calculations. The results are then mapped to the continuum Schroedinger equation, and an approximate analytical expression for the localization phase diagram and the energy dependent mobility edges in the ground band is obtained.

Aluminum oxide slag accumulation and expulsion was previously shown to be a player in various solid rocket motor phenomena, including the Space Shuttle's Reusable Solid Rocket Motor (RSRM) pressure perturbation, or "blip," and phantom moment. In the latter case, such un ]commanded side accelerations near the end of burn have also been identified in several other motor systems. However, efforts to estimate the mass expelled during a given event have come up short. Either bulk calculations are performed without enough physics present, or multiphase, multidimensional Computational Fluid Dynamic analyses are performed that give a snapshot in time and space but do not always aid in grasping the general principle. One ]dimensional, two ]phase compressible flow calculations yield an analytical result for nozzle flow under certain assumptions. This can be carried further to relate the bulk motor parameters of pressure, thrust, and mass flow rate under the different exhaust conditions driven by the addition of condensed phase mass flow. An unknown parameter is correlated to airflow testing with water injection where mass flow rates and pressure are known. Comparison is also made to full ]scale static test motor data where thrust and pressure changes are known and similar behavior is shown. The end goal is to be able to include the accumulation and flow of slag in internal ballistics predictions. This will allow better prediction of the tailoff when much slag is ejected and of mass retained versus time, believed to be a contributor to the widely-observed "flight knockdown" parameter.

When an untapered high-gain free electron laser (FEL) reaches saturation, the exponential growth ceases and the radiation power starts to oscillate about an equilibrium. The FEL radiation power or efficiency can be increased by undulator tapering. For a high-gain tapered FEL, although the power is enhanced after the first saturation, it is known that there is a so-called second saturation where the FEL power growth stops even with a tapered undulator system. The sideband instability is one of the primary reasons leading to this second saturation. In this paper, we provide a quantitative analysis on how the gradient of undulator tapering can mitigate the sideband growth. The study is carried out semianalytically and compared with one-dimensional numerical simulations. The physical parameters are taken from Linac Coherent Light Source-like electron bunch and undulator systems. The sideband field gain and the evolution of the radiation spectra for different gradients of undulator tapering are examined. It is found that a strong undulator tapering (˜10 %) provides effective suppression of the sideband instability in the postsaturation regime.

A general one-dimensional (1-D), three-region model for a bubbling fluidized-bed adsorber with internal heat exchangers has been developed. The model can predict the hydrodynamics of the bed and provides axial profiles for all temperatures, concentrations, and velocities. The model is computationally fast and flexible and allows for any system of adsorption and desorption reactions to be modeled, making the model applicable to any adsorption process. The model has been implemented in both gPROMS and Aspen Custom Modeler, and the behavior of the model has been verified.

Full Text Available We define a class of Y(sl_{(m|n} Yangian invariant Haldane-Shastry (HS like spin chains, by assuming that their partition functions can be written in a particular form in terms of the super Schur polynomials. Using some properties of the super Schur polynomials, we show that the partition functions of this class of spin chains are equivalent to the partition functions of a class of one-dimensional vertex models with appropriately defined energy functions. We also establish a boson-fermion duality relation for the partition functions of this class of supersymmetric HS like spin chains by using their correspondence with one-dimensional vertex models.

conditions. As the pollutant load on the estuary increases, the. water quality may deteriorate rapidly and therefore the scientific interests are centered on the analysis of water quality. The pollutants will be subjected to a number of physical, chemical... study we have applied one-dimensional advection-diffusion model for the waters of Gauthami Godavari estuary to determine the axial diffusion coefficients and thereby to predict the impact assessment. The study area (Fig. 1) is the lower most 32 km...

The long Middle Route of the South to North Water Transfer Project is composed of complex hydraulic structures (aqueduct, tunnel, control gate, diversion, culvert, and diverted siphon), which generate complex flow patterns. It is vital to simulate the flow patterns through hydraulic structures, but it is a challenging work to protect water quality and maintain continuous water transfer. A one-dimensional hydrodynamic and water quality model was built to understand the flow and pollutant movem...

We present some estimates for the memory kernel function in the generalized Langevin equation, derived using the Mori-Zwanzig formalism from a one-dimensional lattice model, in which the particles interactions are through nearest and second nearest neighbors. The kernel function can be explicitly expressed in a matrix form. The analysis focuses on the decay properties, both spatially and temporally, revealing a power-law behavior in both cases. The dependence on the level of coarse-graining is also studied.

Adsorption property of granulated Prussian blue adsorbent on radioactive cesium was evaluated for efficient decontamination in Fukushima area. The adsorbent was found to show an inflective adsorption isotherm, which was expressed by extended Langmuir formula with three adsorption sites. Adsorption speeds of each site were evaluated by time-dependent batch experiment. The simulation using derived parameters and one-dimensional adsorption model successfully reproduced the experimental data of cesium decontamination by small and large columns. (author)

Various proteins are known to exhibit one-dimensional Brownian motion along charged rodlike polymers, such as microtubules (MTs), actin, and DNA. The electrostatic interaction between the proteins and the rodlike polymers appears to be crucial for one-dimensional Brownian motion, although the underlying mechanism has not been fully clarified. We examined the interactions of positively-charged nanoparticles composed of polyacrylamide gels with MTs. These hydrophilic nanoparticles bound to MTs and displayed one-dimensional Brownian motion in a charge-dependent manner, which indicates that nonspecific electrostatic interaction is sufficient for one-dimensional Brownian motion. The diffusion coefficient decreased exponentially with an increasing particle charge (with the exponent being 0.10 kBT per charge), whereas the duration of the interaction increased exponentially (exponent of 0.22 kBT per charge). These results can be explained semiquantitatively if one assumes that a particle repeats a cycle of binding to and movement along an MT until it finally dissociates from the MT. During the movement, a particle is still electrostatically constrained in the potential valley surrounding the MT. This entire process can be described by a three-state model analogous to the Michaelis-Menten scheme, in which the two parameters of the equilibrium constant between binding and movement, and the rate of dissociation from the MT, are derived as a function of the particle charge density. This study highlights the possibility that the weak binding interactions between proteins and rodlike polymers, e.g., MTs, are mediated by a similar, nonspecific charge-dependent mechanism. Copyright 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.

An one-dimensionalmodel to analyze the neutral atoms penetration into a hot plasma has been introduced in order to get the ionic temperature from the energy distribution of the charge exchange neutrals, which is obtained following a Montecarlo procedure. The model enhances the influence of the non homogeneous charge-exchange and temperature profiles over the energy distribution. It also shows haw the inner neutrals are screened by the plasma external layers and the dependence of the effective temperature on the charge-exchange cross section. Results agree with experimental data and with obtained through some others more elaborated models. (author) [es

Torrefaction reactor model is required for the development of reactor and process design for biomass torrefaction. In this study, a one-dimensional reactor model is developed based on the kinetic model describing volatiles components and solid evolution and the existing thermochemical model considering the heat and mass balance. The developed reactor model used the temperature and flow rate of the recycled gas as the practical manipulated variables instead of the torrefaction temperature. The temperature profiles of the gas and solid phase were generated, depending on the practical thermal conditions, using developed model. Moreover, the effect of each selected operating variables on the parameters of the torrefaction process and the effect of whole operating variables with particular energy yield were analyzed. Through the results of sensitivity analysis, it is shown that the residence time insignificantly influenced the energy yield when the flow rate of recycled gas is low. Moreover, higher temperature of recycled gas with low flow rate and residence time produces the attractive properties, including HHV and grindability, of torrefied biomass when the energy yield is specified. - Highlights: • A one-dimensional reactor model for biomass torrefaction is developed considering the heat and mass balance. • The developed reactor model uses the temperature and flow rate of the recycled gas as the practical manipulated variables. • The effect of operating variables on the parameters of the torrefaction process is analyzed. • The results of sensitivity analysis represent notable discussions which were not done by the previous researches

OTEQ is a mathematical simulation model used to characterize the fate and transport of waterborne solutes in streams and rivers. The model is formed by coupling a solute transport model with a chemical equilibrium submodel. The solute transport model is based on OTIS, a model that considers the physical processes of advection, dispersion, lateral inflow, and transient storage. The equilibrium submodel is based on MINTEQ, a model that considers the speciation and complexation of aqueous species, acid-base reactions, precipitation/dissolution, and sorption. Within OTEQ, reactions in the water column may result in the formation of solid phases (precipitates and sorbed species) that are subject to downstream transport and settling processes. Solid phases on the streambed may also interact with the water column through dissolution and sorption/desorption reactions. Consideration of both mobile (waterborne) and immobile (streambed) solid phases requires a unique set of governing differential equations and solution techniques that are developed herein. The partial differential equations describing physical transport and the algebraic equations describing chemical equilibria are coupled using the sequential iteration approach. The model's ability to simulate pH, precipitation/dissolution, and pH-dependent sorption provides a means of evaluating the complex interactions between instream chemistry and hydrologic transport at the field scale. This report details the development and application of OTEQ. Sections of the report describe model theory, input/output specifications, model applications, and installation instructions. OTEQ may be obtained over the Internet at http://water.usgs.gov/software/OTEQ.

We study the configurational structure of the point-island model for epitaxial growth in one dimension. In particular, we calculate the island gap and capture zone distributions. Our model is based on an approximate description of nucleation inside the gaps. Nucleation is described by the joint probability density pnXY(x,y), which represents the probability density to have nucleation at position x within a gap of size y. Our proposed functional form for pnXY(x,y) describes excellently the statistical behavior of the system. We compare our analytical model with extensive numerical simulations. Our model retains the most relevant physical properties of the system.

, the model is based on SARA representation of a petroleum mixture (saturates–aromatics–resins–asphaltenes), which may react differently with oxygen and produce other components (for example, light oils and coke). In total, the model contains 14 components, which may undergo 15 chemical reactions. The set...

In this article, a model for GAEs in a screw pinch plasma geometry is presented. The wave equations are derived from an ideal MHD model with corrections for finite frequency. Analytical and numerical solutions of these equations, applied to parameter sets approximating the TORTUS Tokamak and the Wendelstein WVII-AS advanced stellarator, are presented and discussed. (orig.)

A quantum chain model of multiple molecule motors is proposed as a mathematical physics theory for the microscopic modeling of classical force-velocity relation and tension transients in muscle fibers. The proposed model was a quantum many-particle Hamiltonian to predict the force-velocity relation for the slow release of muscle fibers, which has not yet been empirically defined and was much more complicated than the hyperbolic relationships. Using the same Hamiltonian model, a mathematical force-velocity relationship was proposed to explain the tension observed when the muscle was stimulated with an alternative electric current. The discrepancy between input electric frequency and the muscle oscillation frequency could be explained physically by the Doppler effect in this quantum chain model. Further more, quantum physics phenomena were applied to explore the tension time course of cardiac muscle and insect flight muscle. Most of the experimental tension transient curves were found to correspond to the theoretical output of quantum two- and three-level models. Mathematical modeling electric stimulus as photons exciting a quantum three-level particle reproduced most of the tension transient curves of water bug Lethocerus maximus. (special topic)

Topological expansion for the one-plaquette U(N) gauge model with fermions is investigated in the leading order for the Wilson and Manton actions. It is shown that the introduction of fermions does not change the phase structure

by implementing mobility differences dependent on the type of charged particle, particularly between ions and electrons; we performed experiments to substantiate the model. Results showed that the sub-saturated current and local field intensity were significantly

This paper considers an averaging procedure for the description of a particles arrangement in a Hubbard model with antiferromagnetic interactions. The arrangements are described by the devil's staircase. Completeness of the staircase is proved

We demonstrate a striking connection between the universal parametric correlations of the spectra of quantum chaotic systems and a class of integrable quantum hamiltonians. We begin by deriving a non-perturbative expression for the universal m-point correlation function of the spectra of random matrix ensembles in terms of a non-linear supermatrix σ-model. These results are shown to coincide with those from previous studies of weakly disordered metallic systems. We then introduce a continuous matrix model which describes the quantum mechanics of the Sutherland hamiltonian describing particles interacting through an inverse-square pairwise potential. We demonstrate that a field theoretic approach can be employed to determine exact analytical expressions for correlations of the quantum hamiltonian. The results, which are expressed in terms of a non-linear σ-model, are shown to coincide with those for analogous correlation functions of random matrix ensembles after an appropriate change of variables. We also discuss possible generalizations of the matrix model to higher dimensions. These results reveal a common mathematical structure which underlies branches of theoretical physics ranging from continuous matrix models to strongly interacting quantum hamiltonians, and universalities in the spectra of quantum chaotic systems. (orig.)

Presented is an analysis of an open quantum model of the time-dependent evolution of a flow field governed by the nonlinear Burgers equation in one spatial dimension. The quantum model is a system of qubits where there exists a minimum time interval in the time-dependent dynamics. Each temporally discrete unitary quantum-mechanical evolution is followed by state reduction of the quantum state. The mesoscopic behavior of this quantum model is described by a quantum Boltzmann equation with a naturally emergent entropy function and H theorem and the model obeys the detailed balance principle. The macroscopic-scale effective field theory for the quantum model is derived using a perturbative Chapman-Enskog expansion applied to the linearized quantum Boltzmann equation. The entropy function is consistent with the quantum-mechanical collision process and a Fermi-Dirac single-particle distribution function for the occupation probabilities of the qubit's energy eigenstates. Comparisons are presented between analytical predictions and numerical predictions and the agreement is excellent, indicating that the nonlinear Burgers equation with a tunable shear viscosity is the operative macroscopic scale effective field theory

We study the thermodynamic properties of spin chains of Haldane–Shastry type associated with the A N−1 root system in the presence of a uniform external magnetic field. To this end, we exactly compute the partition function of these models for an arbitrary finite number of spins. We then show that these chains are equivalent to a suitable inhomogeneous classical Ising model in a spatially dependent magnetic field, generalizing the results of Basu-Mallick et al. for the zero magnetic field case. Using the standard transfer matrix approach, we are able to compute in closed form the free energy per site in the thermodynamic limit. We perform a detailed analysis of the chains’ thermodynamics in a unified way, with special emphasis on the zero field and zero temperature limits. Finally, we provide a novel interpretation of the thermodynamic quantities of spin chains of Haldane–Shastry type as weighted averages of the analogous quantities over an ensemble of classical Ising models. - Highlights: ► Partition function of spin chains of Haldane–Shastry type in magnetic field. ► Equivalence to classical inhomogeneous Ising models. ► Free energy per site, other thermodynamic quantities in thermodynamic limit. ► Zero field, zero temperature limits. ► Thermodynamic equivalence with ensemble of classical Ising models.

A method for modeling piezoelectric-based ultrasonic acoustic-electric power and data transmission channels is presented. These channels employ piezoelectric disk transducers to convey signals across a series of physical layers using ultrasonic waves. This model decomposes the mechanical pathway of the signal into individual ultrasonic propagation layers which are generally independent of the layer's adjacent domains. Each layer is represented by a two-by-two traveling pressure wave transfer matrix which relates the forward and reverse pressure waves on one side of the layer to the pressure waves on the opposite face, where each face is assumed to be in contact with a domain of arbitrary reference acoustic impedance. A rigorous implementation of ultrasonic beam spreading is introduced and implemented within applicable domains. Compatible pressure-wave models for piezoelectric transducers are given, which relate the electric voltage and current interface of the transducer to the pressure waves on one mechanical interface while also allowing for passive acoustic loading of the secondary mechanical interface. It is also shown that the piezoelectric model's electrical interface is compatible with transmission line parameters (ABCD-parameters), allowing for connection of electronic components and networks. The model is shown to be capable of reproducing the behavior of realistic physical channels.

A numerical 1-D model of the chemical kinetics related hydrodynamics of room pressure N 2 plasma at 25 degrees C is reported. This generic discharge is assumed to take place between two cylindrical concentric electrodes, coated in a dielectric material, biased between 1 kV and 10 kV at 60Hz - 3kHz. The model includes the integration of particles conservation and the momentum equations as well as the local field approximation and the Poisson equations for the sake of completeness. The outcome shows that an accumulation of electrons takes place in the close vicinity of the higher voltage electrode, due to the electric field convergence to the internal electrode. Thus, this is a region of intense ionization whereas the generation of free radicals would occur away from the internal electrode. The model predicts no significant influence of the electric field on the heavier particles whose density remains practically constant.

Full Text Available We study the momentum distribution of the electrons in an extended periodic Anderson model, where the interaction, Ucf, between itinerant and localized electrons is taken into account. In the symmetric half-filled model, due to the increase of the interorbital interaction, the f electrons become more and more delocalized, while the itinerancy of conduction electrons decreases. Above a certain value of Ucf the f electrons become again localized together with the conduction electrons. In the less than half-filled case, we observe that Ucf causes strong correlations between the f electrons in the mixed valence regime.

The extended Bose-Hubbard model in a quadratic trap potential is studied using a finite-size density-matrix renormalization group method (DMRG). We compute the boson density profiles, the local compressibility and the hopping correlation functions. We observe the phase separation induced by the trap in all the quantities studied and conclude that the local density approximation is valid in the extended Bose-Hubbard model. From the plateaus obtained in the local compressibility it was possible to obtain the phase diagram of the homogeneous system which is in agreement with previous results.

As a first step in understanding what ground support equipment (GSE) is required to provide external cooling during the loading of 5,000 kg of xenon into 4 aluminum lined composite overwrapped pressure vessels (COPVs), a modeling analysis was performed using Microsoft Excel. The goals of the analysis were to predict xenon temperature and pressure throughout loading at the launch facility, estimate the time required to load one tank, and to get an early estimate of what provisions for cooling xenon might be needed while the tanks are being filled. The model uses the governing thermodynamic and heat transfer equations to achieve these goals. Results indicate that a single tank can be loaded in about 15 hours with reasonable external coolant requirements. The model developed in this study was successfully validated against flight and test data. The first data set is from the Dawn mission which also utilizes solar electric propulsion with xenon propellant, and the second is test data from the rapid loading of a hydrogen cylindrical COPV. The main benefit of this type of model is that the governing physical equations using bulk fluid solid temperatures can provide a quick and accurate estimate of the state of the propellant throughout loading which is much cheaper in terms of computational time and licensing costs than a Computation Fluid Dynamics (CFD) analysis while capturing the majority of the thermodynamics and heat transfer.

of Bardos et al. (1979) [7]. We use BV estimates on the density ρ and stability estimates on the potential Π in order to prove uniqueness. Furthermore, we analyze the evolution of characteristics for the original Hughes' model in one space dimension

Numerical models for the simulation of soil water processes can be used to evaluate the spatial and temporal variations of crop water requirements; this information can support the irrigation management in a rationale usage of water resources. This latter objective requires the knowledge of

Full Text Available In the current study, a numerical technique for solving one-dimensional fractional nonsteady heat transfer model is presented. We construct the second kind Chebyshev wavelet and then derive the operational matrix of fractional-order integration. The operational matrix of fractional-order integration is utilized to reduce the original problem to a system of linear algebraic equations, and then the numerical solutions obtained by our method are compared with those obtained by CAS wavelet method. Lastly, illustrated examples are included to demonstrate the validity and applicability of the technique.

The experimental advances in cold fermion gases motivates the investigation of a one-dimensional (1D) correlated electronic system by incorporating a four-body coupling. Using the low-energy field theory scheme and focusing on the weak-coupling regime, we extend the 1D Penson-Kolb-Hubbard (PKH) model at half filling. It is found that the additional four-body interaction may significantly modify the quantum phase diagram, favoring the presence of the superconducting phase even in the case of two-body repulsions.

are accounted for through both friction and acceleration as in a conventional formulation. However, in this analysis the acceleration term is both attributed geometrical effects through the area change and fluid dynamic effects through the expansion of the two-phase flow. The comparison of numerical...... is a onedimensional formulation in space and the equations incorporates the change in tubes and orifice diameter as formulated in (S. Madsen et.al., Dynamic Modeling of Phase Crossings in Two-Phase Flow, Communications in Computational Physics 12 (4), 1129-1147). The pressure changes in the flow...

We consider biased random walks on the infinite cluster of a conditional bond percolation model on the infinite ladder graph. Axelson-Fisk and Häggström established for this model a phase transition for the asymptotic linear speed \\overline{v} of the walk. Namely, there exists some critical value λ c>0 such that \\overline{v}>0 if λ \\in (0,λ c) and \\overline{v}=0 if λ ≥ λ c. We show that the speed \\overline{v} is continuous in λ on (0,∞) and differentiable on (0,λ c/2). Moreover, we characterize the derivative as a covariance. For the proof of the differentiability of \\overline{v} on (0,λ c/2), we require and prove a central limit theorem for the biased random walk. Additionally, we prove that the central limit theorem fails to hold for λ ≥ λ c/2.

We computed the actions for the 1D N = 5 σ-models with respect to the two inequivalent (2, 8, 6) multiplets. 4 supersymmetry generators are manifest, while the constraint originated by imposing the 5-th supersymmetry automatically induces a full N = 8 off-shell invariance. The resulting action coincides in the two cases and corresponds to a conformally flat 2D target satisfying a special geometry of rigid type. To obtain these results we developed a computational method (for Maple 11) which does not require the notion of superfields and is instead based on the nowadays available list of the inequivalent representations of the 1D N-extended supersymmetry. Its application to systematically analyze the σ-models off-shell invariant actions for the remaining N = 5, 6, 7, 8 (k, 8, 8 - k) multiplets, as well as for the N > 8 representations, only requires more cumbersome computations. (author)

Full Text Available The paper is focused on a numerical simulation of unsteady flow in a pipeline. The special attention is paid to a numerical model of an air valve, which has to include all possible regimes: critical/subcritical inflow and critical/subcritical outflow of air. Thermodynamic equation of subcritical mass flow was simplified to get more friendly shape of relevant equations, which enables easier solution of the problem.

''An integrated momentum model'' obtained by Meyer-Rose and widely applicable in calculations of dynamics of the thermal power systems is generalized for a case of flow of a vapour-liquid mixture with phase creep and pressure variation in the heated channel. Pressure distribution along the channel length is shown for a number of cases to be negligible. The obtained equations are found as well applicable in case pressure greatly though slowly varies in the system

Using exact Lanczos diagonalizations we have studied in detail the transition into the η-superconducting state in the Penson-Kolb model. We have shown that the transition occurs at any band filling and the critical value W c varies between W c≅-1.8t at half-filling to W c≅-2t for two particles on the lattice. The transition corresponds to abrupt and drastic change in the ground-state structure. After the transition the contribution of the one-electron band to the ground-state energy is almost suppressed. (orig.)

We show that the formal perturbation expansion of the invariant measure for the Anderson model in one dimension has singularities at all energies E 0 = 2 cos π(p/q); we derive a modified expansion near these energies that we show to have finite coefficients to all orders. Moreover, we show that the first q - 3 of them coincide with those of the naive expansion, while there is an anomaly in the (q - 2)th term. This also gives a weak disorder expansion for the Liapunov exponent and for the density of states. This generalizes previous results of Kappus and Wegner and of Derrida and Gardner

Some generalization are made on the spontaneous emission by a plane of excited atoms, described by two level atom-model, in the Δ1=1, Δm=1, transition and using the semiclassical radiation approximation -both discussed in the text. Initially, the radiation rate of an infinite plane of excited atoms is investigated, using Δ1=0, Δm=0, transition. It is shown that we can observe a limit solution depending on the coupling between field and matter. (author)

The use of shape memory alloys (SMAs) in engineering applications has increased the interest of the accuracy analysis of their thermomechanical description. This work presents an uncertainty analysis related to experimental tensile tests conducted with shape memory alloy wires. Experimental data...... are compared with numerical simulations obtained from a constitutive model with internal constraints employed to describe the thermomechanical behavior of SMAs. The idea is to evaluate if the numerical simulations are within the uncertainty range of the experimental data. Parametric analysis is also developed...

Critical indices, depending on continuous parameters in Bose-gas quantum models and Heisenberg 1/2 spin antiferromagnetic in two-dimensional space-time at zero temperature, have been calculated by means of finite-dimensional effects. In this case the long-wave asymptotics of the correlation functions is of a power character. Derivation of man asymptotics terms is reduced to the determination of a central charge in the appropriate Virassoro algebra representation and the anomalous dimension-operator spectrum in this representation. The finite-dimensional effects allow to find these values

Subcooled flow boiling occurs in many industrial applications and it is characterized by large heat transfer coefficients. However, this efficient heat transfer mechanism is limited by the critical heat flux, where the heat transfer coefficient decreases leading to a fast heater temperature excursion, potentially leading to heater melting and destruction. Subcooled flow boiling is especially important in water-cooled nuclear power reactors, where the presence of vapor bubbles in the core influences the reactor system behavior at operating and accident conditions. With the aim of verifying the subcooled flow boiling calculation models of the most important nuclear reactor thermal-hydraulic computer codes, such as RELAP5, COBRA-EN and COTHA-2tp, the main purpose of this work is to compare experimental data with results from these codes in the pressure range between 15 and 45 bar. For the pressure of 45 bar the results are in good agreement, while for low pressures (15 and 30 bar) the results start to become conflicting. Besides, as a sub-product of this analysis, a comparison among the models is also presented. (author)

Using the open-source software openfoam as the solver, a novel approach to calculate microsphere transport and deposition in a 1D human lung-equivalent trumpet model (TM) is presented. Specifically, for particle deposition in a nonlinear trumpetlike configuration a new radial force has been developed which, along with the regular drag force, generates particle trajectories toward the wall. The new semi-empirical force is a function of any given inlet volumetric flow rate, micron-particle diameter, and lung volume. Particle-deposition fractions (DFs) in the size range from 2 μm to 10 μm are in agreement with experimental datasets for different laminar and turbulent inhalation flow rates as well as total volumes. Typical run times on a single processor workstation to obtain actual total deposition results at comparable accuracy are 200 times less than that for an idealized whole-lung geometry (i.e., a 3D-1D model with airways up to 23rd generation in single-path only).

Context. Numerical simulations of star formation are becoming ever more sophisticated, incorporating new physical processes in increasingly realistic set-ups. These models are being compared to the latest observations through state-of-the-art synthetic renderings that trace the different chemical species present in the protostellar systems. The chemical evolution of the interstellar and protostellar matter is very topical, with more and more chemical databases and reaction solvers available online to the community. Aims: The current study was developed to provide a database of relatively simple numerical simulations of protostellar collapse as a template library for observations of cores and very young protostars, and for researchers who wish to test their chemical modelling under dynamic astrophysical conditions. It was also designed to identify statistical trends that may appear when running many models of the formation of low-mass stars by varying the initial conditions. Methods: A large set of 143 calculations of the gravitational collapse of an isolated sphere of gas with uniform temperature and a Bonnor-Ebert-like density profile was undertaken using a 1D fully implicit Lagrangian radiation hydrodynamics code. The parameter space covered initial masses from 0.2 to 8 M⊙, temperatures of 5-30 K, and radii 3000 ≤ R0 ≤ 30 000 AU. Results: A spread due to differing initial conditions and optical depths, was found in the thermal evolutionary tracks of the runs. Within less than an order of magnitude, all first and second Larson cores had masses and radii essentially independent of the initial conditions. Radial profiles of the gas density, velocity, and temperature were found to vary much more outside of the first core than inside. The time elapsed between the formation of the first and second cores was found to strongly depend on the first core mass accretion rate, and no first core in our grid of models lived for longer than 2000 years before the onset of

We investigate the forming of gamma-ray burst pulses with a simple onedimensional relativistic shock model. The mechanism is that a "central engine" drives forward the nearby plasma inside the fireball to generate a series of pressure waves. We give a relativistic geometric recurrence formula that connects the time when the pressure waves are produced and the time when the corresponding shocks occurred. This relation enables us to relate the pulse magnitude with the observation time. Our analysis shows that the evolution of the pressure waves leads to a fast rise and an exponential decay pulses. In determining the width of the pulses, the acceleration time is more important than that of the deceleration.

We use the transfer matrix formalism to examine the behavior of some anisotropic hard-core fluids, the centers of whose particles are constrained to a line. At large elongation and pressure, the compressibility factor βp/rho is higher than that for a system with complete alignment by a factor 1+ν that depends upon the molecular geometry. For molecules with a finite radius of curvature, e.g., ellipses, ν = d/2, while for objects with flat sides ν = d; here d is the number of orientational degrees of freedom. A connection is made to some recent studies of hard ellipsoid fluids. We also model the effect of an external field on physical adsorption and show the existence of a phase transition in certain limiting situations

We use the transfer matrix formalism to examine the behavior of some anisotropic hard-core fluids, the centers of whose particles are constrained to a line. At large elongation and pressure, the compressibility factor βp/p is higher than that for a system with complete aligment by a factor 1+υ that depends upon the molecular geometry. For molecules with a finite radius of curvature, e.g., ellipses, υ= d/2, while for objects with flat sides υ= d; here d is the number of orientational degrees of freedom. A connection is made to some recent studies of hard ellipsoid fluids. We also model the effect of an external field on physical adsorption and show the existence of a phase transition in certain limiting situations.

It is shown that the usual cut-off procedure (α cut-off parameter) employed in the boson representation of the fermion field opepators of the one-djmensional two-fermion model (TFM) is an incorrect one as the computator of the hermitean-conjugate field operators at the same space-point fails to fulfil a certain relationship which was pointed out long ago by Jordan. The complete form of the boson representation (including the zero-mode) of a single fermion field and the correct values of the cut-off parameter α is reviewed following Jordan and generalized to the TFM. The cut-off parameter α corresponds to a bandwidth cut-off and Jordan's boson representation is exact only in the limit α → 0. The additional zero-mode terms make the exact solution of the backscattering and umklapp scattering problem to be valid only if a supplementary condition is imposed on the coupling constants. Using the present bosonization technique all the inconsistencies of the TFM are removed. The one-particle Green's function and response functions of the Tomonaga-Luttinger model (TLM) are calculated and found to be identical with those obtained by direct diagram summation. The energy gap appearing in the spectrum of the TFM with backscattering and umklapp scattering for certain values of the coupling constants is shown to be proportional to the momentum transfer cut-off γ -1 which has to be kept finite while α goes to zero. Under such conditions the anticommunication relations and Jordan's commutator are invariant under the canonical transformation on the boson operators that diagonalizes the Hamiltonian of the TLM. The charge-density response function of the TFM with backscattering is perturbationally calculated up to the first order. The cut-off α -1 applies in the same way to terms which differ only by their spin indices in the expression of this response function. The charge-density response function is also evaluated at low frequencies for the exactly soluble TFM with

Using the exact diagonalization technique, we study the properties of the ground state of a spin-1/2 antiferromagnetic Heisenberg model for a zigzag polymer chain with side radicals connected to the even sites. We consider the nearest-neighbour exchange J and the next-nearest-neighbour exchange αJ along the main chain, and J 1 between the even site on the main chain and the radical site. For small α the ground state is ferrimagnetic. For α>α c1 , the ground state is a spiral phase, which is characterized by a peak of the static structure factor S(q) locating at an incommensurate value q max . For α>α c2 , the ground state is antiferromagnetic. With increasing J 1 , α c1 decreases while α c2 has a maximum at about J 1 = 0.5. For very small J 1 and α = 0.5, the spin configuration on the main chain is a product of nearest-neighbour singlets. In the antiferromagnetic phase, if J 1 is large enough the even site and the radical site form a singlet with exchange-decoupling from the odd site while the odd sites approximately form an antiferromagnetic chain

Climate change is expected to cause increases in extreme climatic events such as heavy rainstorms and rising tidal level. Heavy rainstorms are known to be serious causes of flooding problems in big cities. Thus, high density residential and commercial areas along the rivers are facing risks of being flooded. For that reason, inundated area determination is now being considered as one of the most important areas of research focus in flood forecasting. In such a context, this paper presents the development of a floodmap in determining flood-prone areas and its volumes. The areas and volumes of flood are computed by the inundated level using the existing digital elevation model (DEM) of a hypothetical catchment chosen for study. The study focuses on the application of Flood Early Warning System (Delft — FEWS, Deltares), which is designated to work with the SOBEK (Delft) to simulate the extent of stormwater on the ground surface. The results from FEWS consist of time-series of inundation maps in Image file format (PNG) and ASCII format, which are subsequently imported to ArcGIS for further calculations. In addition, FEWS results provide options to export the video clip of water spreading out over the catchment. Consequently, inundated area and volume will be determined by the water level on the ground. Final floodmap is displayed in colors created by ArcGIS. Various flood map results corresponding to climate change scenarios will be displayed in the main part of the paper.

Black carbon (BC) has been identified to play a critical role in aerosol-planetary boundary layer (PBL) interaction and further deterioration of near-surface air pollution in megacities, which has been referred to as the dome effect. However, the impacts of key factors that influence this effect, such as the vertical distribution and aging processes of BC, as well as the underlying land surface, have not been quantitatively explored yet. Here, based on available in situ measurements of meteorology and atmospheric aerosols together with the meteorology-chemistry online coupled model WRF-Chem, we conduct a set of parallel simulations to quantify the roles of these factors in influencing the BC dome effect and surface haze pollution. Furthermore, we discuss the main implications of the results to air pollution mitigation in China. We found that the impact of BC on the PBL is very sensitive to the altitude of aerosol layer. The upper-level BC, especially that near the capping inversion, is more essential in suppressing the PBL height and weakening the turbulent mixing. The dome effect of BC tends to be significantly intensified as BC mixed with scattering aerosols during winter haze events, resulting in a decrease in PBL height by more than 15 %. In addition, the dome effect is more substantial (up to 15 %) in rural areas than that in the urban areas with the same BC loading, indicating an unexpected regional impact of such an effect to air quality in countryside. This study indicates that China's regional air pollution would greatly benefit from BC emission reductions, especially those from elevated sources from chimneys and also domestic combustion in rural areas, through weakening the aerosol-boundary layer interactions that are triggered by BC.

We investigate the nonequilibrium roughening transition of a one-dimensional restricted solid-on-solid model by directly sampling the stationary probability density of a suitable order parameter as the surface adsorption rate varies. The shapes of the probability density histograms suggest a typical Ginzburg-Landau scenario for the phase transition of the model, and estimates of the "magnetic" exponent seem to confirm its mean-field critical behavior. We also found that the flipping times between the metastable phases of the model scale exponentially with the system size, signaling the breaking of ergodicity in the thermodynamic limit. Incidentally, we discovered that a closely related model not considered before also displays a phase transition with the same critical behavior as the original model. Our results support the usefulness of off-critical histogram techniques in the investigation of nonequilibrium phase transitions. We also briefly discuss in the appendix a good and simple pseudo-random number generator used in our simulations.

To evaluate the inability of a one-dimensional ground-water model to interact continuously with surrounding hydraulic head gradients, simulations using one-dimensional and two-dimensional ground-water flow models were compared. This approach used two types of models: flow-conserving one-and-two dimensional models, and one-dimensional and two-dimensional models designed to yield two-dimensional solutions. The hydraulic conductivities of controlling features were varied and model comparison was based on the travel times of marker particles. The solutions within each of the two model types compare reasonably well, but a three-dimensional solution is required to quantify the comparison

We consider the distribution function P(|ψ| 2 ) of the eigenfunction amplitude at the center-of-band (E = 0) anomaly in the one-dimensional tight-binding chain with weak uncorrelated on-site disorder (the one-dimensional Anderson model). The special emphasis is on the probability of the anomalously localized states (ALS) with |ψ| 2 much larger than the inverse typical localization length ℓ 0 . Using the recently found solution for the generating function Φ an (u, ϕ) we obtain the ALS probability distribution P(|ψ| 2 ) at |ψ| 2 ℓ 0 ≫ 1. As an auxiliary preliminary step, we found the asymptotic form of the generating function Φ an (u, ϕ) at u ≫ 1 which can be used to compute other statistical properties at the center-of-band anomaly. We show that at moderately large values of |ψ| 2 ℓ 0 , the probability of ALS at E = 0 is smaller than at energies away from the anomaly. However, at very large values of |ψ| 2 ℓ 0 , the tendency is inverted: it is exponentially easier to create a very strongly localized state at E = 0 than at energies away from the anomaly. We also found the leading term in the behavior of P(|ψ| 2 ) at small |ψ| 2 ≪ ℓ −1 0 and show that it is consistent with the exponential localization corresponding to the Lyapunov exponent found earlier by Kappus and Wegner. (paper)

This note examines the two-dimensional unsteady isothermal free surface flow of an incompressible fluid in a non-deformable, homogeneous, isotropic, and saturated porous medium (with zero recharge and neglecting capillary effects). Coupling a Boussinesq-type model for nonlinear water waves with Darcy's law, the two-dimensional flow problem is solved using one-dimensionalmodel equations including vertical effects and seepage face. In order to take into account the seepage face development, the system equations (given by the continuity and momentum equations) are completed by an integral relation (deduced from the Cauchy theorem). After testing the model against data sets available in the literature, some numerical simulations, concerning the unsteady flow through a rectangular dam (with an impermeable horizontal bottom), are presented and discussed.

A one-dimensional electrostatic sheet model of a coaxial geometry Virtual Cathode Oscillator (VCO) is presented. The cathode is centrally located and connects to a peripherally located plate electrode to form a resonant cavity, and is thus considered to be a novel design. Charge is modelled as concentric sheets about the cathode whose absolute position and velocity are determined as a function of time by solving the relativistic equations of motion. The model predicts the formation of a virtual cathode between the grid and plate electrodes for the case of a space-charge limited current. Setting the electron reflexing frequency (as a function of the grid potential) comparable with the cavity resonant frequency is predicted to improve the efficiency of microwave emission

Highlights: • An improved 1-D model of 3-D particle transport in ducts is studied. • The cases of isotropic and directional incidence are treated with the ADO method. • Accurate numerical results are reported for ducts of circular cross section. • A comparison with results of other authors is included. • The ADO method is found to be very efficient. - Abstract: An analytical discrete-ordinates solution is developed for the problem of particle transport in ducts, as described by a one-dimensionalmodel constructed with two basis functions. Two types of particle incidence are considered: isotropic incidence and incidence described by the Dirac delta distribution. Accurate numerical results are tabulated for the reflection probabilities of semi-infinite ducts and the reflection and transmission probabilities of finite ducts. It is concluded that the developed solution is more efficient than commonly used numerical implementations of the discrete-ordinates method.

The critical dynamics of the isotropic XY-model on the one-dimensional superlattice is considered in the framework of the position space renormalization group theory. The decimation transformation is introduced by considering the equations of motion of the operators associated to the excitations of the system, and it corresponds to an extension of the procedure introduced by Stinchcombe and dos Santos (J. Phys. A18, L597 (1985)) for the homogeneous lattice. The dispersion relation is obtained exactly and the static and dynamic scaling forms are explicitly determined. The dynamic critical exponent is also obtained and it is shown that it is identical to the one of the XY-model on the homogeneous chain.

We construct a two-state one-dimensional reaction-path model for ozone open → cyclic isomerization dynamics. The model is based on the intrinsic reaction coordinate connecting the cyclic and open isomers with the O 2 + O asymptote on the ground-state 1 A ′ potential energy surface obtained with the high-level ab initio method. Using this two-state model time-dependent wave packet optimal control simulations are carried out. Two possible pathways are identified along with their respective band-limited optimal control fields; for pathway 1 the wave packet initially associated with the open isomer is first pumped into a shallow well on the excited electronic state potential curve and then driven back to the ground electronic state to form the cyclic isomer, whereas for pathway 2 the corresponding wave packet is excited directly to the primary well of the excited state potential curve. The simulations reveal that the optimal field for pathway 1 produces a final yield of nearly 100% with substantially smaller intensity than that obtained in a previous study [Y. Kurosaki, M. Artamonov, T.-S. Ho, and H. Rabitz, J. Chem. Phys. 131, 044306 (2009)] using a single-state one-dimensionalmodel. Pathway 2, due to its strong coupling to the dissociation channel, is less effective than pathway 1. The simulations also show that nonlinear field effects due to molecular polarizability and hyperpolarizability are small for pathway 1 but could become significant for pathway 2 because much higher field intensity is involved in the latter. The results suggest that a practical control may be feasible with the aid of a few lowly excited electronic states for ozone isomerization

Low-dimensional iterated map models have been widely used to study action potential dynamics in isolated cardiac cells. Coupled iterated map models have also been widely used to investigate action potential propagation dynamics in one-dimensional (1D) coupled cardiac cells, however, these models are usually empirical and not carefully validated. In this study, we first developed two coupled iterated map models which are the standard forms of diffusively coupled maps and overcome the limitations of the previous models. We then determined the coupling strength and space constant by quantitatively comparing the 1D action potential duration profile from the coupled cardiac cell model described by differential equations with that of the coupled iterated map models. To further validate the coupled iterated map models, we compared the stability conditions of the spatially uniform state of the coupled iterated maps and those of the 1D ionic model and showed that the coupled iterated map model could well recapitulate the stability conditions, i.e. the spatially uniform state is stable unless the state is chaotic. Finally, we combined conduction into the developed coupled iterated map model to study the effects of coupling strength on wave stabilities and showed that the diffusive coupling between cardiac cells tends to suppress instabilities during reentry in a 1D ring and the onset of discordant alternans in a periodically paced 1D cable

Full Text Available Mass migrations in coconut slices during osmotic dehydration and drying are described using a diffusion model with boundary condition of the third kind. The osmotic dehydration experiment was performed at 35°Brix (water and sucrose and 40 °C. The convective drying experiments were performed at 50, 60 and 70 °C. The one-dimensional solution of the diffusion equation for an infinite slab was coupled with an optimizer to determine the effective mass diffusivities D and convective mass transfer coefficients h of the five processes studied. The analyses of the obtained results indicate that there is a good agreement between each experimental dataset and the corresponding simulation using D and h determined by optimization.

In this contribution we show results from numerical simulations carried out with a complex biogeochemical fluxes model coupled with a one-dimensional high-resol ution hydrodynamical model and implemented at three different locations of the north- ern Adriatic shelf . One location is directly affected by Po river influence, one has more open-sea characteristics and one is located in the Gulf of Trieste with an in- termediate behavior; emphasis is put on the comparison with observations and on the functioning of the northern Adriatic ecosystem in the three areas. The work has been performed in a climatological context and has to be considered as preliminary to the development of three-dimensional numerical simulations. Biogeochemical model parameterizations have been ameliorated with a detailed description of bacterial sub- strate utilization associated to the quality of the dissolved organic matter (DOM) in order to improve model skill in capturing the observed DOM dynamics in the basin. The coupled model has been calibrated and validated at the three locations by means of climatological datasets. Results show satisfactory model behavior in simulating local seasonal dynamics in the limit of the available boundary conditions and the one-dimensional implementation. Comparisons with available in situ measurements of primary and bacterial production and bacterial abundances have been performed in all locations. Model simulated rates and bacterial dynamics are in the same order of magnitude of observations and show a qualitatively correct time evolution. The importance of temperature as a factor controlling bacteria efficiency is investigated with sensitivity experiments on the model parameterizations. The different model be- havior and pelagic ecosystem structure developed by the model at the three location can be attributed to the local hydrodynamical features and interactions with external inputs of nutrients. The onset of the winter/spring bloom in the climatological

Large eddy simulation (LES) combined with the one-dimensional turbulence (ODT) model is used to simulate spatially developing turbulent reacting shear layers with high heat release and high Reynolds numbers. The LES-ODT results are compared to results from direct numerical simulations (DNS), for model development and validation purposes. The LES-ODT approach is based on LES solutions for momentum and pressure on a coarse grid and solutions for momentum and reactive scalars on a fine, one-dimensional, but three-dimensionally coupled ODT subgrid, which is embedded into the LES computational domain. Although one-dimensional, all three velocity components are transported along the ODT domain. The low-dimensional spatial and temporal resolution of the subgrid scales describe a new modeling paradigm, referred to as autonomous microstructure evolution (AME) models, which resolve the multiscale nature of turbulence down to the Kolmogorv scales. While this new concept aims to mimic the turbulent cascade and to reduce the number of input parameters, AME enables also regime-independent combustion modeling, capable to simulate multiphysics problems simultaneously. The LES as well as the one-dimensional transport equations are solved using an incompressible, low Mach number approximation, however the effects of heat release are accounted for through variable density computed by the ideal gas equation of state, based on temperature variations. The computations are carried out on a three-dimensional structured mesh, which is stretched in the transverse direction. While the LES momentum equation is integrated with a third-order Runge-Kutta time-integration, the time integration at the ODT level is accomplished with an explicit Forward-Euler method. Spatial finite-difference schemes of third (LES) and first (ODT) order are utilized and a fully consistent fractional-step method at the LES level is used. Turbulence closure at the LES level is achieved by utilizing the Smagorinsky

Using receiver functions, Rayleigh wave phase velocity dispersion determined from ambient noise and teleseismic earthquakes, and Rayleigh wave horizontal to vertical ground motion amplitude ratios from earthquakes observed across the PLUTONS seismic array, we construct a one-dimensional (1-D) S-wave velocity (Vs) seismic model with uncertainties for Uturuncu volcano, Bolivia, located in the central Andes and overlying the eastward-subducting Nazca plate. We find a fast upper crustal lid placed upon a low-velocity zone (LVZ) in the mid-crust. By incorporating all three types of measurements with complimentary sensitivity, we also explore the average density and Vp/Vs (ratio of P-wave to S-wave velocity) structures beneath the young silicic volcanic field. We observe slightly higher Vp/Vs and a decrease in density near the LVZ, which implies a dacitic source of the partially molten magma body. We exploit the impact of the 1-D model on full moment tensor inversion for the two largest local earthquakes recorded (both magnitude ∼3), demonstrating that the 1-D model influences the waveform fits and the estimated source type for the full moment tensor. Our 1-D model can serve as a robust starting point for future efforts to determine a three-dimensional velocity model for Uturuncu volcano.

Using receiver functions, Rayleigh wave phase velocity dispersion determined from ambient noise and teleseismic earthquakes, and Rayleigh wave horizontal to vertical ground motion amplitude ratios from earthquakes observed across the PLUTONS seismic array, we construct a one-dimensional (1-D) S-wave velocity (Vs) seismic model with uncertainties for Uturuncu volcano, Bolivia, located in the central Andes and overlying the eastward-subducting Nazca plate. We find a fast upper crustal lid placed upon a low-velocity zone (LVZ) in the mid-crust. By incorporating all three types of measurements with complimentary sensitivity, we also explore the average density and Vp/Vs (ratio of P-wave to S-wave velocity) structures beneath the young silicic volcanic field. We observe slightly higher Vp/Vs and a decrease in density near the LVZ, which implies a dacitic source of the partially molten magma body. We exploit the impact of the 1-D model on full moment tensor inversion for the two largest local earthquakes recorded (both magnitude ∼3), demonstrating that the 1-D model influences the waveform fits and the estimated source type for the full moment tensor. Our 1-D model can serve as a robust starting point for future efforts to determine a three-dimensional velocity model for Uturuncu volcano.

This report presents the development and verification of a one- dimensional finite element model of water flow through saturated- unsaturated media. 1DFEMWATER is very flexible and capable of modeling a wide range of real-world problems. The model is designed to (1) treat heterogeneous media consisting of many geologic formations; (2) consider distributed and point sources/sinks that are spatially and temporally variable; (3) accept prescribed initial conditions or obtain them from steady state simulations; (4) deal with transient heads distributed over the Dirichlet boundary; (5) handle time-dependent fluxes caused by pressure gradient on the Neumann boundary; (6) treat time-dependent total fluxes (i.e., the sum of gravitational fluxes and pressure-gradient fluxes) on the Cauchy boundary; (7) automatically determine variable boundary conditions of evaporation, infiltration, or seepage on the soil-air interface; (8) provide two options for treating the mass matrix (consistent and lumping); (9) provide three alternatives for approximating the time derivative term (Crank-Nicolson central difference, backward difference, and mid-difference); (10) give three options (exact relaxation, underrelaxation, and overrelaxation) for estimating the nonlinear matrix; (11) automatically reset the time step size when boundary conditions or source/sinks change abruptly; and (12) check mass balance over the entire region for every time step. The model is verified with analytical solutions and other numerical models for three examples

Air layers are regions of air between structural elements than can be found in numerous spacecraft structures. The space between folded solar panels and between antennas and a satellite's body are cases of air layers. For some cases, depending on the flexibility of the contiguous structures, the contribution of air layers can modify noticeably the dynamic response of a spacecraft structure. The analysis of these problems in detailed numerical models as Finite and Boundary Element models are characterised by a very small element size because of the requirements imposed by the thickness of the air layers and the fluid-structure interface. Then, a preliminary assessment of the influence of the air layer allows optimizing the development work flow of these elements. This work presents a methodology to preliminarily assess the influence of air layers in the structural response. The methodology is based on the definition of simplified one-dimensionalmodels for the structure and the air gaps. The study of these simple models can be a useful tool to determine the degree of influence of the air layers in the system. Along with the introduction of the methodology a study on several of the model parameters as the number of degrees of freedom for the air layer or the structure is presented. The performance of the methodology is illustrated with results for several cases including actual spacecraft structures.

Full Text Available The long Middle Route of the South to North Water Transfer Project is composed of complex hydraulic structures (aqueduct, tunnel, control gate, diversion, culvert, and diverted siphon, which generate complex flow patterns. It is vital to simulate the flow patterns through hydraulic structures, but it is a challenging work to protect water quality and maintain continuous water transfer. A one-dimensional hydrodynamic and water quality model was built to understand the flow and pollutant movement in this project. Preissmann four-point partial-node implicit scheme was used to solve the governing equations in this study. Water flow and pollutant movement were appropriately simulated and the results indicated that this water quality model was comparable to MIKE 11 and had a good performance and accuracy. Simulation accuracy and model uncertainty were analyzed. Based on the validated water quality model, six pollution scenarios (Q1 = 10 m3/s, Q2 = 30 m3/s, and Q3 = 60 m3/s for volatile phenol (VOP and contaminant mercury (Hg were simulated for the MRP. Emergent pollution accidents were forecasted and changes of water quality were analyzed according to the simulations results, which helped to guarantee continuously transferring water for a large water transfer project.

This technical report describes the new one-dimensional (1D) hydrodynamic and sediment transport model EFDC1D. This model that can be applied to stream networks. The model code and two sample data sets are included on the distribution CD. EFDC1D can simulate bi-directional unstea...

Various proteins are known to exhibit one-dimensional Brownian motion along charged rodlike polymers, such as microtubules (MTs), actin, and DNA. The electrostatic interaction between the proteins and the rodlike polymers appears to be crucial for one-dimensional Brownian motion, although the underlying mechanism has not been fully clarified. We examined the interactions of positively-charged nanoparticles composed of polyacrylamide gels with MTs. These hydrophilic nanoparticles bound to MTs ...

We revisit the equilibrium one-dimensional ϕ4 model from the dynamical systems point of view. We find an infinite number of periodic orbits which are computationally stable. At the same time some of the orbits are found to exhibit positive Lyapunov exponents! The periodic orbits confine every particle in a periodic chain to trace out either the same or a mirror-image trajectory in its two-dimensional phase space. These ;computationally stable; sets of pairs of single-particle orbits are either symmetric or antisymmetric to the very last computational bit. In such a periodic chain the odd-numbered and even-numbered particles' coordinates and momenta are either identical or differ only in sign. ;Positive Lyapunov exponents; can and do result if an infinitesimal perturbation breaking a perfect two-dimensional antisymmetry is introduced so that the motion expands into a four-dimensional phase space. In that extended space a positive exponent results. We formulate a standard initial condition for the investigation of the microcanonical chaotic number dependence of the model. We speculate on the uniqueness of the model's chaotic sea and on the connection of such collections of deterministic and time-reversible states to the Second Law of Thermodynamics.

A novel one-dimensional (1D) analytical model is proposed for quantifying the breakdown voltage of a reduced surface field (RESURF) lateral power device fabricated on silicon on an insulator (SOI) substrate. We assume that the charges in the depletion region contribute to the lateral PN junctions along the diagonal of the area shared by the lateral and vertical depletion regions. Based on the assumption, the lateral PN junction behaves as a linearly graded junction, thus resulting in a reduced surface electric field and high breakdown voltage. Using the proposed model, the breakdown voltage as a function of device parameters is investigated and compared with the numerical simulation by the TCAD tools. The analytical results are shown to be in fair agreement with the numerical results. Finally, a new RESURF criterion is derived which offers a useful scheme to optimize the structure parameters. This simple 1D model provides a clear physical insight into the RESURF effect and a new explanation on the improvement in breakdown voltage in an SOI RESURF device. (paper)

Large-scale free-piston driven expansion tubes have uniquely high total pressure capabilities which make them an important resource for development of access-to-space scramjet engine technology. However, many aspects of their operation are complex, and their test flows are fundamentally unsteady and difficult to measure. While computational fluid dynamics methods provide an important tool for quantifying these flows, these calculations become very expensive with increasing facility size and therefore have to be carefully constructed to ensure sufficient accuracy is achieved within feasible computational times. This study examines modelling strategies for a Mach 10 scramjet test condition developed for The University of Queensland's X3 facility. The present paper outlines the challenges associated with test flow reconstruction, describes the experimental set-up for the X3 experiments, and then details the development of an experimentally tuned quasi-one-dimensional CFD model of the full facility. The 1-D model, which accurately captures longitudinal wave processes, is used to calculate the transient flow history in the shock tube. This becomes the inflow to a higher-fidelity 2-D axisymmetric simulation of the downstream facility, detailed in the Part 2 companion paper, leading to a validated, fully defined nozzle exit test flow.

Large-scale free-piston driven expansion tubes have uniquely high total pressure capabilities which make them an important resource for development of access-to-space scramjet engine technology. However, many aspects of their operation are complex, and their test flows are fundamentally unsteady and difficult to measure. While computational fluid dynamics methods provide an important tool for quantifying these flows, these calculations become very expensive with increasing facility size and therefore have to be carefully constructed to ensure sufficient accuracy is achieved within feasible computational times. This study examines modelling strategies for a Mach 10 scramjet test condition developed for The University of Queensland's X3 facility. The present paper outlines the challenges associated with test flow reconstruction, describes the experimental set-up for the X3 experiments, and then details the development of an experimentally tuned quasi-one-dimensional CFD model of the full facility. The 1-D model, which accurately captures longitudinal wave processes, is used to calculate the transient flow history in the shock tube. This becomes the inflow to a higher-fidelity 2-D axisymmetric simulation of the downstream facility, detailed in the Part 2 companion paper, leading to a validated, fully defined nozzle exit test flow.

Cilia are elastic hairlike protuberances of the cell membrane found in various unicellular organisms and in several tissues of most living organisms. In some tissues such as the airway tissues of the lung, the coordinated beating of cilia induces a fluid flow of crucial importance as it allows the continuous cleaning of our bronchia, known as mucociliary clearance. While most of the models addressing the question of collective dynamics and metachronal wave consider homogeneous carpets of cilia, experimental observations rather show that cilia clusters are heterogeneously distributed over the tissue surface. The purpose of this paper is to investigate the role of spatial heterogeneity on the coherent beating of cilia using a very simple one-dimensionalmodel for cilia known as the rower model. We systematically study systems consisting of a few rowers to hundreds of rowers and we investigate the conditions for the emergence of collective beating. When considering a small number of rowers, a phase drift occurs, hence, a bifurcation in beating frequency is observed as the distance between rower clusters is changed. In the case of many rowers, a distribution of frequencies is observed. We found in particular the pattern of the patchy structure that shows the best robustness in collective beating behavior, as the density of cilia is varied over a wide range.

In this paper we investigate flow localization in viscoplastic slender bars subjected to dynamic tension. We explore loading rates above the critical impact velocity: the wave initiated in the impacted end by the applied velocity is the trigger for the localization of plastic deformation. The problem has been addressed using two kinds of numerical simulations: (1) one-dimensional finite difference calculations and (2) axisymmetric finite element computations. The latter calculations have been used to validate the capacity of the finite difference model to describe plastic flow localization at high impact velocities. The finite difference model, which highlights due to its simplicity, allows to obtain insights into the role played by the strain rate and temperature sensitivities of the material in the process of dynamic flow localization. Specifically, we have shown that viscosity can stabilize the material behavior to the point of preventing the appearance of the critical impact velocity. This is a key outcome of our investigation, which, to the best of the authors' knowledge, has not been previously reported in the literature.

We study analytically and numerically a one-dimensionalmodel of interacting line defects (steps) fluctuating on a vicinal crystal. Our goal is to formulate and validate analytical techniques for approximately solving systems of coupled nonlinear stochastic differential equations (SDEs) governing fluctuations in surface motion. In our analytical approach, the starting point is the Burton-Cabrera-Frank (BCF) model by which step motion is driven by diffusion of adsorbed atoms on terraces and atom attachment-detachment at steps. The step energy accounts for entropic and nearest-neighbor elastic-dipole interactions. By including Gaussian white noise to the equations of motion for terrace widths, we formulate large systems of SDEs under different choices of diffusion coefficients for the noise. We simplify this description via (i) perturbation theory and linearization of the step interactions and, alternatively, (ii) a mean-field (MF) approximation whereby widths of adjacent terraces are replaced by a self-consistent field but nonlinearities in step interactions are retained. We derive simplified formulas for the time-dependent terrace-width distribution (TWD) and its steady-state limit. Our MF analytical predictions for the TWD compare favorably with kinetic Monte Carlo simulations under the addition of a suitably conservative white noise in the BCF equations.

Understanding the mechanism of plasma build-up in vacuum arcs is essential in many fields of physics. A one-dimensional particle-in-cell computer simulation model is presented, which models the plasma developing from a field emitter tip under electrical breakdown conditions, taking into account the relevant physical phenomena. As a starting point, only an external electric field and an initial enhancement factor of the tip are assumed. General requirements for plasma formation have been identified and formulated in terms of the initial local field and a critical neutral density. The dependence of plasma build-up on tip melting current, the evaporation rate of neutrals and external circuit time constant has been investigated for copper and simulations imply that arcing involves melting currents around 0.5-1 A/mu m(2),evaporation of neutrals to electron field emission ratios in the regime 0.01 - 0.05, plasma build-up timescales in the order of similar to 1 - 10 ns and two different regimes depending on initial ...

Highlights: • Seven heat transfer mechanisms are studied numerically by the model. • A semi-empirical method is proposed to account for the transition boiling effect. • The parametric effects on the heat transfer mechanisms are investigated. • The thermal non-equilibrium phenomenon between vapor and droplets is investigated. - Abstract: The objective of this paper is to develop a one-dimensional semi-empirical model for the dispersed flow film boiling considering transition boiling effects. The proposed model consists of conservation equations, i.e., vapor mass, vapor energy, droplet mass and droplet momentum conservation, and a set of closure relations to address the interactions among wall, vapor and droplets. The results show that the transition boiling effect is of vital importance in the dispersed flow film boiling regime, since the flowing situation in the downstream would be influenced by the conditions in the upstream. In addition, the present paper, through evaluating the vapor temperature and the amount of heat transferred to droplets, investigates the thermal non-equilibrium phenomenon under different flowing conditions. Comparison of the wall temperature predictions with the 1394 experimental data in the literature, the present model ranging from system pressure of 30–140 bar, heat flux of 204–1837 kW/m{sup 2} and mass flux of 380–5180 kg/m{sup 2} s, shows very good agreement with RMS of 8.80% and standard deviation of 8.81%. Moreover, the model well depicts the thermal non-equilibrium phenomenon for the dispersed flow film boiling.

The impact of confinement on self-assembly of particles interacting with short-range attraction and long-range repulsion potential is studied for thermodynamic states corresponding to local ordering of clusters or layers in the bulk. Exact and asymptotic expressions for the local density and for the effective potential between the confining surfaces are obtained for a one-dimensional lattice model introduced by J. Pȩkalski et al. [J. Chem. Phys. 138, 144903 (2013)]. The simple asymptotic formulas are shown to be in good quantitative agreement with exact results for slits containing at least 5 layers. We observe that the incommensurability of the system size and the average distance between the clusters or layers in the bulk leads to structural deformations that are different for different values of the chemical potential μ. The change of the type of defects is reflected in the dependence of density on μ that has a shape characteristic for phase transitions. Our results may help to avoid misinterpretation of the change of the type of defects as a phase transition in simulations of inhomogeneous systems. Finally, we show that a system confined by soft elastic walls may exhibit bistability such that two system sizes that differ approximately by the average distance between the clusters or layers are almost equally probable. This may happen when the equilibrium separation between the soft boundaries of an empty slit corresponds to the largest stress in the confined self-assembling system.

A tuning method for FLake, a one-dimensional (1-D) freshwater lake model, is applied for the individual tuning of 244 globally distributed large lakes using observed lake surface water temperatures (LSWTs) derived from along-track scanning radiometers (ATSRs). The model, which was tuned using only three lake properties (lake depth, snow and ice albedo and light extinction coefficient), substantially improves the measured mean differences in various features of the LSWT annual cycle, including the LSWTs of saline and high altitude lakes, when compared to the observed LSWTs. Lakes whose lake-mean LSWT persists below 1 °C for part of the annual cycle are considered to be seasonally ice-covered. For trial seasonally ice-covered lakes (21 lakes), the daily mean and standard deviation (2σ) of absolute differences between the modelled and observed LSWTs are reduced from 3.07 °C ± 2.25 °C to 0.84 °C ± 0.51 °C by tuning the model. For all other trial lakes (14 non-ice-covered lakes), the improvement is from 3.55 °C ± 3.20 °C to 0.96 °C ± 0.63 °C. The post tuning results for the 35 trial lakes (21 seasonally ice-covered lakes and 14 non-ice-covered lakes) are highly representative of the post-tuning results of the 244 lakes. For the 21 seasonally ice-covered lakes, the modelled response of the summer LSWTs to changes in snow and ice albedo is found to be statistically related to lake depth and latitude, which together explain 0.50 (R2adj, p = 0.001) of the inter-lake variance in summer LSWTs. Lake depth alone explains 0.35 (p = 0.003) of the variance. Lake characteristic information (snow and ice albedo and light extinction coefficient) is not available for many lakes. The approach taken to tune the model, bypasses the need to acquire detailed lake characteristic values. Furthermore, the tuned values for lake depth, snow and ice albedo and light extinction coefficient for the 244 lakes provide some guidance on improving FLake LSWT modelling.

Concentrated solar power plants have attracted increasing interest from researchers and governments all over the world in recent years. An important part of these plants is the storage system which improves dispatchability and makes the plant more reliable. In this paper, a one-dimensional transi...

This user's guide describes a mathematical model for predicting the transport of mixed sizes of sediment by flow in networks of one-dimensional open channels. The simulation package is useful for general sediment routing problems, prediction of erosion and deposition following dam removal, and scour in channels at road embankment crossings or other artificial structures. The model treats input hydrographs as stepwise steady-state, and the flow computation algorithm automatically switches between sub- and supercritical flow as dictated by channel geometry and discharge. A variety of boundary conditions including weirs and rating curves may be applied both external and internal to the flow network. The model may be used to compute flow around islands and through multiple openings in embankments, but the network must be 'simple' in the sense that the flow directions in all channels can be specified before simulation commences. The location and shape of channel banks are user specified, and all bedelevation changes take place between these banks and above a user-specified bedrock elevation. Computation of sediment-transport emphasizes the sand-size range (0.0625-2.0 millimeter) but the user may select any desired range of particle diameters including silt and finer (user may set the original bed-sediment composition of any number of layers of known thickness. The model computes the time evolution of total transport and the size composition of bed- and suspended-load sand through any cross section of interest. It also tracks bed -surface elevation and size composition. The model is written in the FORTRAN programming language for implementation on personal computers using the WINDOWS operating system and, along with certain graphical output display capability, is accessed from a graphical user interface (GUI). The GUI provides a framework for selecting input files and parameters of a number of components of the sediment-transport process. There are no restrictions in the

The classical Garman-Kohlhagen model for the currency exchange assumes that the domestic and foreign currency risk-free interest rates are constant and the exchange rate follows a log-normal diffusion process.In this paper we consider the general case, when exchange rate evolves according to arbitrary one-dimensional diffusion process with local volatility that is the function of time and the current exchange rate and where the domestic and foreign currency risk-free interest rates may be arbitrary continuous functions of time. First non-trivial problem we encounter in time-dependent case is the continuity in time argument of the value function of the American put option and the regularity properties of the optimal exercise boundary. We establish these properties based on systematic use of the monotonicity in volatility for the value functions of the American as well as European options with convex payoffs together with the Dynamic Programming Principle and we obtain certain type of comparison result for the value functions and corresponding exercise boundaries for the American puts with different strikes, maturities and volatilities.Starting from the latter fact that the optimal exercise boundary curve is left continuous with right-hand limits we give a mathematically rigorous and transparent derivation of the significant early exercise premium representation for the value function of the American foreign exchange put option as the sum of the European put option value function and the early exercise premium.The proof essentially relies on the particular property of the stochastic integral with respect to arbitrary continuous semimartingale over the predictable subsets of its zeros. We derive from the latter the nonlinear integral equation for the optimal exercise boundary which can be studied by numerical methods

We discuss the exit probability of the one-dimensional q-voter model and present tools to obtain estimates about this probability, both through simulations in large networks (around 107 sites) and analytically in the limit where the network is infinitely large. We argue that the result E(ρ )=ρq/ρq+(1-ρ)q, that was found in three previous works [F. Slanina, K. Sznajd-Weron, and P. Przybyła, Europhys. Lett. 82, 18006 (2008), 10.1209/0295-5075/82/18006; R. Lambiotte and S. Redner, Europhys. Lett. 82, 18007 (2008), 10.1209/0295-5075/82/18007, for the case q =2; and P. Przybyła, K. Sznajd-Weron, and M. Tabiszewski, Phys. Rev. E 84, 031117 (2011), 10.1103/PhysRevE.84.031117, for q >2] using small networks (around 103 sites), is a good approximation, but there are noticeable deviations that appear even for small systems and that do not disappear when the system size is increased (with the notable exception of the case q =2). We also show that, under some simple and intuitive hypotheses, the exit probability must obey the inequality ρq/ρq+(1-ρ)≤E(ρ)≤ρ/ρ +(1-ρ)q in the infinite size limit. We believe this settles in the negative the suggestion made [S. Galam and A. C. R. Martins, Europhys. Lett. 95, 48005 (2001), 10.1209/0295-5075/95/48005] that this result would be a finite size effect, with the exit probability actually being a step function. We also show how the result that the exit probability cannot be a step function can be reconciled with the Galam unified frame, which was also a source of controversy.

A modification of the one-dimensional Fermi accelerator model is considered in this work. The dynamics of a classical particle of mass m, confined to bounce elastically between two rigid walls where one is described by a nonlinear van der Pol type oscillator while the other one is fixed, working as a reinjection mechanism of the particle for a next collision, is carefully made by the use of a two-dimensional nonlinear mapping. Two cases are considered: (i) the situation where the particle has mass negligible as compared to the mass of the moving wall and does not affect the motion of it; and (ii) the case where collisions of the particle do affect the movement of the moving wall. For case (i) the phase space is of mixed type leading us to observe a scaling of the average velocity as a function of the parameter (χ) controlling the nonlinearity of the moving wall. For large χ, a diffusion on the velocity is observed leading to the conclusion that Fermi acceleration is taking place. On the other hand, for case (ii), the motion of the moving wall is affected by collisions with the particle. However, due to the properties of the van der Pol oscillator, the moving wall relaxes again to a limit cycle. Such kind of motion absorbs part of the energy of the particle leading to a suppression of the unlimited energy gain as observed in case (i). The phase space shows a set of attractors of different periods whose basin of attraction has a complicated organization.

The paper presents a model of fixed and variable geometry turbines. The aim of this model is to provide an efficient boundary condition to model turbocharged internal combustion engines with zero- and one-dimensional gas dynamic codes. The model is based from its very conception on the measured characteristics of the turbine. Nevertheless, it is capable of extrapolating operating conditions that differ from those included in the turbine maps, since the engines usually work within these zones. The presented model has been implemented in a one-dimensional gas dynamic code and has been used to calculate unsteady operating conditions for several turbines. The results obtained have been compared with success against pressure-time histories measured upstream and downstream of the turbine during on-engine operation. (author)

The paper presents a model of fixed and variable geometry turbines. The aim of this model is to provide an efficient boundary condition to model turbocharged internal combustion engines with zero- and one-dimensional gas dynamic codes. The model is based from its very conception on the measured characteristics of the turbine. Nevertheless, it is capable of extrapolating operating conditions that differ from those included in the turbine maps, since the engines usually work within these zones. The presented model has been implemented in a one-dimensional gas dynamic code and has been used to calculate unsteady operating conditions for several turbines. The results obtained have been compared with success against pressure-time histories measured upstream and downstream of the turbine during on-engine operation

A one-dimensional water quality model, Gulf of Mexico Dissolved Oxygen Model (GoMDOM-1D), was developed to simulate phytoplankton, carbon, nutrients, and dissolved oxygen in Gulf of Mexico. The model was calibrated and corroborated against a comprehensive set of field observation...

Highlights: → CCFL in the hot leg of a PWR with ECC Injection. → Three-Field Model of counter flowing water film and entrained droplets. → Flow of steam can cause a hydraulic jump in the supercritical flow of water. → Condensation of steam on subcooled water increases the required flow for hydraulic jump. → Better agreement with UPTF experimental data than Wallis-type correlation. - Abstract: A one-dimensional three-field model was developed to predict the flow of liquid and vapor that results from countercurrent flow of water injected into the hot leg of a PWR and the oncoming steam flowing from the upper plenum. The model solves the conservation equations for mass, momentum, and energy in a continuous-vapor field, a continuous-liquid field, and a dispersed-liquid (entrained-droplet) field. Single-effect experiments performed in the upper plenum test facility (UPTF) of the former SIEMENS KWU (now AREVA) at Mannheim, Germany, were used to validate the countercurrent flow limitation (CCFL) model in case of emergency core cooling water injection into the hot legs. Subcooled water and saturated steam flowed countercurrent in a horizontal pipe with an inside diameter of 0.75 m. The flow of injected water was varied from 150 kg/s to 400 kg/s, and the flow of steam varied from 13 kg/s to 178 kg/s. The subcooling of the liquid ranged from 0 K to 104 K. The velocity of the water at the injection point was supercritical (greater than the celerity of a gravity wave) for all the experiments. The three-field model was successfully used to predict the experimental data, and the results from the model provide insight into the mechanisms that influence the flows of liquid and vapor during countercurrent flow in a hot leg. When the injected water was saturated and the flow of steam was small, all or most of the injected water flowed to the upper plenum. Because the velocity of the liquid remained supercritical, entrainment of droplets was suppressed. When the injected

A solution to the extended Hubbard Hamiltonian for the case of two-particles in an infinite one-dimensional lattice is presented, using a real-space mapping method and the Green function technique. This Hamiltonian considers the on-site (U) and the nearest-neighbor (V) interactions. The method is based on mapping the correlated many-body problem onto an equivalent site-impurity tight-binding one in a higher dimensional space. In this new space we obtained the analytical solution for the ground state binding energy. Results are in agreement with the numerical solution obtained previously [1], and with those obtained in the reciprocal space [2]. (Author)

This article is a theoretical discussion that links Marcuse's concept of one-dimensional society and the Great Refusal with critical race theory in order to achieve a more robust interrogation of whiteness. The author argues that in the context of the United States, the one-dimensionality that Marcuse condemns in "One-Dimensional Man" is best…

In WWTP models, the accurate assessment of solids inventory in bioreactors equipped with solidliquid separators, mostly described using one-dimensional (1-D) secondary settling tank (SST) models, is the most fundamental requirement of any calibration procedure. Scientific knowledge...... of the solids settling behaviour is investigated. It is found that the settler behaviour, simulated by the hyperbolic model, can introduce significant errors into the approximation of the solids retention time and thus solids inventory of the system. We demonstrate that these impacts can potentially cause...

We study the time evolution of quantum one-dimensional gapless systems evolving from initial states with a domain wall. We generalize the path integral imaginary time approach that together with boundary conformal field theory allows us to derive the time and space dependence of general correlation functions. The latter are explicitly obtained for the Ising universality class, and the typical behavior of one- and two-point functions is derived for the general case. Possible connections with the stochastic Loewner evolution are discussed and explicit results for one-point time dependent averages are obtained for generic κ for boundary conditions corresponding to stochastic Loewner evolution. We use this set of results to predict the time evolution of the entanglement entropy and obtain the universal constant shift due to the presence of a domain wall in the initial state

Full Text Available Purpose: The main objective of this study is to develop a model for solving the onedimensional cutting stock problem in the wood working industry, and develop a program for its implementation. Design/methodology/approach: This study adopts the pattern oriented approach in the formulation of the cutting stock model. A pattern generation algorithm was developed and coded using Visual basic.NET language. The cutting stock model developed is a Linear Programming (LP Model constrained by numerous feasible patterns. A LP solver was integrated with the pattern generation algorithm program to develop a one - dimensional cutting stock model application named GB Cutting Stock Program. Findings and Originality/value: Applying the model to a real life optimization problem significantly reduces material waste (off-cuts and minimizes the total stock used. The result yielded about 30.7% cost savings for company-I when the total stock materials used is compared with the former cutting plan. Also, to evaluate the efficiency of the application, Case I problem was solved using two top commercial 1D-cutting stock software. The results show that the GB program performs better when related results were compared. Research limitations/implications: This study round up the linear programming solution for the number of pattern to cut. Practical implications: From Managerial perspective, implementing optimized cutting plans increases productivity by eliminating calculating errors and drastically reducing operator mistakes. Also, financial benefits that can annually amount to millions in cost savings can be achieved through significant material waste reduction. Originality/value: This paper developed a linear programming onedimensional cutting stock model based on a pattern generation algorithm to minimize waste in the wood working industry. To implement the model, the algorithm was coded using VisualBasic.net and linear programming solver called lpsolvedll (dynamic

We calculate the electronic structure and exchange interactions in a copper(II)phthalocyanine [Cu(II)Pc] crystal as a one-dimensional molecular chain using hybrid exchange density functional theory (DFT). In addition, the intermolecular exchange interactions are also calculated in a molecular dimer using Green’s function perturbation theory (GFPT) to illustrate the underlying physics. We find that the exchange interactions depend strongly on the stacking angle, but weakly on the sliding angle (defined in the text). The hybrid DFT calculations also provide an insight into the electronic structure of the Cu(II)Pc molecular chain and demonstrate that on-site electron correlations have a significant effect on the nature of the ground state, the band gap, and magnetic excitations. The exchange interactions predicted by our DFT calculations and GFPT calculations agree qualitatively with the recent experimental results on newly found η-Cu(II)Pc and the previous results for the α and β phases. This work provides a reliable theoretical basis for the further application of Cu(II)Pc to molecular spintronics and organic-based quantum information processing.

We present a computationally efficient framework to compute the neutral flux in high aspect ratio structures during three-dimensional plasma etching simulations. The framework is based on a one-dimensional radiosity approach and is applicable to simulations of convex rotationally symmetric holes and convex symmetric trenches with a constant cross-section. The framework is intended to replace the full three-dimensional simulation step required to calculate the neutral flux during plasma etching simulations. Especially for high aspect ratio structures, the computational effort, required to perform the full three-dimensional simulation of the neutral flux at the desired spatial resolution, conflicts with practical simulation time constraints. Our results are in agreement with those obtained by three-dimensional Monte Carlo based ray tracing simulations for various aspect ratios and convex geometries. With this framework we present a comprehensive analysis of the influence of the geometrical properties of high aspect ratio structures as well as of the particle sticking probability on the neutral particle flux.

A new physical cause for a temperature-dependent double peak in exciton systems is put forward within a thermal equilibrium approach for the calculation of optical properties of exciton systems. Indeed, it is found that one-dimensional exciton systems with only one molecule per unit cell can have an absorption spectrum characterized by a double peak provided that the coupling between excitations in different molecules is positive. The two peaks, whose relative intensities vary with temperature, are located around the exciton band edges, being separated by an energy of approximately 4V, where V is the average coupling between nearest neighbours. For small amounts of diagonal and off-diagonal disorder, the contributions from the intermediate states in the band are also visible as intermediate structure between the two peaks, this being enhanced for systems with periodic boundary conditions. At a qualitative level, these results correlate well with experimental observations in the molecular aggregates of the thiacarbocyanine dye THIATS and in the organic crystals of acetanilide and N-methylacetamide.

A new physical cause for a temperature-dependent double peak in exciton systems is put forward within a thermal equilibrium approach for the calculation of optical properties of exciton systems. Indeed, it is found that one-dimensional exciton systems with only one molecule per unit cell can have an absorption spectrum characterized by a double peak provided that the coupling between excitations in different molecules is positive. The two peaks, whose relative intensities vary with temperature, are located around the exciton band edges, being separated by an energy of approximately 4V, where V is the average coupling between nearest neighbours. For small amounts of diagonal and off-diagonal disorder, the contributions from the intermediate states in the band are also visible as intermediate structure between the two peaks, this being enhanced for systems with periodic boundary conditions. At a qualitative level, these results correlate well with experimental observations in the molecular aggregates of the thiacarbocyanine dye THIATS and in the organic crystals of acetanilide and N-methylacetamide

The main functions of spent fuel pools are to remove the residual heat from spent fuel assemblies and to perform the function of biological shielding. In the case of loss of heat removal from spent fuel pool, the fuel rods and pool water temperatures would increase continuously. After the saturated temperature is reached, due to evaporation of water the pool water level would drop, eventually causing the uncover of spent fuel assemblies, fuel overheating and fuel rods failure. This paper presents an analysis of loss of heat removal accident in spent fuel pool of BWR 4 and a comparison of two different modelling approaches. The one-dimensional system thermal-hydraulic computer code RELAP5 and CFD tool ANSYS Fluent were used for the analysis. The results are similar, but the local effects cannot be simulated using a one-dimensional code. The ANSYS Fluent calculation demonstrated that this three-dimensional treatment allows to avoid the need for many one-dimensionalmodelling assumptions in the pool modelling and enables to reduce the uncertainties associated with natural circulation flow calculation.

Fluorine ions in NaBr have associated large dipole moments with low-lying energy levels. It is well known that the dipoles were found to have equilibrium orientations in the (110) direction. A one-dimensional, double-well harmonic oscillator potential model is assumed for the relaxation rate calculation of this off-centre system. It is possible by superimposing an asymmetric potential which localizes the particle in one potential well and assuming that, the coupling between the particle and the lattice vibrations can lead to the relaxation of the system. These preliminaries theoretical studies are used to determine the height of the potential barrier between the two minima of the off-centre potential in the one-dimensional case approximation. (author). 13 refs

Full Text Available In the context of data modeling and comparisons between different fit models, Bayesian analysis calls that model best which has the largest evidence, the prior-weighted integral over model parameters of the likelihood function. Evidence calculations automatically take into account both the usual chi-squared measure and an Occam factor which quantifies the price for adding extra parameters. Applying Bayesian analysis to projections onto azimuth of 2D angular correlations from 200 GeV AuAu collisions, we consider typical model choices including Fourier series and a Gaussian plus combinations of individual cosine components. We find that models including a Gaussian component are consistently preferred over pure Fourier-series parametrizations, sometimes strongly so. For 0–5% central collisions the Gaussian-plus-dipole model performs better than Fourier Series models or any other combination of Gaussian-plus-multipoles.

The subject of this thesis is a onedimensionalmodel for a quantum system of fermions with attractive or repulsive interaction. The eigenvalues and eigenfunctions of the Hamiltonian with periodic boundary conditions are exactly determined. The knowledge of the spectrum is essentially applied on the study of the attractive gas, characterized by the presence of 'pairs' or two particles bound states. This system can be described as a gas of 'onedimensional deuterons', which has some analogy with a boson gas. Some extensive properties of the ground state have been discussed for example energy as a function of the density and magnetization, for all the values of the coupling constant. The analytic properties of the energy function are studied, but not completely resolved. Finally the elementary excitations of the phonon type are considered and the dispersion curves are given. (author) [French] On etudie un modele a une dimension pour un systeme quantique de fermions en interaction attractive ou repulsive dans un volume donne. L'ensemble des niveaux d'energie et des etats propres du systeme est determine exactement. La connaissance du spectre est surtout appliquee a l'etude du gaz attractif, interessant par la presence de 'paires' ou etats lies a deux particules. On peut decrire ce systeme comme un gaz de 'deuterons a une dimension' qui possede quelque ressemblance avec un systeme de bosons. Quelques proprietes extensives de l'etat fondamental sont donnees, comme l'energie en fonction de la densite et de la magnetisation totale, pour toute valeur de la constante de couplage. Les proprietes analytiques de la fonction energie sont etudiees sans etre completement elucidees. On aborde enfin les excitations elementaires du systeme et on etablit la courbe de dispersion d'une excitation de type phonon. (auteur)

We present an algorithm for automated simplification of 1D pipe network models. The impact of the simplifications on the flooding simulated by coupled 1D-2D models is evaluated in an Australian case study. Significant reductions of the simulation time of the coupled model are achieved by reducing...... the 1D network model. The simplifications lead to an underestimation of flooded area because interaction points between network and surface are removed and because water is transported downstream faster. These effects can be mitigated by maintaining nodes in flood-prone areas in the simplification...... and by adjusting pipe roughness to increase transport times....

The possibility of interpreting multifragmentation data obtained from heavy-ion collisions at intermediate energies, by a new type of model: the DLNA (Dynamical Limited Nuclear Aggregation) is discussed. This model is connected to a more general class of models presenting Self-Organization Criticality (SOC). It is shown that the fragment size distributions exhibit a power-law dependence comparable to those obtained in second-order phase transition or percolation models. Fluctuations in term of scaled-factorial moments and cumulants are also studied: no signal of intermittency is seen. (K.A.)

Full Text Available This article describes the development of a suite of programs that is capable of simulating the radiation properties of a random rough surface (RRS. The fundamental approach involves the generation, by fast Fourier transform (FFT built with rigorous finite difference time domain (FDTD, as the theoretical basis for the simulation of a bidirectional reflectance distribution function (BRDF of the RRS. The results are compared with the measurements and modeling of existing work to verify the feasibility of customized programming. It was found that the results of this study were a better match to the measurement data than those achieved in other modeling work.

Side channel construction is a common intervention applied to increase a river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: (1) a difference in channel slope between the side channel

Side channel construction is a common intervention applied to increase a river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: (1) a difference in channel slope between the side channel

Side channel construction is a common intervention applied to increase the river's conveyance capacity and to increase its ecological value. Past modelling efforts suggest two mechanisms affecting the morphodynamic change of a side channel: 1) a difference in channel slope between the side channel

The Kolmogorov-Johnson-Mehl-Avrami (KJMA) growth model is considered on a one-dimensional (1D) lattice. Cells can grow with constant speed and continuously nucleate on the empty sites. We offer an alternative mean-field-like approach for describing theoretically the dynamics and derive an analytical cell-size distribution function. Our method reproduces the same scaling laws as the KJMA theory and has the advantage that it leads to a simple closed form for the cell-size distribution function. It is shown that a Weibull distribution is appropriate for describing the final cell-size distribution. The results are discussed in comparison with Monte Carlo simulation data.

Highlights: • Quasiperiodic lattice models with next-nearest-neighbor hopping are studied. • Shannon information entropies are used to reflect state localization properties. • Phase diagrams are obtained for the inverse bronze and golden means, respectively. • Our studies present a more complete picture than existing works. - Abstract: We explore the reduced relative Shannon information entropies SR for a quasiperiodic lattice model with nearest- and next-nearest-neighbor hopping, where an irrational number is in the mathematical expression of incommensurate on-site potentials. Based on SR, we respectively unveil the phase diagrams for two irrationalities, i.e., the inverse bronze mean and the inverse golden mean. The corresponding phase diagrams include regions of purely localized phase, purely delocalized phase, pure critical phase, and regions with mobility edges. The boundaries of different regions depend on the values of irrational number. These studies present a more complete picture than existing works.

Highlights: • Quasiperiodic lattice models with next-nearest-neighbor hopping are studied. • Shannon information entropies are used to reflect state localization properties. • Phase diagrams are obtained for the inverse bronze and golden means, respectively. • Our studies present a more complete picture than existing works. - Abstract: We explore the reduced relative Shannon information entropies SR for a quasiperiodic lattice model with nearest- and next-nearest-neighbor hopping, where an irrational number is in the mathematical expression of incommensurate on-site potentials. Based on SR, we respectively unveil the phase diagrams for two irrationalities, i.e., the inverse bronze mean and the inverse golden mean. The corresponding phase diagrams include regions of purely localized phase, purely delocalized phase, pure critical phase, and regions with mobility edges. The boundaries of different regions depend on the values of irrational number. These studies present a more complete picture than existing works.

In a Belt-Pinch device a hot, high-β-plasma is produced by the method of fast shock-compression. For this heating method the radial equilibrium depends upon the time evolution of the fast rising magnetic fields outside the plasma. A simple mathematical model for plasma compression and relaxation toward radial equilibrium in the case β = 1 including the external electric circuit is presented in this paper. The numerical solution for different experimental parameters leads to information on values of these parameters for reaching a radial equilibrium. For the plasma compression a snowplow-bounce-model is used which gives the initial values for the equations of the relaxation phase. In this phase the plasma is described by three thin sheaths containing the total mass of the plasma, whose motions are damped by viscosity terms. (orig.) 891 HT/orig. 892 HIS

This thesis reports on three separate investigations in solid state physics. The first is electron paramagnetic resonance in the spin glass Ag:Mn. EPR measurements were performed at two resonance frequencies, concentrating on temperatures above the glass transition temperature. The measured linewidth appears to diverge at T/sub g/ for low resonance frequencies. These results will be compared with recently proposed phenomenological and microscopic theories. The second topic reported in this thesis is the superconducting transition of thin aluminum films. These films were investigated as a function of grain size and thickness. The transition temperature was enhanced over the bulk value, in agreement with many previous investigations of granular aluminum. The third topic reported in this thesis is an extension of the variable rate hopping theory applied in one dimension to N-ME-Qn(TCNQ) 2 . This model is a classical one used to explain both the dc and ac electrical conductivity of organic conductors. The temperature dependence of the model does not agree with experiment at low temperatures. Tunneling has been added to the hopping. This increases the conductivity at low temperatures, and results in excellent agreement with the experimental conductivity over the measured temperature range. The model also predicts that the frequency dependence of the conductivity varies as ω/sup .5/ at low frequencies. This long time tail prediction agrees with the measured dielectric constant of N-Me-iso-Qn(TCNQ) 2

A time-dependent two-fluid model has been developed to understand axial variations in the plasma parameters in a very high density (peak ne≳ 5 ×1019 m-3 ) argon inductively coupled discharge in a long 1.1 cm radius tube. The model equations are written in 1D with radial losses to the tube walls accounted for by the inclusion of effective particle and energy sink terms. The ambipolar diffusion equation and electron energy equation are solved to find the electron density ne(z ,t ) and temperature Te(z ,t ) , and the populations of the neutral argon 4s metastable, 4s resonant, and 4p excited state manifolds are calculated to determine the stepwise ionization rate and calculate radiative energy losses. The model has been validated through comparisons with Langmuir probe ion saturation current measurements; close agreement between the simulated and measured axial plasma density profiles and the initial density rise rate at each location was obtained at pA r=30 -60 mTorr . We present detailed results from calculations at 60 mTorr, including the time-dependent electron temperature, excited state populations, and energy budget within and downstream of the radiofrequency antenna.

We develop a biomorphodynamic model to investigate sediment and vegetation dynamics on a schematic intertidal flat characterized by an initially well-mixed sand-mud mixture. Major interactions between tides, wind waves, salt marshes, sediment transport and sea level rise (SLR) are taken into account. For a bare flat under only tidal action, the model predicts a convex cross-shore profile with the surficial distribution of mud and sand on the upper and lower part of the intertidal flat, respectively. When wind waves are strong, the intertidal flat is highly eroded resulting in a concave profile near the high water mark. This behavior is pronouncedly altered when the intertidal flat is vegetated with the presence of salt marshes. Numerical results suggest that a considerable amount of mud can still remain in the vegetated region even when wave action is strong. A steeper transition zone forms at the boundary between salt marshes and bare flats because of the differential sediment deposition in the two neighboring regions. The inclusion of wind waves is found to considerably enhance the size of the marsh-edge transition zone. For the numerical experiments designed in this study, the profile shape and sediment sorting behavior of tidal flats are not significantly modified by a gradual rising sea level. However, the impacts of SLR on vegetated tidal flats are still manifold: (a) driving the landward migration of intertidal zone and salt marshes; (b) enhancing sediment erosion on intertidal flats; and (c) drowning salt marshes under limited sediment supply with the constrain of seawalls. Finally, model results suggest that organic carbon accumulation on marshlands may be enhanced with an increasing SLR rate provided that salt marshes are not drowned.

We study a 1+1D, stochastic, Burton-Cabrera-Frank (BCF) model of interacting steps fluctuating on a vicinal crystal. The step energy accounts for entropic and nearest-neighbor elastic-dipole interactions. Our goal is to formulate and validate a self-consistent mean-field (MF) formalism to approximately solve the system of coupled, nonlinear stochastic differential equations (SDEs) governing fluctuations in surface motion. We derive formulas for the time-dependent terrace width distribution (TWD) and its steady-state limit. By comparison with kinetic Monte-Carlo simulations, we show that our MF formalism improves upon models in which step interactions are linearized. We also indicate how fitting parameters of our steady state MF TWD may be used to determine the mass transport regime and step interaction energy of certain experimental systems. PP and TLE supported by NSF MRSEC under Grant DMR 05-20471 at U. of Maryland; DM supported by NSF under Grant DMS 08-47587.

The radiation emitted by a plane of excited two-level atoms, Δl=1, Δm=0 transition, is exactly calculated using a semiclassical model in which the electromagnetic field is treated classically (Maxwell's equations) and the coupling between matter and field is described in the electric dipole approximation. The influence of the plane density in the radiation rate is investigated in both limits of weak and strong coupling (defined in the text). It is shown that, in the second case, we can observe infinitely many solutions of the problem, depending on the initial value of the phase difference appearing in the definition of the excited state. Cases of phase choices leading to enhanced or attenuated emission rates are also discussed [pt

Full Text Available The article deals with the development of CFD (Computational Fluid Dynamics model of natural draft wet-cooling tower flow, heat and mass transfer. The moist air flow is described by the system of conservation laws along with additional equations. Moist air is assumed to be homogeneous mixture of dry air and water vapour. Liquid phase in the fill zone is described by the system of ordinary differential equations. Boundary value problem for the system of conservation laws is discretized in space using Kurganov-Tadmor central scheme and in time using strong stability preserving Runge-Kutta scheme. Initial value problems in the fill zone is solved by using standard fourth order Runge-Kutta scheme. The interaction between liquid water and moist air is done by source terms in governing equations.

A 1D TFM numerical simulation of near horizontal stratified two-phase flow is performed where the TFM, including surface tension and viscous stresses, is simplified to a two-equation model using the fixed-flux approximation. As the angle of inclination of the channel increases so does the driving body force, so the flow becomes KH unstable, and waves grow and develop nonlinearities. It is shown that these waves grow until they reach a limit cycle due to viscous dissipation at wave fronts. Upon further inclination of the channel, chaos is observed. The appearance of chaos in a 1D TFM implies a nonlinear process that transfers energy intermittently from long wavelengths where energy is produced to short wavelengths where energy is dissipated by viscosity, so that an averaged energy equilibrium in frequency space is attained. This is comparable to the well-known turbulent stability mechanism of the multi-dimensional Navier–Stokes equations, i.e., chaos implies Lyapunov stability, but in this case it is strictly a two-phase phenomenon.

The dimensional crossover from a 1D fluctuating state at high temperatures to a 3D phase coherent state in the low temperature regime in two coaxial weakly-coupled cylindrical surfaces formed by two-dimensional arrays of parallel nanowires is studied via an 8-state 3D-XY model. This system serves as a model for quasi-one-dimensional superconductors in the form of bundles of weakly-coupled superconducting nanowires. A periodic variation of the dimensional crossover temperature T{sub DC} is observed when the inner superconducting cylindrical surface is rotated in the angular plane. T{sub DC} reaches a maximum when the relative angle between the cylinders is 2.81°, which corresponds to the maximum separation of nanowires between the two cylindrical surfaces. We demonstrate that the relative strength of phase fluctuations in this system is controllable by the rotational angle between the two surfaces with a strong suppression of the fluctuation strength at 2.81°. The phase fluctuations are suppressed gradually upon cooling, before they abruptly vanish below T{sub DC}. Our model thus allows us to study how phase fluctuations can be suppressed in quasi-one-dimensional superconductors in order to achieve a global phase coherent state throughout the nanowire array with zero electric resistance.

In WWTP models, the accurate assessment of solids inventory in bioreactors equipped with solid-liquid separators, mostly described using one-dimensional (1-D) secondary settling tank (SST) models, is the most fundamental requirement of any calibration procedure. Scientific knowledge on characterising particulate organics in wastewater and on bacteria growth is well-established, whereas 1-D SST models and their impact on biomass concentration predictions are still poorly understood. A rigorous assessment of two 1-DSST models is thus presented: one based on hyperbolic (the widely used Takács-model) and one based on parabolic (the more recently presented Plósz-model) partial differential equations. The former model, using numerical approximation to yield realistic behaviour, is currently the most widely used by wastewater treatment process modellers. The latter is a convection-dispersion model that is solved in a numerically sound way. First, the explicit dispersion in the convection-dispersion model and the numerical dispersion for both SST models are calculated. Second, simulation results of effluent suspended solids concentration (XTSS,Eff), sludge recirculation stream (XTSS,RAS) and sludge blanket height (SBH) are used to demonstrate the distinct behaviour of the models. A thorough scenario analysis is carried out using SST feed flow rate, solids concentration, and overflow rate as degrees of freedom, spanning a broad loading spectrum. A comparison between the measurements and the simulation results demonstrates a considerably improved 1-D model realism using the convection-dispersion model in terms of SBH, XTSS,RAS and XTSS,Eff. Third, to assess the propagation of uncertainty derived from settler model structure to the biokinetic model, the impact of the SST model as sub-model in a plant-wide model on the general model performance is evaluated. A long-term simulation of a bulking event is conducted that spans temperature evolution throughout a summer

Full Text Available The African great lakes are of utmost importance for the local economy (fishing, as well as being essential to the survival of the local people. During the past decades, these lakes experienced fast changes in ecosystem structure and functioning, and their future evolution is a major concern. In this study, for the first time a set of one-dimensional lake models are evaluated for Lake Kivu (2.28°S; 28.98°E, East Africa. The unique limnology of this meromictic lake, with the importance of salinity and subsurface springs in a tropical high-altitude climate, presents a worthy challenge to the seven models involved in the Lake Model Intercomparison Project (LakeMIP. Meteorological observations from two automatic weather stations are used to drive the models, whereas a unique dataset, containing over 150 temperature profiles recorded since 2002, is used to assess the model's performance. Simulations are performed over the freshwater layer only (60 m and over the average lake depth (240 m, since salinity increases with depth below 60 m in Lake Kivu and some lake models do not account for the influence of salinity upon lake stratification. All models are able to reproduce the mixing seasonality in Lake Kivu, as well as the magnitude and seasonal cycle of the lake enthalpy change. Differences between the models can be ascribed to variations in the treatment of the radiative forcing and the computation of the turbulent heat fluxes. Fluctuations in wind velocity and solar radiation explain inter-annual variability of observed water column temperatures. The good agreement between the deep simulations and the observed meromictic stratification also shows that a subset of models is able to account for the salinity- and geothermal-induced effects upon deep-water stratification. Finally, based on the strengths and weaknesses discerned in this study, an informed choice of a one-dimensional lake model for a given research purpose becomes possible.

This report summarizes the contents and sources of the photochemical and radiative segment of the LLNL one-dimensional transport-kinetics model of the troposphere and stratosphere. Data include the solar flux incident at the top of the atmosphere, absorption spectra for O 2 , O 3 and NO 2 , and effective absorption coefficients for about 40 photolytic processes as functions of wavelength and, in a few cases, temperature and pressure. The current data set represents understanding of atmospheric photochemical processes as of late 1982 and relies largely on NASA Evaluation Number 5 of Chemical Kinetics and Photochemical Data for Use in Stratospheric Modeling, JPL Publication 82-57 (DeMore et al., 1982). Implementation in the model, including the treatment of multiple scattering and cloud cover, is discussed in Wuebbles (1981)

We consider the Sinai model, in which a random walker moves in a random quenched potential V, and ask the following questions: 1. how can the quenched potential V be inferred from the observations of one or more realizations of the random motion? 2. how many observations (walks) are required to make a reliable inference, that is, to be able to distinguish between two similar but distinct potentials, V 1 and V 2 ? We show how question 1 can be easily solved within the Bayesian framework. In addition, we show that the answer to question 2 is, in general, intimately connected to the calculation of the survival probability of a fictitious walker in a potential W defined from V 1 and V 2 , with partial absorption at sites where V 1 and V 2 do not coincide. For the one-dimensional Sinai model, this survival probability can be analytically calculated, in excellent agreement with numerical simulations.

Full Text Available Diversity plays critical roles in ecosystem functioning, but it remains challenging to model phytoplankton diversity in order to better understand those roles and reproduce consistently observed diversity patterns in the ocean. In contrast to the typical approach of resolving distinct species or functional groups, we present a ContInuous TRAiT-basEd phytoplankton model (CITRATE that focuses on macroscopic system properties such as total biomass, mean trait values, and trait variance. This phytoplankton component is embedded within a nitrogen–phytoplankton-zooplankton–detritus–iron model that itself is coupled with a simplified one-dimensional ocean model. Size is used as the master trait for phytoplankton. CITRATE also incorporates trait diffusion for sustaining diversity and simple representations of physiological acclimation, i.e., flexible chlorophyll-to-carbon and nitrogen-to-carbon ratios. We have implemented CITRATE at two contrasting stations in the North Pacific where several years of observational data are available. The model is driven by physical forcing including vertical eddy diffusivity imported from three-dimensional general ocean circulation models (GCMs. One common set of model parameters for the two stations is optimized using the Delayed-Rejection Adaptive Metropolis–Hasting Monte Carlo (DRAM algorithm. The model faithfully reproduces most of the observed patterns and gives robust predictions on phytoplankton mean size and size diversity. CITRATE is suitable for applications in GCMs and constitutes a prototype upon which more sophisticated continuous trait-based models can be developed.

Diversity plays critical roles in ecosystem functioning, but it remains challenging to model phytoplankton diversity in order to better understand those roles and reproduce consistently observed diversity patterns in the ocean. In contrast to the typical approach of resolving distinct species or functional groups, we present a ContInuous TRAiT-basEd phytoplankton model (CITRATE) that focuses on macroscopic system properties such as total biomass, mean trait values, and trait variance. This phytoplankton component is embedded within a nitrogen-phytoplankton-zooplankton-detritus-iron model that itself is coupled with a simplified one-dimensional ocean model. Size is used as the master trait for phytoplankton. CITRATE also incorporates trait diffusion for sustaining diversity and simple representations of physiological acclimation, i.e., flexible chlorophyll-to-carbon and nitrogen-to-carbon ratios. We have implemented CITRATE at two contrasting stations in the North Pacific where several years of observational data are available. The model is driven by physical forcing including vertical eddy diffusivity imported from three-dimensional general ocean circulation models (GCMs). One common set of model parameters for the two stations is optimized using the Delayed-Rejection Adaptive Metropolis-Hasting Monte Carlo (DRAM) algorithm. The model faithfully reproduces most of the observed patterns and gives robust predictions on phytoplankton mean size and size diversity. CITRATE is suitable for applications in GCMs and constitutes a prototype upon which more sophisticated continuous trait-based models can be developed.

Full Text Available A striking difference between the folk-narrative genres of legend and folktale is how the human characters respond to supernatural, otherworldly, or uncanny beings such as ghosts, gods, dwarves, giants, trolls, talking animals, witches, and fairies. In legend the human actors respond with fear and awe, whereas in folktale they treat such beings as if they were ordinary and unremarkable. Since folktale humans treat all characters as belonging to a single realm, folklorists have described the world of the folktale as one-dimensional, in contrast to the two-dimensionality of the legend. The present investigation examines dimensionality in the third major genre of folk narrative: myth. Using the Greek and Hebrew myths of primordial paradise as sample narratives, the present essay finds—surprisingly—that the humans in these stories respond to the otherworldly one-dimensionally, as folktale characters do, and suggests an explanation for their behavior that is peculiar to the world of myth.

We present a one-dimensional Fickian model that predicts the formation of a double Ga gradient during the fabrication of Cu(In,Ga)Se{sub 2} thin films by three-stage thermal co-evaporation. The model is based on chemical reaction equations, structural data, and effective Ga diffusivities. In the model, the Cu(In,Ga)Se{sub 2} surface is depleted from Ga during the deposition of Cu-Se in the second deposition stage, leading to an accumulation of Ga near the back contact. During the third deposition stage, where In-Ga-Se is deposited at the surface, the atomic fluxes within the growing layer are inverted. This results in the formation of a double Ga gradient within the Cu(In,Ga)Se{sub 2} layer and reproduces experimentally observed Ga distributions. The final shape of the Ga depth profile strongly depends on the temperatures, times and deposition rates used. The model is used to evaluate possible paths to flatten the marked Ga depth profile that is obtained when depositing at low substrate temperatures. We conclude that inserting Ga during the second deposition stage is an effective way to achieve this.

We present a one-dimensional Fickian model that predicts the formation of a double Ga gradient during the fabrication of Cu(In,Ga)Se 2 thin films by three-stage thermal co-evaporation. The model is based on chemical reaction equations, structural data, and effective Ga diffusivities. In the model, the Cu(In,Ga)Se 2 surface is depleted from Ga during the deposition of Cu-Se in the second deposition stage, leading to an accumulation of Ga near the back contact. During the third deposition stage, where In-Ga-Se is deposited at the surface, the atomic fluxes within the growing layer are inverted. This results in the formation of a double Ga gradient within the Cu(In,Ga)Se 2 layer and reproduces experimentally observed Ga distributions. The final shape of the Ga depth profile strongly depends on the temperatures, times and deposition rates used. The model is used to evaluate possible paths to flatten the marked Ga depth profile that is obtained when depositing at low substrate temperatures. We conclude that inserting Ga during the second deposition stage is an effective way to achieve this.

The theory of linear and non-linear infrared response of vibrational Holstein polarons in one-dimensional lattices is presented in order to identify the spectral signatures of self-trapping phenomena. Using a canonical transformation, the optical response is computed from the small polaron point of view which is valid in the anti-adiabatic limit. Two types of phonon baths are considered: optical phonons and acoustical phonons, and simple expressions are derived for the infrared response. It is shown that for the case of optical phonons, the linear response can directly probe the polaron density of states. The model is used to interpret the experimental spectrum of crystalline acetanilide in the C=O range. For the case of acoustical phonons, it is shown that two bound states can be observed in the two-dimensional infrared spectrum at low temperature. At high temperature, analysis of the time-dependence of the two-dimensional infrared spectrum indicates that bath mediated correlations slow down spectral diffusion. The model is used to interpret the experimental linear-spectroscopy of model α-helix and β-sheet polypeptides. This work shows that the Davydov Hamiltonian cannot explain the observations in the NH stretching range.

A radiative-convective equilibrium model of the atmosphere has been coupled with a mixed layer model of the ocean to investigate the response of this one-dimensional system to increasing carbon dioxide amounts in the atmosphere. For global mean conditions a surface temperature rise of about 2 0 K was obtained for a doubling of the carbon dioxide amount, in reasonable agreement with the commonly accepted results of Manabe and Wetherald. This temperature rise was essentially invariant with season and indicates that including a shallow (300 m) ocean slab in this problem does not basically alter previous assessments. While the mixed layer depth of the ocean was only very slightly changed by the temperature increase, which extended throughout the depth of the mixed layer, the impact of this increase on the overall behavior of the ocean warrants further study. A calculation was also made of the temporal variation of the sea surface temperature for three possible carbon dioxide growth rates starting from an initial carbon dioxide content of 300 ppm. This indicated that the thermal inertia of the slab ocean provides a time lag of 8 years in the sea surface temperature response compared to a land situation. This is not considered to be of great significance as regards the likely future climatic impact of carbon dioxide increase

Two-phase flow simulation has been an extended research topic over the years due to the importance of predicting with accuracy the flow behavior within different installations, including nuclear power plants. Some of them are low pressure events, like low water pressure injection, nuclear refueling or natural circulation. This work is devoted to investigate the level of accuracy of the results when a two-phase flow experiment, which has been carried out at low pressure, is performed in a one-dimensional simulation code. In particular, the codes that have been selected to represent the experiment are the best-estimate system codes RELAP5/MOD3 and TRACE v5.0 patch4. The experiment consists in a long vertical pipe along which an air-water fluid in bubbly regime moves upwards in adiabatic conditions and atmospheric pressure. The simulations have been first performed in both codes with their original correlations, which are based on the drift flux model for the case of bubbly regime in vertical pipes. Then, a different implementation for the drag force has been undertaken, in order to perform a simulation with equivalent bubble diameter to the experiment. Results show that the calculation obtained from the codes are within the ranges of validity of the experiment with some discrepancies, which leads to the conclusion that the use of a drag correlation approach is more realistic than drift flux model. (author)

Two-phase flow simulation has been an extended research topic over the years due to the importance of predicting with accuracy the flow behavior within different installations, including nuclear power plants. Some of them are low pressure events, like low water pressure injection, nuclear refueling or natural circulation. This work is devoted to investigate the level of accuracy of the results when a two-phase flow experiment, which has been carried out at low pressure, is performed in a one-dimensional simulation code. In particular, the codes that have been selected to represent the experiment are the best-estimate system codes RELAP5/MOD3 and TRACE v5.0 patch4. The experiment consists in a long vertical pipe along which an air-water fluid in bubbly regime moves upwards in adiabatic conditions and atmospheric pressure. The simulations have been first performed in both codes with their original correlations, which are based on the drift flux model for the case of bubbly regime in vertical pipes. Then, a different implementation for the drag force has been undertaken, in order to perform a simulation with equivalent bubble diameter to the experiment. Results show that the calculation obtained from the codes are within the ranges of validity of the experiment with some discrepancies, which leads to the conclusion that the use of a drag correlation approach is more realistic than drift flux model. (author)

A model of CO2 atmospheric transport in vegetated canopies is tested against measurements of the flow, as well as CO2 concentrations at the Norunda research station located inside a mixed pine-spruce forest. We present the results of simulations of wind-speed profiles and CO2 concentrations inside and above the forest canopy with a one-dimensionalmodel of profiles of the turbulent diffusion coefficient above the canopy accounting for the influence of the roughness sub-layer on turbulent mixing according to Harman and Finnigan (Boundary-Layer Meteorol 129:323-351, 2008; hereafter HF08). Different modelling approaches are used to define the turbulent exchange coefficients for momentum and concentration inside the canopy: (1) the modified HF08 theory—numerical solution of the momentum and concentration equations with a non-constant distribution of leaf area per unit volume; (2) empirical parametrization of the turbulent diffusion coefficient using empirical data concerning the vertical profiles of the Lagrangian time scale and root-mean-square deviation of the vertical velocity component. For neutral, daytime conditions, the second-order turbulence model is also used. The flexibility of the empirical model enables the best fit of the simulated CO2 concentrations inside the canopy to the observations, with the results of simulations for daytime conditions inside the canopy layer only successful provided the respiration fluxes are properly considered. The application of the developed model for radiocarbon atmospheric transport released in the form of ^{14}CO2 is presented and discussed.

Full Text Available This work investigates the influence of meteoric smoke particles (MSP on the charge balance in the D-region ionosphere. Both experimental in situ measurements and a one-dimensional ionospheric model reveal a clear impact of MSP on the ionospheric composition of the D-region. The study reviews rocket-borne in situ measurements of electron and positive ion density, which show a distinct deficit of electrons in comparison to positive ions between 80 and 95 km. This deficit can be explained by the ambient negatively charged MSP measured simultaneously with a Faraday cup. The influence of MSP on the D-region charge balance is addressed with a simplified ionospheric model with only six components, i.e. electrons, positive and negative ions and neutral and charged MSP (both signs. The scheme includes reactions of plasma captured by MSP and MSP photo reactions as well as the standard ionospheric processes, e.g. ion-ion recombination. The model shows that the capture of plasma constituents by MSP is an important process leading to scavenging of electrons. Since Faraday cup measurements are biased towards heavy MSP because of aerodynamical filtering, we have applied an estimate of this filter on the modelled MSP densities. By doing that, we find good qualitative agreement between the experimental data and our model results. In addition, the model study reveals an increase of positive ions in the presence of MSP. That is primarily caused by the reduced dissociative recombination with electrons which have been removed from the gas phase by the MSP.

The present study aims at using statistically designed computational fluid dynamics (CFD) simulations as numerical experiments for the identification of one-dimensional (1-D) advection-dispersion models - computationally light tools, used e.g., as sub-models in systems analysis. The objective is to develop a new 1-D framework, referred to as interpreted CFD (iCFD) models, in which statistical meta-models are used to calculate the pseudo-dispersion coefficient (D) as a function of design and flow boundary conditions. The method - presented in a straightforward and transparent way - is illustrated using the example of a circular secondary settling tank (SST). First, the significant design and flow factors are screened out by applying the statistical method of two-level fractional factorial design of experiments. Second, based on the number of significant factors identified through the factor screening study and system understanding, 50 different sets of design and flow conditions are selected using Latin Hypercube Sampling (LHS). The boundary condition sets are imposed on a 2-D axi-symmetrical CFD simulation model of the SST. In the framework, to degenerate the 2-D model structure, CFD model outputs are approximated by the 1-D model through the calibration of three different model structures for D. Correlation equations for the D parameter then are identified as a function of the selected design and flow boundary conditions (meta-models), and their accuracy is evaluated against D values estimated in each numerical experiment. The evaluation and validation of the iCFD model structure is carried out using scenario simulation results obtained with parameters sampled from the corners of the LHS experimental region. For the studied SST, additional iCFD model development was carried out in terms of (i) assessing different density current sub-models; (ii) implementation of a combined flocculation, hindered, transient and compression settling velocity function; and (iii

The post-arc dielectric recovery process has a decisive effect on the current interruption performance in a vacuum circuit breaker. The dissipation of residual plasma at the moment of current zero under the transient recovery voltage, which is the first stage of the post-arc dielectric recovery process and forms the post-arc current, has attracted many concerns. A one-dimensional particle-in-cell model is developed to simulate the measured post-arc current in the vacuum circuit breaker in this paper. At first, the parameters of the residual plasma are estimated roughly by the waveform of the post-arc current which is taken from measurements. After that, different components of the post-arc current, which are formed by the movement of charged particles in the residual plasma, are discussed. Then, the residual plasma density is adjusted according to the proportion of electrons and ions absorbed by the post-arc anode derived from the particle-in-cell simulation. After this adjustment, the post-arc current waveform obtained from the simulation is closer to that obtained from measurements.

Potential formation in a bounded plasma system that contains electrons with a two-temperature velocity distribution and is terminated by a floating, electron emitting electrode (collector) is studied by a one-dimensional kinetic model. A method on how to determine the boundary conditions at the collector for the numerical solution of the Poisson equation is presented. The difference between the regular and the irregular numerical solutions of the Poisson equation is explained. The regular numerical solution of the Poisson equation fulfills the boundary conditions at the source and can be computed for any distance from the collector. The irregular solution does not fulfill the source boundary conditions and the computation breaks down at some distance from the collector. An excellent agreement of the values of the potential at the inflection point found from the numerical solution of the Poisson equation with the values predicted by the analytical model is obtained. Potential, electric field, and particle density profiles found by the numerical solution of the Poisson equation are compared to the profiles obtained with the particle in cell computer simulation. A very good quantitative agreement of the potential and electric field profiles is obtained. For certain values of the parameters the analytical model predicts three possible values of the potential at the inflection point. In such cases always only one of the corresponding numerical solutions of the Poisson equation is regular, while the other two are irregular. The regular numerical solution of the Poisson equation always corresponds to the solution of the model that predicts the largest ion flux to the collector

A one-dimensional, two-fluid, steady state model is used for the analysis of ion temperature effects to the plasma-wall transition. In this paper, the model is solved for a finite ratio ɛ between the Debye and the ionization length, while in Part II [T. Gyergyek and J. Kovačič, Phys Plasmas 24, 063506 (2017)], the solutions for ɛ = 0 are presented. Ion temperature is treated as a given, independent parameter and it is included in the model as a boundary condition. It is shown that when the ion temperature larger than zero is selected, the ion flow velocity and the electric field at the boundary must be consistent with the selected ion temperature. A numerical procedure, how to determine such "consistent boundary conditions," is proposed, and a simple relation between the ion temperature and ion velocity at the boundary of the system is found. The effects of the ion temperature to the pre-sheath length, potential, ion temperature, and ion density drops in the pre-sheath and in the sheath are investigated. It is concluded that larger ion temperature results in a better shielding of the plasma from the wall. An attempt is made to include the ion heat flux qi into the model in its simplest form q i = - K ' /d T i d x , where K ' is a constant heat conduction coefficient. It is shown that inclusion of such a term into the energy transfer equation introduces an additional ion heating mechanism into the system and the ion flow then becomes isothermal instead of adiabatic even in the sheath.

Full Text Available Existing theoretical formulations for the size-resolved scavenging coefficient Λ(d for atmospheric aerosol particles scavenged by rain predict values lower by one to two orders of magnitude than those estimated from field measurements of particle-concentration changes for particles smaller than 3 μm in diameter. Vertical turbulence is not accounted for in the theoretical formulations of Λ(d but does contribute to the field-derived estimates of Λ(d due to its influence on the overall concentration changes of aerosol particles in the layers undergoing impaction scavenging. A detailed one-dimensional cloud microphysics model has been used to simulate rain production and below-cloud particle scavenging, and to quantify the contribution of turbulent diffusion to the overall Λ(d values calculated from particle concentration changes. The relative contribution of vertical diffusion to below-cloud scavenging is found to be largest for submicron particles under weak precipitation conditions. The discrepancies between theoretical and field-derived Λ(d values can largely be explained by the contribution of vertical diffusion to below-cloud particle scavenging for all particles larger than 0.01 μm in diameter for which field data are available. The results presented here suggest that the current theoretical framework for Λ(d can provide a reasonable approximation of below-cloud aerosol particle scavenging by rain in size-resolved aerosol transport models if vertical diffusion is also considered by the models.

Since the inception of the Antarctic ice sheet at the Eocene-Oligocene transition (˜ 34 Myr ago), land ice has played a crucial role in Earth's climate. Through feedbacks in the climate system, land ice variability modifies atmospheric temperature changes induced by orbital, topographical, and greenhouse gas variations. Quantification of these feedbacks on long timescales has hitherto scarcely been undertaken. In this study, we use a zonally averaged energy balance climate model bidirectionally coupled to a one-dimensional ice sheet model, capturing the ice-albedo and surface-height-temperature feedbacks. Potentially important transient changes in topographic boundary conditions by tectonics and erosion are not taken into account but are briefly discussed. The relative simplicity of the coupled model allows us to perform integrations over the past 38 Myr in a fully transient fashion using a benthic oxygen isotope record as forcing to inversely simulate CO2. Firstly, we find that the results of the simulations over the past 5 Myr are dependent on whether the model run is started at 5 or 38 Myr ago. This is because the relation between CO2 and temperature is subject to hysteresis. When the climate cools from very high CO2 levels, as in the longer transient 38 Myr run, temperatures in the lower CO2 range of the past 5 Myr are higher than when the climate is initialised at low temperatures. Consequently, the modelled CO2 concentrations depend on the initial state. Taking the realistic warm initialisation into account, we come to a best estimate of CO2, temperature, ice-volume-equivalent sea level, and benthic δ18O over the past 38 Myr. Secondly, we study the influence of ice sheets on the evolution of global temperature and polar amplification by comparing runs with ice sheet-climate interaction switched on and off. By passing only albedo or surface height changes to the climate model, we can distinguish the separate effects of the ice-albedo and surface

Surface stones affect erosion rates by reducing raindrop-driven detachment and protecting the original soil against overland flow induced-hydraulic stress. Numerous studies have shown that the effect of surface stones on erosion depends on both the stone characteristics (e.g., size, distribution) and the soil properties. The aim of this study was (i) to quantify how the stone characteristics can affect the total sediment concentration and the concentrations of the individual size classes, (ii) to test if stones affect preferentially a particular size class within the eroded sediment and (iii) to determine whether the 1D Hairsine-Rose (H-R) erosion model can represent the experimental data. A series of laboratory experiments were conducted using the 2 m × 6 m EPFL erosion flume for a high rainfall intensity (60 mm/h) event on a gentle slope (2.2%). The flume was divided into two identical 1-m wide flumes. This separation was done to allow simultaneous replicate experiments. Experiments were conducted with different configurations and scenarios (stone coverage, size and emplacement). Three coverage proportions (20%, 40%, and 70%), two stone diameters (3-4 and 6-7 cm) and two emplacement types (topsoil and partially embedded) were tested. For each experiment, the total sediment concentration, the concentration for the individual size classes, and the flume discharge were measured. Infiltration rates were measured at different depths and locations. A high resolution laser scanner provided details of the surface change due to erosion during the experiments. This technique allowed us to quantify the spatial distribution of eroded soil and to understand better if sediment transport is 1D or rather 2D over the flumes. The one-dimensional Hairsine-Rose (H-R) erosion model was used to fit the integrated data and to provide estimates of the parameters. The ability of the 1D H-R model to predict the measured sediment concentrations in the presence of stones in the soil matrix

We study the quantum-statistical mechanics, at low temperatures, of a one-dimensional antiferromagnetic Heisenberg model with two anisotropies. In the weak-coupling limit we determine the temperature dependences of the soliton energy and the soliton density. We have found that the leading correction to the sine-Gordon (SG) expression for the soliton density and the quantum soliton energy comes from the out-of-plane magnon mode, not present in the pure SG model. We also show that when an external magnetic field is applied, the chain supports a new type of kink, where the sublattices rotate in opposite directions

The Full EQuations (FEQ) model is a computer program for solution of the full, dynamic equations of motion for one-dimensional unsteady flow in open channels and through control structures. A stream system that is simulated by application of FEQ is subdivided into stream reaches (branches), parts of the stream system for which complete information on flow and depth are not required (dummy branches), and level-pool reservoirs. These components are connected by special features; that is, hydraulic control structures, including junctions, bridges, culverts, dams, waterfalls, spillways, weirs, side weirs, and pumps. The principles of conservation of mass and conservation of momentum are used to calculate the flow and depth throughout the stream system resulting from known initial and boundary conditions by means of an implicit finite-difference approximation at fixed points (computational nodes). The hydraulic characteristics of (1) branches including top width, area, first moment of area with respect to the water surface, conveyance, and flux coefficients and (2) special features (relations between flow and headwater and (or) tail-water elevations, including the operation of variable-geometry structures) are stored in function tables calculated in the companion program, Full EQuations UTiLities (FEQUTL). Function tables containing other information used in unsteady-flow simulation (boundary conditions, tributary inflows or outflows, gate settings, correction factors, characteristics of dummy branches and level-pool reservoirs, and wind speed and direction) are prepared by the user as detailed in this report. In the iterative solution scheme for flow and depth throughout the stream system, an interpolation of the function tables corresponding to the computational nodes throughout the stream system is done in the model. FEQ can be applied in the simulation of a wide range of stream configurations (including loops), lateral-inflow conditions, and special features. The

Full Text Available Since the inception of the Antarctic ice sheet at the Eocene–Oligocene transition (∼ 34 Myr ago, land ice has played a crucial role in Earth's climate. Through feedbacks in the climate system, land ice variability modifies atmospheric temperature changes induced by orbital, topographical, and greenhouse gas variations. Quantification of these feedbacks on long timescales has hitherto scarcely been undertaken. In this study, we use a zonally averaged energy balance climate model bidirectionally coupled to a one-dimensional ice sheet model, capturing the ice–albedo and surface–height–temperature feedbacks. Potentially important transient changes in topographic boundary conditions by tectonics and erosion are not taken into account but are briefly discussed. The relative simplicity of the coupled model allows us to perform integrations over the past 38 Myr in a fully transient fashion using a benthic oxygen isotope record as forcing to inversely simulate CO2. Firstly, we find that the results of the simulations over the past 5 Myr are dependent on whether the model run is started at 5 or 38 Myr ago. This is because the relation between CO2 and temperature is subject to hysteresis. When the climate cools from very high CO2 levels, as in the longer transient 38 Myr run, temperatures in the lower CO2 range of the past 5 Myr are higher than when the climate is initialised at low temperatures. Consequently, the modelled CO2 concentrations depend on the initial state. Taking the realistic warm initialisation into account, we come to a best estimate of CO2, temperature, ice-volume-equivalent sea level, and benthic δ18O over the past 38 Myr. Secondly, we study the influence of ice sheets on the evolution of global temperature and polar amplification by comparing runs with ice sheet–climate interaction switched on and off. By passing only albedo or surface height changes to the climate model, we can distinguish the

This paper discusses some onedimensional quantum models with singular perturbations. Eventually, a mass discontinuity is added at the points that support the singular perturbations. The simplest model includes an attractive singular potential with a mass jump both located at the origin. We study the form of the only bound state. Another model exhibits a hard core at the origin plus one or more repulsive deltas with mass jumps at the points supporting these deltas. We study the location and the multiplicity of these resonances for the case of one or two deltas and settle the basis for a generalization. Finally, we consider the harmonic oscillator and the infinite square well plus a singular potential at the origin. We see how the energy of bound states is affected by the singular perturbation.

In this work we generalize the solution of the one-dimensional neutron transport equation to a multi- group approach in planar geometry. The basic idea of this work consists in consider the hierarchical construction of a solution for a generic number G of energy groups, starting from a mono-energetic solution. The hierarchical method follows the reasoning of the decomposition method. More specifically, the additional terms from adding energy groups is incorporated into the recursive scheme as source terms. This procedure leads to an analytical representation for the solution with G energy groups. The recursion depth is related to the accuracy of the solution, that may be evaluated after each recursion step. The authors present a heuristic analysis of stability for the results. Numerical simulations for a specific example with four energy groups and a localized pulsed source. (author)

A one-dimensional rule-based model for flocking, which combines velocity alignment and long-range centering interactions, is presented and studied. The induced cohesion in the collective motion of the self-propelled agents leads to unique group behavior that contrasts with previous studies. Our results show that the largest cluster of particles, in the condensed states, develops a mean velocity slower than the preferred one in the absence of noise. For strong noise, the system also develops a non-vanishing mean velocity, alternating its direction of motion stochastically. This allows us to address the directional switching phenomenon. The effects of different sources of stochasticity on the system are also discussed. (paper)

Largely nonmathematical qualitative lectures are given on the basic physics of nearly one-dimensional conductors. The main emphasis is placed on the properties of a purely one-dimensional electron gas. The effects of a real system having interchain coupling, impurities, a compressible lattice, lattice distortions and phonon anomalies are discussed

There is currently a great interest in the physics of degenerate quantum gases and low-energy few-body scattering due to the recent experimental advances in manipulation of ultracold atoms by light. In particular, almost perfect periodic potentials, called optical lattices, can be generated. The lattice spacing is fixed by the wavelength of the laser field employed and the angle betwen the pair of laser beams; the lattice depth, defining the magnitude of the different band gaps, is tunable within a large interval of values. This flexibility permits the exploration of different regimes, ranging from the ''free-electron'' picture, modified by the effective mass for shallow optical lattices, to the tight-binding regime of a very deep periodic potential. In the latter case, effective single-band theories, widely used in condensed matter physics, can be implemented with unprecedent accuracy. The tunability of the lattice depth is nowadays complemented by the use of magnetic Feshbach resonances which, at very low temperatures, can vary the relevant atom-atom scattering properties at will. Moreover, optical lattices loaded with gases of effectively reduced dimensionality are experimentally accessible. This is especially important for one spatial dimension, since most of the exactly solvable models in many-body quantum mechanics deal with particles on a line; therefore, experiments with one-dimensional gases serve as a testing ground for many old and new theories which were regarded as purely academic not so long ago. The physics of few quantum particles on a one-dimensional lattice is the topic of this thesis. Most of the results are obtained in the tight-binding approximation, which is amenable to exact numerical or analytical treatment. For the two-body problem, theoretical methods for calculating the stationary scattering and bound states are developed. These are used to obtain, in closed form, the two-particle solutions of both the Hubbard and extended Hubbard models

There is currently a great interest in the physics of degenerate quantum gases and low-energy few-body scattering due to the recent experimental advances in manipulation of ultracold atoms by light. In particular, almost perfect periodic potentials, called optical lattices, can be generated. The lattice spacing is fixed by the wavelength of the laser field employed and the angle betwen the pair of laser beams; the lattice depth, defining the magnitude of the different band gaps, is tunable within a large interval of values. This flexibility permits the exploration of different regimes, ranging from the ''free-electron'' picture, modified by the effective mass for shallow optical lattices, to the tight-binding regime of a very deep periodic potential. In the latter case, effective single-band theories, widely used in condensed matter physics, can be implemented with unprecedent accuracy. The tunability of the lattice depth is nowadays complemented by the use of magnetic Feshbach resonances which, at very low temperatures, can vary the relevant atom-atom scattering properties at will. Moreover, optical lattices loaded with gases of effectively reduced dimensionality are experimentally accessible. This is especially important for one spatial dimension, since most of the exactly solvable models in many-body quantum mechanics deal with particles on a line; therefore, experiments with one-dimensional gases serve as a testing ground for many old and new theories which were regarded as purely academic not so long ago. The physics of few quantum particles on a one-dimensional lattice is the topic of this thesis. Most of the results are obtained in the tight-binding approximation, which is amenable to exact numerical or analytical treatment. For the two-body problem, theoretical methods for calculating the stationary scattering and bound states are developed. These are used to obtain, in closed form, the two-particle solutions of both the Hubbard and

The physics of one-dimensional interacting bosonic systems is reviewed. Beginning with results from exactly solvable models and computational approaches, the concept of bosonic Tomonaga-Luttinger liquids relevant for one-dimensional Bose fluids is introduced, and compared with Bose-Einstein condensates existing in dimensions higher than one. The effects of various perturbations on the Tomonaga-Luttinger liquid state are discussed as well as extensions to multicomponent and out of equilibrium ...

The notion of an island has surfaced in recent algebra and coding theory research. Discrete versions provide interesting combinatorial problems. This paper presents the one-dimensional case with finitely many heights, a topic convenient for student research.

A class of one-dimensional classical systems is characterized from an algebraic point of view. The Hamiltonians of these systems are factorized in terms of two functions that together with the Hamiltonian itself close a Poisson algebra. These two functions lead directly to two time-dependent integrals of motion from which the phase motions are derived algebraically. The systems so obtained constitute the classical analogues of the well known factorizable one-dimensional quantum mechanical systems

In this article we present a method to determine the band spectrum, band gaps, and discrete energy levels, of a one-dimensional photonic crystal with localized impurities. For one-dimensional crystals with piecewise constant refractive indices we develop an algorithm to recover the refractive index distribution from the period map. Finally, we derive the relationship between the period map and the scattering matrix containing the information on the localized modes.

Here, the diffusiophoretic motion of suspended colloidal particles under one-dimensional solute gradients is solved using numerical and analytical techniques. Similarity solutions are developed for the injection and withdrawal dynamics of particles into semi-infinite pores. Furthermore, a method of characteristics formulation of the diffusion-free particle transport model is presented and integrated to realize particle trajectories. Analytical solutions are presented for the limit of small particle diffusiophoretic mobility Γp relative to the solute diffusivity Ds for particle motions in both semi-infinite and finite domains. Results confirm the build up of local maxima and minima in the propagating particle front dynamics. The method of characteristics is shown to successfully predict particle motions and the position of the particle front, although it fails to accurately predict suspended particle concentrations in the vicinity of sharp gradients, such as at the particle front peak seen in some injection cases, where particle diffusion inevitably plays an important role. Results inform the design of applications in which the use of applied solute gradients can greatly enhance particle injection into and withdrawal from pores.

Here, the diffusiophoretic motion of suspended colloidal particles under one-dimensional solute gradients is solved using numerical and analytical techniques. Similarity solutions are developed for the injection and withdrawal dynamics of particles into semi-infinite pores. Furthermore, a method of characteristics formulation of the diffusion-free particle transport model is presented and integrated to realize particle trajectories. Analytical solutions are presented for the limit of small particle diffusiophoretic mobility Γ p relative to the solute diffusivity D s for particle motions in both semi-infinite and finite domains. Results confirm the build up of local maxima and minima in the propagating particle front dynamics. The method of characteristics is shown to successfully predict particle motions and the position of the particle front, although it fails to accurately predict suspended particle concentrations in the vicinity of sharp gradients, such as at the particle front peak seen in some injection cases, where particle diffusion inevitably plays an important role. Results inform the design of applications in which the use of applied solute gradients can greatly enhance particle injection into and withdrawal from pores.

In this article we have used the 2D fluid turbulence numerical model, ESEL, to simulate turbulent transport in edge tokamak plasma. Basic plasma parameters from the ASDEX Upgrade and COMPASS tokamaks are used as input for the model, and the output is compared with experimental observations obtain...... for an extension of the ESEL model from 2D to 3D to fully resolve the parallel dynamics, and the coupling from the plasma to the sheath....

Full Text Available The helicopter-borne electromagnetic (HEM frequency-domain exploration method is an airborne electromagnetic (AEM technique that is widely used for vast and rough areas for resistivity imaging. The vast amount of digitized data flowing from the HEM method requires an efficient and accurate inversion algorithm. Generally, the inverse modelling of HEM data in the first step requires a precise and efficient technique provided by a forward modelling algorithm. The exact calculation of the sensitivity matrix or Jacobian is also of the utmost importance. As such, the main objective of this study is to design an efficient algorithm for the forward modelling of HEM frequency-domain data for the configuration of horizontal coplanar (HCP coils using fast Hankel transforms (FHTs. An attempt is also made to use an analytical approach to derive the required equations for the Jacobian matrix. To achieve these goals, an elaborated algorithm for the simultaneous calculation of the forward computation and sensitivity matrix is provided. Finally, using two synthetic models, the accuracy of the calculations of the proposed algorithm is verified. A comparison indicates that the obtained results of forward modelling are highly consistent with those reported in Simon et al. (2009 for a four-layer model. Furthermore, the comparison of the results for the sensitivity matrix for a two-layer model with those obtained from software is being used by the BGR Centre in Germany, showing that the proposed algorithm enjoys a high degree of accuracy in calculating this matrix.

using the example of a circular secondary settling tank (SST). First, the significant design and flow factors are screened out by applying the statistical method of two-level fractional factorial design of experiments. Second, based on the number of significant factors identified through the factor...... of (i) assessing different density current sub-models; (ii) implementation of a combined flocculation, hindered, transient and compression settling velocity function; and (iii) assessment of modelling the onset of transient and compression settling. Furthermore, the optimal level of model discretization...

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

Motivated by possible generalizations to more complex multiphase multicomponent systems in higher dimensions, we develop an Eulerian-Lagrangian numerical approximation for a system of two conservation laws in one space dimension modeling a

The paper is devoted to a new branch in the theory of one-dimensional variational problems with branching extremals, the investigation of one-dimensional minimal fillings introduced by the authors. On the one hand, this problem is a one-dimensional version of a generalization of Gromov's minimal fillings problem to the case of stratified manifolds. On the other hand, this problem is interesting in itself and also can be considered as a generalization of another classical problem, the Steiner problem on the construction of a shortest network connecting a given set of terminals. Besides the statement of the problem, we discuss several properties of the minimal fillings and state several conjectures. Bibliography: 38 titles.

Full Text Available The present work was a study on global reaction rate of methanol synthesis. We measured experimentally the global reaction rate in the internal recycle gradientless reactor over catalyst SC309. The diffusion-reaction model of methanol synthesis was suggested. For model we chose the hydrogenation of CO and CO2 as key reaction. CO and CO2 were key components in our model. The internal diffusion effectiveness factors of CO and CO2 in the catalyst were calculated by the numerical integration. A comparison with the experiment showed that all the absolute values of the relative error were less than 10%. The simulation results showed that decreasing reaction temperature and catalyst diameter were conducive to reduce the influence of the internal diffusion on the methanol synthesis.

This paper describes a heat transfer model to be implemented in a global engine 1-D gas-dynamic code to calculate reciprocating internal combustion engine performance in steady and transient operations. A trade off between simplicity and accuracy has been looked for, in order to fit with the stated objective. To validate the model, the temperature of the exhaust manifold wall in a high-speed direct injection (HSDI) turbocharged diesel engine has been measured during a full load transient. In addition, an indirect assessment of the exhaust gas temperature during this transient process has been carried out. The results show good agreement between the measured and modelled data with good accuracy to predict the engine performance. A dual-walled air gap exhaust manifold has been tested in order to quantify the potential of exhaust gas thermal energy saving on engine transient performance. The experimental results together with the heat transfer model have been used to analyse the influence of thermal energy saving on dynamic performance during the load transient of an HSDI turbocharged diesel engine. (author)

From both the theoretical and the practical points of view, the problem of constitutive laws is a part and parcel of the modeling problem. In particular, the necessity to restore in the model, through topological laws, some of the information lost during the usual averaging process is emphasized. It is shown that the customary 'void fraction' topological law Psub(V)=Psub(L) should be proscribed whenever propagation phenomena are involved. A new void fraction topological law is proposed. The limitations of the current assumption of constant pressure within any phase in any cross section are also illustrated. The importance of proximity effects (neighborhood and history effects, related to characteristic lengths and times) is brought out. It results in the importance of the mathematical form of the constitutive laws. Various approaches to the constitutive law problem and possible mathematical forms for the transfer laws are reviewed. The simplest form (transfert terms as functions of the dependent variables only) may have some usefulness if interpretation of the results in terms of propagation phenomena is banned. A good compromise between the necessity to take proximity effects into account and to obtain a tractable set of equations is carried out when so called 'differential terms' are introduced in the transfer laws. The last part of the paper is devoted to some restrictions, which are imposed to the transfer terms because of some basic principles: indifference to Galilean changes of frame and to some changes of origins, second law of thermodynamics and assumption of local thermodynamic equilibrium, closure constraints. Practical recommendations are formulated [fr

The temperature variations of first-, second-, and third-sound velocity and attenuation coefficients in one-dimensional superfluid helium are evaluated explicitly for very low temperatures and frequencies (ω/sub s/tau 2 , and the ratio of second sound to first sound becomes unity as the temperature decreases to absolute zero

Full Text Available We study the equilibrium statistical mechanics of simple fluids in narrow pores. A systematic expansion is made about a one-dimensional limit of this system. It starts with a density functional, constructed from projected densities, which depends upon projected one and two-body potentials. The nature of higher order corrections is discussed.

Although the problem of a metal in one dimension has long been known to solid-state physicists, it was not until the synthesis of real one-dimensional or quasi-one-dimensional systems that this subject began to attract considerable attention. This has been due in part to the search for high­ temperature superconductivity and the possibility of reaching this goal with quasi-one-dimensional substances. A period of intense activity began in 1973 with the report of a measurement of an apparently divergent conduc­ tivity peak in TfF-TCNQ. Since then a great deal has been learned about quasi-one-dimensional conductors. The emphasis now has shifted from trying to find materials of very high conductivity to the many interesting problems of physics and chemistry involved. But many questions remain open and are still under active investigation. This book gives a review of the experimental as well as theoretical progress made in this field over the last years. All the chapters have been written by scientists who have ...

Full Text Available In this paper we study a class of one-dimensional Dirichlet boundary value problems involving the Caputo fractional derivatives. The existence of infinitely many solutions for this equations is obtained by exploiting a recent abstract result. Concrete examples of applications are presented.

An approximate method for determining the characteristics associated with one-dimensional particulate two-phase flow models is presented. The method is based on iteration and is valid for small particulate volume fractions. The method is applied to several special cases involving incompressible particles suspended in a gas. The influences of certain changes in the physical model are investigated

The Collaborative Modeling for Parametric Assessment of Space Systems (COMPASS) team at Glenn Research Center has performed integrated system analysis of conceptual spacecraft mission designs since 2006 using a multidisciplinary concurrent engineering process. The set of completed designs was archived in a database, to allow for the study of relationships between design parameters. Although COMPASS uses a parametric spacecraft costing model, this research investigated the possibility of using a top-down approach to rapidly estimate the overall vehicle costs. This paper presents the relationships between significant design variables, including breakdowns of dry mass, wet mass, and cost. It also develops a model for a broad estimate of these parameters through basic mission characteristics, including the target location distance, the payload mass, the duration, the delta-v requirement, and the type of mission, propulsion, and electrical power. Finally, this paper examines the accuracy of this model in regards to past COMPASS designs, with an assessment of outlying spacecraft, and compares the results to historical data of completed NASA missions.

We present a first-principles band calculation for the quasi-one-dimensional (Q1D) organic superconductor (TMTSF) 2ClO4 . An effective tight-binding model with the TMTSF molecule to be regarded as the site is derived from a calculation based on maximally localized Wannier orbitals. We apply a two-particle self-consistent (TPSC) analysis by using a four-site Hubbard model, which is composed of the tight-binding model and an onsite (intramolecular) repulsive interaction, which serves as a variable parameter. We assume that the pairing mechanism is mediated by the spin fluctuation, and the sign of the superconducting gap changes between the inner and outer Fermi surfaces, which correspond to a d -wave gap function in a simplified Q1D model. With the parameters we adopt, the critical temperature for superconductivity estimated by the TPSC approach is approximately 1 K, which is consistent with experiment.

experiments show strong signatures of beam steering that are dependent upon graphene doping. This seed grant has allowed our team to establish the essential...based, one-dimensional reflectarrays. Several immediate improvements to the device design and process flow are essential to suppress specular...beam steering that are dependent upon graphene doping. This seed grant has allowed our team to establish the essential operating procedures (i.e

The zero-temperature correlation functions of the one-dimensional attractive Bose gas with a delta-function interaction are calculated analytically for any value of the interaction parameter and number of particles, directly from the integrability of the model. We point out a number of interesting features, including zero recoil energy for a large number of particles, analogous to the Moessbauer effect

This paper brings into focus the results of numerous studies that have addressed issues associated with the validity of assumptions which are used to justify reducing the dimensionality of numerical calculations of water flow through Yucca Mountain, NV. It is shown that in many cases, one-dimensionalmodeling is more rigorous than previously assumed

This report brings into focus the results of numerous studies that have addressed issues associated with the validity of assumptions which are used to justify reducing the dimensionality of numerical calculations of water flow through Yucca Mountain, NV. it is shown that, in many cases, one-dimensionalmodeling is more rigorous than previously assumed

In the companion to the present paper, the one-dimensional k-ɛ lake model SIMSTRAT is coupled to a single-column atmospheric model, nicknamed FIZC, and an application of the coupled model to the deep Lake Geneva, Switzerland, is described. In this paper, the response of Lake Geneva to global warming caused by an increase in atmospheric carbon dioxide concentration (i.e., 2 × CO2) is investigated. Coupling the models allowed for feedbacks between the lake surface and the atmosphere and produced changes in atmospheric moisture and cloud cover that further modified the downward radiation fluxes. The time evolution of atmospheric variables as well as those of the lake's thermal profile could be reproduced realistically by devising a set of adjustable parameters. In a "control" 1 × CO2 climate experiment, the coupled FIZC-SIMSTRAT model demonstrated genuine skills in reproducing epilimnetic and hypolimnetic temperatures, with annual mean errors and standard deviations of 0.25°C ± 0.25°C and 0.3°C ± 0.15°C, respectively. Doubling the CO2 concentration induced an atmospheric warming that impacted the lake's thermal structure, increasing the stability of the water column and extending the stratified period by 3 weeks. Epilimnetic temperatures were seen to increase by 2.6°C to 4.2°C, while hypolimnion temperatures increased by 2.2°C. Climate change modified components of the surface energy budget through changes mainly in air temperature, moisture, and cloud cover. During summer, reduced cloud cover resulted in an increase in the annual net solar radiation budget. A larger water vapor deficit at the air-water interface induced a cooling effect in the lake.

Formulation was made on a one-dimensional beam finite element which is effective in analyzing structural response of very large floating structures by modeling them on beams on an elastic foundation. This element allows strict solution of vibration response in the beams on the elastic foundation to be calculated efficiently for a case where mass and rigidity change in the longitudinal direction. This analysis method was used to analyze structural response of a large pontoon-type floating structure to investigate mass in the end part for the structural response and the effect of decay while passing the structure. With a pontoon-type floating structure, reduction in bends and bending stress in the end part of the floating structure is important in designing the structure. Reducing the mass in the end part is effective as a means to avoid resonance in these responses and reduce the responses. Increase in rigidity of a floating structure shifts the peak in quasi-static response to lower frequency side, and reduces response in resonance, hence it is advantageous for improving the response. Since incident waves decay while passing through the floating structure, response in the lower wave side decreases. The peak frequency in the quasi-static response also decreases at the end part of the structure in the upper wave side due to decay in wave force. 7 refs., 11 figs., 1 tab.

Formulation was made on a one-dimensional beam finite element which is effective in analyzing structural response of very large floating structures by modeling them on beams on an elastic foundation. This element allows strict solution of vibration response in the beams on the elastic foundation to be calculated efficiently for a case where mass and rigidity change in the longitudinal direction. This analysis method was used to analyze structural response of a large pontoon-type floating structure to investigate mass in the end part for the structural response and the effect of decay while passing the structure. With a pontoon-type floating structure, reduction in bends and bending stress in the end part of the floating structure is important in designing the structure. Reducing the mass in the end part is effective as a means to avoid resonance in these responses and reduce the responses. Increase in rigidity of a floating structure shifts the peak in quasi-static response to lower frequency side, and reduces response in resonance, hence it is advantageous for improving the response. Since incident waves decay while passing through the floating structure, response in the lower wave side decreases. The peak frequency in the quasi-static response also decreases at the end part of the structure in the upper wave side due to decay in wave force. 7 refs., 11 figs., 1 tab.

It is shown that in the framework of the boundary condition models (BCM) for the two-particle interaction the Schroedinger equation for the system of three identical bosons can be reduced to the one-dimensional integral equation in an exact way. The method used for obtaining such an equation is based on a special consideration of the two-particle off-shell wave functions. The binding energy of the simple three-particle system is calculated. It is indicated that by means of the equation obtained it is possible to change the off-shell behaviour of the two-particle t-matrix and therefore to simulate three particle effects. (Auth.)

Some basic plasma phenomena are studied by a one-dimensional electrostatic simulation code. A brief description of the code and its application to a test problem is given. The experiments carried out include Landau damping of an excited wave, particle retardation by smoothed field and beam-plasma instability. In each case, the set-up of the experiment is described and the results are compared with theoretical predictions. In the theoretical discussions, the oscillatory behaviour found in the Landau damping experiment is explained, an explicit formula for the particle retardation rate is derived and a rudimentary picture of the beam-plasma instability in terms of quasilinear theory is given. (author)

We report experimental observation of a normal incidence phononic band gap in one-dimensional periodic (SiO(2)/poly(methyl methacrylate)) multilayer film at gigahertz frequencies using Brillouin spectroscopy. The band gap to midgap ratio of 0.30 occurs for elastic wave propagation along the periodicity direction, whereas for inplane propagation the system displays an effective medium behavior. The phononic properties are well captured by numerical simulations. The porosity in the silica layers presents a structural scaffold for the introduction of secondary active media for potential coupling between phonons and other excitations, such as photons and electrons.

One-dimensional dynamics of a plane slab of cold (β ≪ 1) isothermal plasma accelerated by a magnetic field is studied in terms of the MHD equations with a finite constant conductivity. The passage to the limit β → 0 is analyzed in detail. It is shown that, at β = 0, the character of the solution depends substantially on the boundary condition for the electric field at the inner plasma boundary. The relationship between the boundary condition for the pressure at β > 0 and the conditions for the electric field at β = 0 is found. The stability of the solution against one-dimensional longitudinal perturbations is analyzed. It is shown that, in the limit β → 0, the stationary solution is unstable if the time during which the acoustic wave propagates across the slab is longer than the time of magnetic field diffusion. The growth rate and threshold of instability are determined, and results of numerical simulation of its nonlinear stage are presented.

The search for higher energy density, safer, and longer cycling-life energy storage systems is progressing quickly. One-dimensional (1D) nanomaterials have a large length-to-diameter ratio, resulting in their unique electrical, mechanical, magnetic and chemical properties, and have wide applications as electrode materials in different systems. This article reviews the latest hot topics in applying 1D nanomaterials, covering both their synthesis and their applications. 1D nanomaterials can be grouped into the categories: carbon, silicon, metal oxides, and conducting polymers, and we structure our discussion accordingly. Then, we survey the unique properties and application of 1D nanomaterials in batteries and supercapacitors, and provide comments on the progress and advantages of those systems, paving the way for a better understanding of employing 1D nanomaterials for energy storage.

A one-dimensional crystal model is constructed with a complex periodic potential. A wave function solution for the crystal model is derived without relying on Bloch functions. The new wave function solution of this model is shown to correspond to the solution for the probability amplitude of a two-level system. The energy discriminant is evaluated using an analytic formula derived from the probability amplitude solution, and based on an expansion parameter related to the energy and potential amplitude. From the wave function energy discriminant the crystal band structure is derived and related to standard energy bands and gaps. It is also shown that several of the properties of the two-level system apply to the one-dimensional crystal model. The two-level system solution which evolves in time is shown to manifest as a spatial configuration of the one-dimensional crystal model. The sensitivity of the wave function probability density is interpreted in the context of the new solution. The spatial configuration of the wave function, and the appearance of a long wavelength in the wave function probability density is explained in terms of the properties of Bessel functions

that in the weakly-coupled limit the inter-tube interaction is similar to a zero-range term with a suitable rescaled strength. This allows us to address the corresponding many-body physics of the system by constructing a model where bound chains with one molecule in each tube are the effective degrees of freedom......We consider one-dimensional tubes containing bosonic polar molecules. The long-range dipole-dipole interactions act both within a single tube and between different tubes. We consider arbitrary values of the externally aligned dipole moments with respect to the symmetry axis of the tubes. The few....... This model can be mapped onto one-dimensional Hamiltonians for which exact solutions are known....

This paper applies the variational iteration method to solve the Cauchy problem arising in onedimensional nonlinear thermoelasticity. The advantage of this method is to overcome the difficulty of calculation of Adomian's polynomials in the Adomian's decomposition method. The numerical results of this method are compared with the exact solution of an artificial model to show the efficiency of the method. The approximate solutions show that the variational iteration method is a powerful mathematical tool for solving nonlinear problems

A novel type of one-dimensional (1D) photonic crystal formed by the array of periodically located stacks of alternating graphene and dielectric stripes embedded into a background dielectric medium is proposed. The wave equation for the electromagnetic wave propagating in such structure solved in the framework of the Kronig-Penney model. The frequency band structure of 1D graphene-based photonic crystal is obtained analytically as a function of the filling factor and the thickness of the diele...

Although inextricable, the act of leading, the leader, and outcome of leadership are unique entities. Lack of such differentiation may ensnare novice leaders in broad suppositions. This conceptual article introduces a tool for analyzing leadership. Leaders can leverage the model to evaluate the act of leading, in route, via a measurable trajectory…

Automated calibration of complex deterministic water quality models with a large number of biogeochemical parameters can reduce time-consuming iterative simulations involving empirical judgements of model fit. We undertook autocalibration of the one-dimensional hydrodynamic-ecological lake model DYRESM-CAEDYM, using a Monte Carlo sampling (MCS) method, in order to test the applicability of this procedure for shallow, polymictic Lake Rotorua (New Zealand). The calibration procedure involved independently minimizing the root-mean-square error (RMSE), maximizing the Pearson correlation coefficient (r) and Nash-Sutcliffe efficient coefficient (Nr) for comparisons of model state variables against measured data. An assigned number of parameter permutations was used for 10 000 simulation iterations. The "optimal" temperature calibration produced a RMSE of 0.54 °C, Nr value of 0.99, and r value of 0.98 through the whole water column based on comparisons with 540 observed water temperatures collected between 13 July 2007 and 13 January 2009. The modeled bottom dissolved oxygen concentration (20.5 m below surface) was compared with 467 available observations. The calculated RMSE of the simulations compared with the measurements was 1.78 mg L-1, the Nr value was 0.75, and the r value was 0.87. The autocalibrated model was further tested for an independent data set by simulating bottom-water hypoxia events from 15 January 2009 to 8 June 2011 (875 days). This verification produced an accurate simulation of five hypoxic events corresponding to DO < 2 mg L-1 during summer of 2009-2011. The RMSE was 2.07 mg L-1, Nr value 0.62, and r value of 0.81, based on the available data set of 738 days. The autocalibration software of DYRESM-CAEDYM developed here is substantially less time-consuming and more efficient in parameter optimization than traditional manual calibration which has been the standard tool practiced for similar complex water quality models.

Thermal conducitivity of one-dimensional nonlinear systems typically diverges in the thermodynamic limit, whenever the momentum is conserved (i.e. in the absence of interactions with an external substrate). Evidence comes from detailed studies of Fermi-Pasta-Ulam and diatomic Toda chains. Here, we discuss the first example of a one-dimensional system obeying Fourier law : a chain of coupled rotators. Numerical estimates of the thermal conductivity obtained by simulating a chain in contact with two thermal baths at different temperatures are found to be consistent with those ones based on linear response theory. The dynamics of the Fourier modes provides direct evidence of energy diffusion. The finiteness of the conductivity is traced back to the occurrence of phase-jumps. Our conclusions are confirmed by the analysis of two variants of the rotator model.

This is a comparative study by using different onedimensional diffusion codes which are available at our Nuclear Engineering Department. Some modifications have been made in the used codes to fit the problems. One of the codes, DIFFUSE, solves the neutron diffusion equation in slab, cylindrical and spherical geometries by using 'Forward elimination- Backward substitution' technique. DIFFUSE code calculates criticality, critical dimensions and critical material concentrations and adjoint fluxes as well. It is used for the space and energy dependent neutron flux distribution. The whole scattering matrix can be used if desired. Normalisation of the relative flux distributions to the reactor power, plotting of the flux distributions and leakage terms for the other two dimensions have been added. Some modifications also have been made for the code output. Two Benchmark problems have been calculated with the modified version and the results are compared with BBD code which is available at our department and uses same techniques of calculation. Agreements are quite good in results such as k-eff and the flux distributions for the two cases studies. (author)

Highlights: ► Effect of using diesel or biodiesel on injector hydraulic behavior has been analyzed. ► Single and main + post injections have been studied for different injection pressures. ► Higher viscosity affects needle dynamics, especially for low injection pressure. ► The post injection masses are lower for biodiesel fuel despite its higher density. ► Modified injector has been proposed to compensate the differences between the fuels. - Abstract: The influence of using biodiesel fuels on the hydraulic behavior of a solenoid operated common rail injection system has been explored by means of a one-dimensionalmodel. This model has been previously obtained, including a complete characterization of the different components of the injector (mainly the nozzle, the injector holder and the electrovalve), and extensively validated by means of mass flow rate results under different conditions. After that, both single and multiple injection strategies have been analyzed, using a standard diesel fuel and rapeseed methyl ester (RME) as working fluids. Single long injections allowed the characterization of the hydraulic delay of the injector, the needle dynamics and the discharge capability of the couple injector-nozzle for the two fuels considered. Meanwhile, the effect of biodiesel on main plus post injection strategies has been evaluated in several aspects, such as the separation of the two injections or the effect of the main injection on the post injection fueling. Finally, a modification in the injector hardware has been proposed in order to have similar performances using biodiesel as the original injector configuration using standard diesel fuel.

COMPASS is an indigenously developed Chinese global navigation satellite system and will share many features in common with GPS (Global Positioning System). Since the ultra-tight GPS/INS (Inertial Navigation System) integration shows its advantage over independent GPS receivers in many scenarios, the federated ultra-tight COMPASS/INS integration has been investigated in this paper, particularly, by proposing a simplified prefilter model. Compared with a traditional prefilter model, the state space of this simplified system contains only carrier phase, carrier frequency and carrier frequency rate tracking errors. A two-quadrant arctangent discriminator output is used as a measurement. Since the code tracking error related parameters were excluded from the state space of traditional prefilter models, the code/carrier divergence would destroy the carrier tracking process, and therefore an adaptive Kalman filter algorithm tuning process noise covariance matrix based on state correction sequence was incorporated to compensate for the divergence. The federated ultra-tight COMPASS/INS integration was implemented with a hardware COMPASS intermediate frequency (IF), and INS's accelerometers and gyroscopes signal sampling system. Field and simulation test results showed almost similar tracking and navigation performances for both the traditional prefilter model and the proposed system; however, the latter largely decreased the computational load. PMID:23012564

Full Text Available COMPASS is an indigenously developed Chinese global navigation satellite system and will share many features in common with GPS (Global Positioning System. Since the ultra-tight GPS/INS (Inertial Navigation System integration shows its advantage over independent GPS receivers in many scenarios, the federated ultra-tight COMPASS/INS integration has been investigated in this paper, particularly, by proposing a simplified prefilter model. Compared with a traditional prefilter model, the state space of this simplified system contains only carrier phase, carrier frequency and carrier frequency rate tracking errors. A two-quadrant arctangent discriminator output is used as a measurement. Since the code tracking error related parameters were excluded from the state space of traditional prefilter models, the code/carrier divergence would destroy the carrier tracking process, and therefore an adaptive Kalman filter algorithm tuning process noise covariance matrix based on state correction sequence was incorporated to compensate for the divergence. The federated ultra-tight COMPASS/INS integration was implemented with a hardware COMPASS intermediate frequency (IF, and INS’s accelerometers and gyroscopes signal sampling system. Field and simulation test results showed almost similar tracking and navigation performances for both the traditional prefilter model and the proposed system; however, the latter largely decreased the computational load.

Full Text Available Automated calibration of complex deterministic water quality models with a large number of biogeochemical parameters can reduce time-consuming iterative simulations involving empirical judgements of model fit. We undertook autocalibration of the one-dimensional hydrodynamic-ecological lake model DYRESM-CAEDYM, using a Monte Carlo sampling (MCS method, in order to test the applicability of this procedure for shallow, polymictic Lake Rotorua (New Zealand. The calibration procedure involved independently minimizing the root-mean-square error (RMSE, maximizing the Pearson correlation coefficient (r and Nash–Sutcliffe efficient coefficient (Nr for comparisons of model state variables against measured data. An assigned number of parameter permutations was used for 10 000 simulation iterations. The "optimal" temperature calibration produced a RMSE of 0.54 °C, Nr value of 0.99, and r value of 0.98 through the whole water column based on comparisons with 540 observed water temperatures collected between 13 July 2007 and 13 January 2009. The modeled bottom dissolved oxygen concentration (20.5 m below surface was compared with 467 available observations. The calculated RMSE of the simulations compared with the measurements was 1.78 mg L−1, the Nr value was 0.75, and the r value was 0.87. The autocalibrated model was further tested for an independent data set by simulating bottom-water hypoxia events from 15 January 2009 to 8 June 2011 (875 days. This verification produced an accurate simulation of five hypoxic events corresponding to DO

Exchange symmetry in acceleration partitions the configuration space of an N particle one-dimensional gravitational system (OGS) into N! equivalent cells. We take advantage of the resulting small angular separation between the forces in neighboring cells to construct a related integrable version of the system that takes the form of a central force problem in N-1 dimensions. The properties of the latter, including the construction of trajectories and possible continuum limits, are developed. Dynamical simulation is employed to compare the two models. For some initial conditions, excellent agreement is observed

We study the acoustic and electronic properties of one-dimensional quasicrystals. Both numerical (nonperturbative) and analytical (perturbative) results are shown. The phonon and electronic spectra exhibit a self-similar hierarchy of gaps and many localized states in the gaps. We study quasiperiodic structures with any number of layers and several types of boundary conditions. We discuss the connection between our phonon model and recent experiments on quasiperiodic GaAs-AlAs superlattices. We predict the existence of many gap states localized at the surfaces

Full Text Available Paraphrase the title of the famous essay by Herbert Marcuse, since the image has traditionally been generated of man, masculinity, has been one-dimensional. I mean, the man was characterized by traits and behaviors established and entrenched since ancient time, considering all other distinguishing signs as mere deviations from the normative improper. But observe that this undeniable reality, as analyzed various researchers through what has come to be called Men's studies, has proven to be a fallacy difficult to maintain throughout history and today turns into fallacious and ineffective against changes in our current existing corporate models.

Full Text Available Nowadays, Marcuse’s main book One-Dimensional Man is almost obsolete, or rather passé. However, there are reasons to renew the reading of his book because of “the crisis of capitalism,” and the prevailing framework of technological domination in “advanced industrial society” in which we live today. “The new forms of control” in “advanced industrial societies” have replaced traditional methods of political and economic administration. The dominant structural element of “advanced industrial society” has become a technical and scientific apparatus of production and distribution of technology and administrative practice based on application of impersonal rules by a hierarchy of associating authorities. Technology has been liberated from the control of particular interests, and it has become the factor of domination in itself. Technological domination stems from the technical development of the productive apparatus that reproduces its ability into all spheres of social life (cultural, political, and economic. Based upon this consideration, in this paper, I will examine Marcuse’s ideas of “the new forms of control,” which creates a one–dimensional society. Marcuse’s fundamental thesis in One-Dimensional Man is that technological rationality is the most dominant factor in an “advanced industrial society,” which unites two earlier opposing forces of dissent: the bourgeoisie and the proletariat.

In this article,the qualities of Wigner function and the corresponding stationary perturbation theory are introduced and applied to one-dimensional infinite potential well and one-dimensional harmonic oscillator, and then the particular Wigner function of one-dimensional infinite potential well is specified and a special constriction effect in its pure state Wigner function is discovered, to which,simultaneously, a detailed and reasonable explanation is elaborated from the perspective of uncertainty principle. Ultimately, the amendment of Wigner function and energy of one-dimensional infinite potential well and one-dimensional harmonic oscillator under perturbation are calculated according to stationary phase space perturbation theory. (authors)

In this article we discuss charge and spin separation and quantum interference in one-dimensionalmodels. After a short introduction we briefly present the Hubbard and Luttinger models and discuss some of the known exact results. We study numerically the charge and spin separation in the Hubbard model. The time evolution of a wave packet is obtained and the charge and spin densities are evaluated for different times. The charge and spin wave packets propagate with different velocities. The results are interpreted in terms of the Bethe-ansatz solution. In section IV we study the effect of charge and spin separation on the quantum interference in a Aharonov-Bohm experiment. By calculating the one-particle propagators of the Luttinger model for a mesoscopic ring with a magnetic field we calculate the Aharonov-Bohm conductance. The conductance oscillates with the magnetic field with a characteristic frequency that depends on the charge and spin velocities. (author)

While the general theory for the terminal-initial value problem for mean-field games (MFGs) has achieved a substantial progress, the corresponding forward–forward problem is still poorly understood—even in the one-dimensional setting. Here, we consider one-dimensional forward–forward MFGs, study the existence of solutions and their long-time convergence. First, we discuss the relation between these models and systems of conservation laws. In particular, we identify new conserved quantities and study some qualitative properties of these systems. Next, we introduce a class of wave-like equations that are equivalent to forward–forward MFGs, and we derive a novel formulation as a system of conservation laws. For first-order logarithmic forward–forward MFG, we establish the existence of a global solution. Then, we consider a class of explicit solutions and show the existence of shocks. Finally, we examine parabolic forward–forward MFGs and establish the long-time convergence of the solutions.

While the general theory for the terminal-initial value problem for mean-field games (MFGs) has achieved a substantial progress, the corresponding forward–forward problem is still poorly understood—even in the one-dimensional setting. Here, we consider one-dimensional forward–forward MFGs, study the existence of solutions and their long-time convergence. First, we discuss the relation between these models and systems of conservation laws. In particular, we identify new conserved quantities and study some qualitative properties of these systems. Next, we introduce a class of wave-like equations that are equivalent to forward–forward MFGs, and we derive a novel formulation as a system of conservation laws. For first-order logarithmic forward–forward MFG, we establish the existence of a global solution. Then, we consider a class of explicit solutions and show the existence of shocks. Finally, we examine parabolic forward–forward MFGs and establish the long-time convergence of the solutions.

While the general theory for the terminal-initial value problem for mean-field games (MFGs) has achieved a substantial progress, the corresponding forward–forward problem is still poorly understood—even in the one-dimensional setting. Here, we consider one-dimensional forward–forward MFGs, study the existence of solutions and their long-time convergence. First, we discuss the relation between these models and systems of conservation laws. In particular, we identify new conserved quantities and study some qualitative properties of these systems. Next, we introduce a class of wave-like equations that are equivalent to forward–forward MFGs, and we derive a novel formulation as a system of conservation laws. For first-order logarithmic forward–forward MFG, we establish the existence of a global solution. Then, we consider a class of explicit solutions and show the existence of shocks. Finally, we examine parabolic forward–forward MFGs and establish the long-time convergence of the solutions.

In a companion paper [Phys. Rev. E 97, 043115 (2018), 10.1103/PhysRevE.97.043115], a formalism allowing to describe viscous fibers as one-dimensional objects was developed. We apply it to the special case of a viscous fluid torus. This allows to highlight the differences with the basic viscous string model and with its viscous rod model extension. In particular, an elliptic deformation of the torus section appears because of surface tension effects, and this cannot be described by viscous string nor viscous rod models. Furthermore, we study the Rayleigh-Plateau instability for periodic deformations around the perfect torus, and we show that the instability is not sufficient to lead to the torus breakup in several droplets before it collapses to a single spherical drop. Conversely, a rotating torus is dynamically attracted toward a stationary solution, around which the instability can develop freely and split the torus in multiple droplets.

A lattice-gas model with heterogeneity is developed for the description of fluid condensation in finite sized one-dimensional pores of arbitrary shape. Mapping to the random-field Ising model allows an exact solution of the model to be obtained at zero-temperature, reproducing the experimentally observed dependence of the amount of fluid adsorbed in the pore on external pressure. It is demonstrated that the disorder controls the sorption for long pores and can result in H2-type hysteresis. Finite-temperature Metropolis dynamics simulations support analytical findings in the limit of low temperatures. The proposed framework is viewed as a fundamental building block of the theory of capillary condensation necessary for reliable structural analysis of complex porous media from adsorption-desorption data.

The original Kronig-Penney lattice, which had delta function interactions at the end of each of the equal segments, seems a good model for the motion of neutrons in a linear lattice if the strength b of the δ functions depends of the energy of the neutrons, i.e., b(E). We derive the equation for the transmission bands and consider the relations of b(E) with the R(E) function discussed in a previous paper. We note the great difference in the behavior of the bands when b(E) is constant and when it is related with a single resonance of the R function. (Author)

We compare results of different numerical approaches to compute the NMR relaxation rate 1 /T1 in quasi one-dimensional (1d) antiferromagnets. In the purely 1d regime, recent numerical simulations using DMRG have provided the full crossover behavior from classical regime at high temperature to universal Tomonaga-Luttinger liquid at low-energy (in the gapless case) or activated behavior (in the gapped case). For quasi 1d models, we can use mean-field approaches to reduce the problem to a 1d one that can be studied using DMRG. But in some cases, we can also simulate the full microscopic model using quantum Monte-Carlo techniques. This allows to compute dynamical correlations in imaginary time and we will discuss recent advances to perform stochastic analytic continuation to get real frequency spectra. Finally, we connect our results to experiments on various quasi 1d materials.

We present results of transport measurements of individual suspended electrospun nanofibers Poly(methyl methacrylate)-multiwalled carbon nanotubes. The nanofiber is comprised of highly aligned consecutive multiwalled carbon nanotubes. We have confirmed that at the range temperature from room temperature down to ∼60 K, the conductance behaves as power-law of temperature with an exponent of α ∼ 2.9−10.2. The current also behaves as power law of voltage with an exponent of β ∼ 2.3−8.6. The power-law behavior is a footprint for onedimensional transport. The possible models of this confined system are discussed. Using the model of Luttinger liquid states in series, we calculated the exponent for tunneling into the bulk of a single multiwalled carbon nanotube α{sub bulk} ∼ 0.06 which agrees with theoretical predictions.

We build a general formalism to describe thin viscous jets as one-dimensional objects with an internal structure. We present in full generality the steps needed to describe the viscous jets around their central line, and we argue that the Taylor expansion of all fields around that line is conveniently expressed in terms of symmetric trace-free tensors living in the two dimensions of the fiber sections. We recover the standard results of axisymmetric jets and we report the first and second corrections to the lowest order description, also allowing for a rotational component around the axis of symmetry. When applied to generally curved fibers, the lowest order description corresponds to a viscous string model whose sections are circular. However, when including the first corrections, we find that curved jets generically develop elliptic sections. Several subtle effects imply that the first corrections cannot be described by a rod model since it amounts to selectively discard some corrections. However, in a fast rotating frame, we find that the dominant effects induced by inertial and Coriolis forces should be correctly described by rod models. For completeness, we also recover the constitutive relations for forces and torques in rod models and exhibit a missing term in the lowest order expression of viscous torque. Given that our method is based on tensors, the complexity of all computations has been beaten down by using an appropriate tensor algebra package such as xAct, allowing us to obtain a one-dimensional description of curved viscous jets with all the first order corrections consistently included. Finally, we find a description for straight fibers with elliptic sections as a special case of these results, and recover that ellipticity is dynamically damped by surface tension. An application to toroidal viscous fibers is presented in the companion paper [Pitrou, Phys. Rev. E 97, 043116 (2018), 10.1103/PhysRevE.97.043116].

This study analyzes thermoelectric properties of a one-dimensional random conductor which shows localization effects and simultaneously includes resonant scatterers yielding sharp conductance resonances. These sharp features give rise to a distinct behavior of the Seebeck coefficient in finite systems and incorporate the degree of localization as a means to enhance thermoelectric performance, in principle. The model for noninteracting electrons is discussed within the Landauer-Büttiker formalism such that analytical treatment is possible for a wide range of properties, if a special averaging scheme is applied. The approximations in the averaging procedure are tested with numerical evaluations showing good qualitative agreement, with some limited quantitative disagreement. The validity of low-temperature Mott's formula is determined and a good approximation is developed for the intermediate temperature range. In both regimes the intricate interplay between Anderson localization due to disorder and conductance resonances of the disorder potential is analyzed.

We have studied the magnon transmission through of one-dimensional magnonic k-component Fibonacci structures, where k different materials are arranged in accordance with the following substitution rule: S{sub n}{sup (k)}=S{sub n−1}{sup (k)}S{sub n−k}{sup (k)} (n≥k=0,1,2,…), where S{sub n}{sup (k)} is the nth stage of the sequence. The calculations were carried out in exchange dominated regime within the framework of the Heisenberg model and taking into account the RPA approximation. We have considered multilayers composed of simple cubic spin-S Heisenberg ferromagnets, and, by using the powerful transfer-matrix method, the spin wave transmission is obtained. It is demonstrated that the transmission coefficient has a rich and interesting magnonic pass- and stop-bands structures, which depends on the frequency of magnons and the k values.

Incidents caused by fire and combat operations can heat energetic materials that may lead to thermal explosion and result in structural damage and casualty. Some explosives may thermally explode at fairly low temperatures (< 100 C) and the violence from thermal explosion may cause a significant damage. Thus it is important to understand the response of energetic materials to thermal insults. The OneDimensional Time to Explosion (ODTX) system at the Lawrence Livermore National Laboratory has been used for decades to measure times to explosion, threshold thermal explosion temperature, and determine kinetic parameters of energetic materials. Samples of different configurations (pressed part, powder, paste, and liquid) can be tested in the system. The ODTX testing can also provide useful data for assessing the thermal explosion violence of energetic materials. This report summarizes the recent ODTX experimental data and modeling results for 2,6-diamino-3,5-dintropyrazine (ANPZ).

The authors study a class of periodic Schrödinger operators, which in distinguished cases can be proved to have linear band-crossings or "Dirac points". They then show that the introduction of an "edge", via adiabatic modulation of these periodic potentials by a domain wall, results in the bifurcation of spatially localized "edge states". These bound states are associated with the topologically protected zero-energy mode of an asymptotic one-dimensional Dirac operator. The authors' model captures many aspects of the phenomenon of topologically protected edge states for two-dimensional bulk structures such as the honeycomb structure of graphene. The states the authors construct can be realized as highly robust TM-electromagnetic modes for a class of photonic waveguides with a phase-defect.

We present a study of one-dimensional magnetic system using field theory methods. We studied the discreteness effects in a classical anisotropic onedimensional antiferromagnet in an external magnetic field. It is shown that for TMMC, at the temperatures and magnetic fields where most experiments have been done, the corrections are small and can be neglected. (author)

Understanding the non-equilibrium dynamics of isolated quantum many-body systems is an open problem on vastly different energy, length, and time scales. Examples range from the dynamics of the early universe and heavy-ion collisions to the subtle coherence and transport properties in condensed matter physics. However, realizations of such quantum many-body systems, which are both well isolated from the environment and accessible to experimental study are scarce. This thesis presents a series of experiments with ultracold one-dimensional Bose gases. These gases combine a nearly perfect isolation from the environment with many well-established methods to manipulate and probe their quantum states. This makes them an ideal model system to explore the physics of quantum many body systems out of equilibrium. In the experiments, a well-defined non-equilibrium state is created by splitting a single one-dimensional gas coherently into two parts. The relaxation of this state is probed using matter-wave interferometry. The Observations reveal the emergence of a prethermalized steady state which differs strongly from thermal equilibrium. Such thermal-like states had previously been predicted for a large variety of systems, but never been observed directly. Studying the relaxation process in further detail shows that the thermal correlations of the prethermalized state emerge locally in their final form and propagate through the system in a light-cone-like evolution. This provides first experimental evidence for the local relaxation conjecture, which links relaxation processes in quantum many-body systems to the propagation of correlations. Furthermore, engineering the initial state of the evolution demonstrates that the prethermalized state is described by a generalized Gibbs ensemble, an observation which substantiates the importance of this ensemble as an extension of standard statistical mechanics. Finally, an experiment is presented, where pairs of gases with an atom

A phenomenological picture of ac hopping in the symmetric hopping model (regular lattice, equal site energies, random energy barriers) is proposed according to which conduction in the extreme disorder limit is dominated by essentially one-dimensional "percolation paths." Modeling a percolation path...... as strictly onedimensional with a sharp jump rate cutoff leads to an expression for the universal ac conductivity that fits computer simulations in two and three dimensions better than the effective medium approximation....

Spin glasses are paradigmatic models that deliver concepts relevant for a variety of systems. However, rigorous analytical results are difficult to obtain for spin-glass models, in particular for realistic short-range models. Therefore large-scale numerical simulations are the tool of choice. Concepts and algorithms derived from the study of spin glasses have been applied to diverse fields in computer science and physics. In this work a one-dimensional long-range spin-glass model with power-law interactions is discussed. The model has the advantage over conventional systems in that by tuning the power-law exponent of the interactions the effective space dimension can be changed thus effectively allowing the study of large high-dimensional spin-glass systems to address questions as diverse as the existence of an Almeida-Thouless line, ultrametricity and chaos in short range spin glasses. Furthermore, because the range of interactions can be changed, the model is a formidable test-bed for optimization algorithms

This course is part of the ENSTA 3rd year thermal hydraulics program (nuclear power option). Its purpose is to provide the theoretical basis and main physical notions pertaining to two-phase flow, mainly focussed on water-steam flows. The introduction describes the physical specificities of these flows, emphasizing their complexity. The mathematical bases are then presented (partial derivative equations), leading to a one-dimensional type, simplified description. Balances drawn up for a pipe length volume are used to introduce the mass conservation. motion and energy equations for each phase. Various postulates used to simplify two-phase models are presented, culminating in homogeneous model definitions and equations, several common examples of which are given. The model is then applied to the calculation of pressure drops in two-phase flows. This involves presenting the models most frequently used to represent pressure drops by friction or due to pipe irregularities, without giving details (numerical values of parameters). This chapter terminates with a brief description of static and dynamic instabilities in two-phase flows. Finally, heat transfer conditions frequently encountered in liquid-steam flows are described, still in the context of a 1D model. This chapter notably includes reference to under-saturated boiling conditions and the various forms of DNB. The empirical heat transfer laws are not discussed in detail. Additional material is appended, some of which is in the form of corrected exercises. (author). 6 appends

Both point- and finite-source stochastic one-dimensional ground motion models, coupled to vertically propagating equivalent-linear shear-wave site response models are validated using an extensive set of strong motion data as part of the Yucca Mountain Project. The validation and comparison exercises are presented entirely in terms of 5% damped pseudo absolute response spectra. The study consists of a quantitative analyses involving modeling nineteen well-recorded earthquakes, M 5.6 to 7.4 at over 600 sites. The sites range in distance from about 1 to about 200 km in the western US (460 km for central-eastern US). In general, this validation demonstrates that the stochastic point- and finite-source models produce accurate predictions of strong ground motions over the range of 0 to 100 km and for magnitudes M 5.0 to 7.4. The stochastic finite-source model appears to be broadband, producing near zero bias from about 0.3 Hz (low frequency limit of the analyses) to the high frequency limit of the data (100 and 25 Hz for response and Fourier amplitude spectra, respectively)

Quantum walks are interesting models with potential applications to quantum algorithms and physical processes such as photosynthesis. In this paper, we study two models of one-dimensional quantum walks, namely, quantum walks with a moving absorbing wall and quantum walks with one stationary and one moving absorbing wall. For the former, we calculate numerically the survival probability, the rate of change of average position, and the rate of change of standard deviation of the particle's position in the long time limit for different wall velocities. Moreover, we also study the asymptotic behavior and the dependence of the survival probability on the initial particle's state. While for the latter, we compute the absorption probability of the right stationary wall for different velocities and initial positions of the left wall boundary. The results for these two models are compared with those obtained for the classical model. The difference between the results obtained for the quantum and classical models can be attributed to the difference in the probability distributions.

Both point- and finite-source stochastic one-dimensional ground motion models, coupled to vertically propagating equivalent-linear shear-wave site response models are validated using an extensive set of strong motion data as part of the Yucca Mountain Project. The validation and comparison exercises are presented entirely in terms of 5% damped pseudo absolute response spectra. The study consists of a quantitative analyses involving modeling nineteen well-recorded earthquakes, M 5.6 to 7.4 at over 600 sites. The sites range in distance from about 1 to about 200 km in the western US (460 km for central-eastern US). In general, this validation demonstrates that the stochastic point- and finite-source models produce accurate predictions of strong ground motions over the range of 0 to 100 km and for magnitudes M 5.0 to 7.4. The stochastic finite-source model appears to be broadband, producing near zero bias from about 0.3 Hz (low frequency limit of the analyses) to the high frequency limit of the data (100 and 25 Hz for response and Fourier amplitude spectra, respectively).

The plasma properties of the quasi-one-dimensional ring in the threshold cases of low and high frequencies, corresponding to the plasma oscillations and dielectric relaxation are studied within the frames of the classical approach. The plasma oscillations spectrum and the electron dielectric relaxation frequency in the quasi-one-dimensional ring are calculated. The plasmons spectrum equidistance is identified. It is shown , that in contrast to the three-dimensional case there takes place the dielectric relaxation dispersion, wherefrom there follows the possibility of studying the carriers distribution in the quasi-one-dimensional rings through the method of the dielectric relaxation spectroscopy

The effects of the Dzyaloshinskii—Moriya (DM) and the Kaplan—Shekhtman—Entinwohlman—Aharony (KSEA) superexchange interactions on the ground state properties of the one-dimensional spin-Peilers system in open chain are studied by using the Lanczos numerical method. The study concentrates mainly on the influence of systemic dimerisation in open chain. The results show that systemic ground state energy density varies with dimerisation parameter δ in different DM interactions, and there exists a special point δ c where the DM interaction has no influence on the systemic dimerisation, no matter whether the DM interaction is relative or irrelative to systemic dimerisation (η = 1 or η = 0). The KSEA interaction has no fixed special point, but the points of intersection are dense relatively in a certain numberical value range, and sparse in other numberical value ranges. So we can conclude that the antisymmetric anisotropy DM interaction differs from the symmetric anisotropy KSEA interaction, but they are analogous in the sense of the influence of systemic dimerisation in open chain

In this thesis, we consider stationary one-dimensional mean-field games (MFGs) with or without congestion. Our aim is to understand the qualitative features of these games through the analysis of explicit solutions. We are particularly interested

voltage characteristics of a one-dimensional molecular wire with odd number of ... lem, although interesting both from a fundamental point of view and in terms of ..... SKP acknowledges the DST, Government of India, for financial support.

Models of superconductors having a quasi-one-dimensional crystal structure based on the convoluted into a tube Ginzburg sandwich, which comprises a layered dielectric–metal–dielectric structure, have been suggested. The critical crystal chemistry parameters of the Ginzburg sandwich determining the possibility of the emergence of superconductivity and the T c value in layered high-T c cuprates, which could have the same functions in quasi-one-dimensional fragments (sandwich-type tubes), have been examined. The crystal structures of known low-temperature superconductors, in which one can mark out similar quasi-one-dimensional fragments, have been analyzed. Five compounds with quasi-one-dimensional structures, which can be considered as potential parents of new superconductor families, possibly with high transition temperatures, have been suggested. The methods of doping and modification of these compounds are provided. (paper)

This paper is concerned with the periodic solutions for the onedimensional nonlinear wave equation with either constant or variable coefficients. The constant coefficient model corresponds to the classical wave equation, while the variable coefficient model arises from the forced vibrations of a nonhomogeneous string and the propagation of seismic waves in nonisotropic media. For finding the periodic solutions of variable coefficient wave equation, it is usually required that the coefficient u (x) satisfies ess infηu (x) > 0 with ηu (x) = 1/2 u″/u - 1/4 (u‧/u)2, which actually excludes the classical constant coefficient model. For the case ηu (x) = 0, it is indicated to remain an open problem by Barbu and Pavel (1997) [6]. In this work, for the periods having the form T = 2p-1/q (p , q are positive integers) and some types of boundary value conditions, we find some fundamental properties for the wave operator with either constant or variable coefficients. Based on these properties, we obtain the existence of periodic solutions when the nonlinearity is monotone and bounded. Such nonlinearity may cross multiple eigenvalues of the corresponding wave operator. In particular, we do not require the condition ess infηu (x) > 0.

Observations of galaxies over large distances reveal the possibility of a fractal distribution of their positions. The source of fractal behavior is the lack of a length scale in the two body gravitational interaction. However, even with new, larger, sample sizes from recent surveys, it is difficult to extract information concerning fractal properties with confidence. Similarly, three-dimensional N-body simulations with a billion particles only provide a thousand particles per dimension, far too small for accurate conclusions. With one-dimensionalmodels these limitations can be overcome by carrying out simulations with on the order of a quarter of a million particles without compromising the computation of the gravitational force. Here the multifractal properties of two of these models that incorporate different features of the dynamical equations governing the evolution of a matter dominated universe are compared. For each model at least two scaling regions are identified. By employing criteria from dynamical systems theory it is shown that only one of them can be geometrically significant. The results share important similarities with galaxy observations, such as hierarchical clustering and apparent bifractal geometry. They also provide insights concerning possible constraints on length and time scales for fractal structure. They clearly demonstrate that fractal geometry evolves in the mu (position, velocity) space. The observed patterns are simply a shadow (projection) of higher-dimensional structure.

of its star tracker, the Advanced Stellar Compass (ASC). The ASC features low cost, low mass, low power, low magnetic disturbance, autonomous operation, a high level of functionality and the high precision. These features are enabled by the use of advanced optical and electronic design which permit...

In this paper we study bond percolation on a one-dimensional chain with power-law bond probability C /rd +σ , where r is the distance length between distinct sites and d =1 . We introduce and test an order-N Monte Carlo algorithm and we determine as a function of σ the critical value Cc at which percolation occurs. The critical exponents in the range 0 values for Cc are compared with a known exact bound, while the critical exponent ν is compared with results from mean-field theory, from an expansion around the point σ =1 and from the ɛ -expansion used with the introduction of a suitably defined effective dimension deff relating the long-range model with a short-range one in dimension deff. We finally present a formulation of our algorithm for bond percolation on general graphs, with order N efficiency on a large class of graphs including short-range percolation and translationally invariant long-range models in any spatial dimension d with σ >0 .

by the complexity of the field. Next is a section on how to move from an idea or problem to a research question by placing a concrete idea or problem within a conceptual, theoretical framework. The following sections are structured around an overview model of approaches to medical education research, 'The research...... compass'. Core to the model is the conceptual, theoretical framework that is the key to any direction. The compass depicts four main categories of research approaches that can be applied when studying medical education phenomena, 'Explorative studies'; 'Experimental studies'; 'Observational studies......This AMEE Guide offers an introduction to research in medical education. It is intended for those who are contemplating conducting research in medical education but are new to the field. The Guide is structured around the process of transforming ideas and problems into researchable questions...

A one-dimensional photon transport code BERMUDA-1DG has been developed for spherical and infinite slab geometries. The purpose of development is to equip the function of gamma rays calculation for the BERMUDA code system, which was developed by 1983 only for neutron transport calculation as a preliminary version. A group constants library has been prepared for 30 nuclides, and it now consists of the 36-group total cross sections and secondary gamma ray yields by the 120-group neutron flux. For the Compton scattering, group-angle transfer matrices are accurately obtained by integrating the Klein-Nishina formula taking into account the energy and scattering angle correlation. The pair production cross sections are now calculated in the code from atomic number and midenergy of each group. To obtain angular flux distribution, the transport equation is solved in the same way as in case of neutron, using the direct integration method in a multigroup model. Both of an independent gamma ray source problem and a neutron-gamma source problem are possible to be solved. This report is written as a user's manual with a brief description of the calculational method. (author)

The Poisson–Boltzmann (PB) equation is conventionally used to model the equilibrium of bulk ionic species in different media and solvents. In this paper we study a new Poisson–Boltzmann type (PB n ) equation with a small dielectric parameter ε 2 and non-local nonlinearity which takes into consideration the preservation of the total amount of each individual ion. This equation can be derived from the original Poisson–Nernst–Planck system. Under Robin-type boundary conditions with various coefficient scales, we demonstrate the asymptotic behaviours of one-dimensional solutions of PB n equations as the parameter ε approaches zero. In particular, we show that in case of electroneutrality, i.e. α = β, solutions of 1D PB n equations have a similar asymptotic behaviour as those of 1D PB equations. However, as α ≠ β (non-electroneutrality), solutions of 1D PB n equations may have blow-up behaviour which cannot be found in 1D PB equations. Such a difference between 1D PB and PB n equations can also be verified by numerical simulations

A Kronig-Penney model with a constant electric filed for a non-interacting electron is used to study the transmission properties of Anderson transition in one-dimensional (1-D) systems with disordered strengths of δ-function potentials. we examined the cases where the potential varies uniformly from O to W (barriers) or from -W to O (wells) for a given disorder W. Mainly, we observe unexpected abrupt transition at the points E + Fx = n 2 π 2 . However, these transitions are related to the small oscillations observed by Soukoulis et al. in the mixed case (wells and barriers). An interesting feature in the wells is that in the presence of a small field the states become more localized and the localization length decrease up to a minimum for a critical value F m . In the end, we have studied the effect of the disorder on the Anderson transition by the mean of the participation ratio and the localization length. (author). 27 refs, 6 figs

A recent theory based on fluctuating hydrodynamics predicts that one-dimensional interacting systems with particle, momentum, and energy conservation exhibit anomalous transport that falls into two main universality classes. The classification is based on behavior of equilibrium dynamical correlations of the conserved quantities. One class is characterized by sound modes with Kardar-Parisi-Zhang scaling, while the second class has diffusive sound modes. The heat mode follows Lévy statistics, with different exponents for the two classes. Here we consider heat current fluctuations in two specific systems, which are expected to be in the above two universality classes, namely, a hard particle gas with Hamiltonian dynamics and a harmonic chain with momentum conserving stochastic dynamics. Numerical simulations show completely different system-size dependence of current cumulants in these two systems. We explain this numerical observation using a phenomenological model of Lévy walkers with inputs from fluctuating hydrodynamics. This consistently explains the system-size dependence of heat current fluctuations. For the latter system, we derive the cumulant-generating function from a more microscopic theory, which also gives the same system-size dependence of cumulants.

A code MARG1D has been developed which computes outer region matching data of the onedimensional Newcomb equation. Matching data play an important role in the resistive (and non ideal) Magneto-hydrodynamic (MHD) stability analysis in a tokamak plasma. The MARG1D code computes matching data by using the boundary value method or by the eigenvalue method. Variational principles are derived for the problems to be solved and a finite element method is applied. Except for the case of marginal stability, the eigenvalue method is equivalent to the boundary value method. However, the eigenvalue method has the several advantages: it is a new method of ideal MHD stability analysis for which the marginally stable state can be identified, and it guarantees numerical stability in computing matching data close to marginal stability. We perform detailed numerical experiments for a model equation with analytical solutions and for the Newcomb equation in the m=1 mode theory. Numerical experiments show that MARG1D code gives the matching data with numerical stability and high accuracy. (author)

The TRAC-BWR code development program at the Idaho National Engineering Laboratory is developing a version of the TRAC code for the U.S. Nuclear Regulatory Commission (USNRC) to provide a best-estimate analysis capability for the simulation of postulated accidents in boiling water reactor (BWR) power systems and related experimental facilities. Recent development efforts in the TRAC-BWR program have focused on improving the computational efficiency through the incorporation of a hybrid Courant- limit-violating numerical solution scheme in the one-dimensional component models and on improving code accuracy through the development of a one-dimensional neutron kinetics model. Many other improvements have been incorporated into TRAC-BWR to improve code portability, accuracy, efficiency, and maintainability. This paper will describe the one- dimensional neutron kinetics model, the generation of the required input data for this model, and present results of the first calculations using the model

Full Text Available In recent years, one-dimensional piezoelectric nanomaterials have become a research topic of interest because of their special morphology and excellent piezoelectric properties. This article presents a short review on onedimensional perovskite piezoelectric materials in different systems including Pb(Zr,TiO3, BaTiO3 and (K,NaNbO3 (KNN. We emphasize KNN as a promising lead-free piezoelectric compound with a high Curie temperature and high piezoelectric properties and describe its synthesis and characterization. In particular, details are presented for nanoscale piezoelectricity characterization of a single KNN nanocrystal by piezoresponse force microscopy. Finally, this review describes recent progress in applications based on onedimensional piezoelectric nanostructures with a focus on energy harvesting composite materials.

According to the Poincare-Bendixson theorem, a minimum of three autonomous equations is required to exhibit deterministic chaos. Because a one-dimensional quantum system is described by only two autonomous equations using de Broglie-Bohm's trajectory interpretation, chaos in one-dimensional quantum systems has long been considered impossible. We will prove in this paper that chaos phenomenon does exist in one-dimensional quantum systems, if the domain of quantum motions is extended to complex space by noting that the quantum world is actually characterized by a four-dimensional complex spacetime according to the E (∞) theory. Furthermore, we point out that the interaction between the real and imaginary parts of complex trajectories produces a new chaos phenomenon unique to quantum systems, called strong chaos, which describes the situation that quantum trajectories may emerge and diverge spontaneously without any perturbation in the initial position

In this work we carry out a theoretical analysis of the spectra of magnons in quasiperiodic magnonic crystals arranged in accordance with generalized Fibonacci sequences in the exchange regime, by using a model based on a transfer-matrix method together random-phase approximation (RPA). The generalized Fibonacci sequences are characterized by an irrational parameter {sigma}(p,q), which rules the physical properties of the system. We discussed the magnonic fractal spectra for first three generalizations, i.e., silver, bronze and nickel mean. By varying the generation number, we have found that the fragmentation process of allowed bands makes possible the emergence of new allowed magnonic bulk bands in spectra regions that were magnonic band gaps before, such as which occurs in doped semiconductor devices. This interesting property arises in one-dimensional magnonic quasicrystals fabricated in accordance to quasiperiodic sequences, without the need to introduce some deferent atomic layer or defect in the system. We also make a qualitative and quantitative investigations on these magnonic spectra by analyzing the distribution and magnitude of allowed bulk bands in function of the generalized Fibonacci number F{sub n} and as well as how they scale as a function of the number of generations of the sequences, respectively. - Highlights: Black-Right-Pointing-Pointer Quasiperiodic magnonic crystals are arranged in accordance with the generalized Fibonacci sequence. Black-Right-Pointing-Pointer Heisenberg model in exchange regime is applied. Black-Right-Pointing-Pointer We use a theoretical model based on a transfer-matrix method together random-phase approximation. Black-Right-Pointing-Pointer Fractal spectra are characterized. Black-Right-Pointing-Pointer We analyze the distribution of allowed bulk bands in function of the generalized Fibonacci number.

We show theoretically that the absorption of one-dimensional metallic-dielectric photonic crystals can be enhanced considerably over the corresponding constituent metal. By properly choosing the structural and material parameters, the absorption of one-dimensional metallic-dielectric photonic crystals can be enhanced by one order of magnitude in the visible and in the near infrared regions. It is found that the absorptance of such photonic crystals increases with increasing number of periods. Rules on how to obtain a absorption enhancement in a certain frequency range are discussed. (letter to the editor)

The purpose of this research is to examine the relationship of self-compassion and internet addiction. Participants were 261 university students who completed a questionnaire package that included the Self-compassion Scale and the Online Cognition Scale. The hypothesis model was tested through structural equation modeling. In correlation analysis,…

In Part I the fundamental form of the hydrodynamic basic equations for a one-dimensional two-phase flow (two-fluid model) is described. Discussions are concentrated on the treatment of phase change inertial force terms in the equations of motion and the author's equations of motion which have a remarkable uniqueness on the following three points. (1) To express force balance of unit mass two-phase fluid instead of that of unit volume two-phase fluid. (2) To pick up the unit existing mass and the unit flowing mass as the unit mass of two-phase fluid. (3) To apply the kinetic energy principle instead of the momentum low in the evaluation of steady inertial force term. In these three, the item (1) is for excluding a part of momentum change or kinetic energy change due to mass change of the examined part of fluid, which is independent of force. The item (2) is not to introduce a phenomenological physical model into the evaluation of phase change inertial force term. And the item (3) is for correctly applying the momentum law taking into account the difference of representative velocities between the main flow fluid (vapor phase or liquid phase) and the phase change part of fluid. In Part II, characteristics of various kinds of high speed two-phase flow are clarified theoretically by the basic equations derived. It is demonstrated that the steam-water two-phase critical flow with violent flashing and the airwater two-phase critical flow without phase change can be described with fundamentally the same basic equations. Furthermore, by comparing the experimental data from the two-phase critical discharge test and the theoretical prediction, the two-phase discharge coefficient, C D , for large sharp-edged orifice is determined as the value which is not affected by the experimental facility characteristics, etc. (author)

We develop a rigorous computational method for estimating the Lyapunov exponents in uniformly expanding regions of the phase space for one-dimensional maps. Our method uses rigorous numerics and graph algorithms to provide results that are mathematically meaningful and can be achieved in an efficient way

The current-voltage (I-V) characteristics of quasi-one-dimensional superconductors were discussed. The I-V characteristics exhibited an unusual S behavior. The dynamics of superconducting condensate and the existence of two different critical currents resulted in such an unusual behavior....

Transport properties of onedimensional Brownian diffusion under the influence of a quenched random force, distributed as a two-level Poisson process is discussed. Large time scaling laws of the position of the Brownian particle, and the probability distribution of the stationary flux going through a sample between two prescribed concentrations are studied. (author) 14 refs.; 3 figs

The development of a one-dimensional position readout system with microchannel plates, is described, for heavy ion detectors for use in a particle time-of-flight telescope and as a position sensitive device in front of an ionisation counter at the Nuclear Structure Facility. (U.K.)

The Ward identity is derived for non-relativistic fermions with two-body spin-independent interaction. Using this identity for the one-dimensional Fermi gas with backward scattering the equations of the perturbation theory are solved for the effective interaction and the collective excitations of the particle density fluctuations are obtained. (author)

The effects of the variation of atomic spacing ratio of a onedimensional quasicrystal material are investigated. The work involves the use of the solid state simulation code, Laue written by Silsbee and Drager. We are able to observe the general features of the diffraction pattern by a quasicrystal. In addition, it has been found ...

We review and analyze the method [U. Leonhardt, M.G. Raymer: Phys. Rev. Lett. 76, 1985 (1996)] for quantum-state reconstruction of one-dimensional non-relativistic wave packets from position observations. We illuminate the theoretical background of the technique and show how to extend the procedure to the continuous part of the spectrum.

The Lagrangian and the Generalized Linear Momentum are given in terms of a constant of motion for a one-dimensional autonomous system. The possibility of having an explicit Hamiltonian expression is also analyzed. The approach is applied to some dissipative systems. Copyright copyright 1996 Academic Press, Inc

In this thesis we study quantum transport in several one-dimensional systems with strong electronic interactions. The first chapter contains an introduction to the concepts treated throughout this thesis, such as the Aharonov-Bohm effect, the Kondo effect, the Fano effect and quantum state transfer.

We calculate the thermal correlation functions of the one-dimensional damped harmonic oscillator in contact with a reservoir, in an exact form by applying Green's function method. In this way the thermal fluctuations are incorporated in the Caldirola-Kanai Hamiltonian

We provide firm convincing evidence that the energy transport in a one-dimensional gas of elastically colliding free particles of unequal masses is anomalous, i.e., the Fourier law does not hold. Our conclusions are confirmed by a theoretical and numerical analysis based on a Green-Kubo-type approach specialized to momentum-conserving lattices.

Conditions for the existence of band gaps in a one-dimensional disordered array of δ-function potentials possessing short range order are developed in a relativistic framework. Both Lorentz vector and scalar type potentials are treated. The relationship between the energy gaps and the transmission properties of the array are also discussed. 20 refs., 2 figs

We have calculated the electromagnetic Brillouin precursor that arises in a one-dimensional photonic crystal that consists of two homogeneous slabs which each have a single electron resonance. This forerunner is compared with the Brillouin precursor that arises in a homogeneous double-electron

The quantisation of one-dimensional MIT bags by expanding the fields as a sum of classical modes and truncating the series after the first term is discussed. The lowest states of a bag in a world containing two scalar quark fields are obtained. Problems associated with the zero-point oscillations of the field are discussed. (Auth.)

The one-dimensional finite well is a textbook problem. We propose approximate approaches to obtain the energy levels of the well. The finite well is also encountered in semiconductor heterostructures where the carrier mass inside the well (m[subscript i]) is taken to be distinct from mass outside (m[subscript o]). A relevant parameter is the mass…

A family of one-dimensional inverse problems are considered with the goal of reconstructing velocity profiles to reasonably high accuracy. The travel-time variable change is used together with an iteration scheme to produce an effective algorithm for computation. Under modest assumptions the scheme is shown to be convergent

@@ By using the measure of von Neumann entropy, we numerically investigate quantum entanglement of an electronmoving in the one-dimensional Harper model and in the one-dimensional slowly varying potential model. Thedelocalized and localized eigenstates can be distinguished by von Neumann entropy of the individual eigenstates.There are drastic decreases in yon Neumann entropy of the individual eigenstates at mobility edges. In the curveof the spectrum averaged yon Neumann entropy as a function of potential parameter λ, a sharp transition existsat the metal-insulator transition point λc = 2. It is found that the yon Neumann entropy is a good quantity toreflect localization and metal-insulator transition.

A one-dimensional system of interacting fermions in an external potential is studied. The problem was for this purpose transformed to two classical models of statistical mechanics in two dimensions in which occasionally results were found in complementary ranges of the interaction constants of the fermion system. The conductivity appeared as a simple correlation function in both classical models. It was shown that the interaction in a one-dimensional polluted fermion system can cause an isolator-metal transition. (orig./HSI) [de

In this paper, we present a model to describe hopping transport and electrical conductivity of one-dimensional systems with off-diagonal disorder, in which electrons are transported via hopping between localized states. We find that off-diagonal disorder leads to delocalization and drastically enhances the electrical conductivity of systems. The model also quantitatively explains the temperature and electrical field dependence of the conductivity in one-dimensional systems with off-diagonal disorder. In addition, we also show the dependence of the conductivity on the strength of off-diagonal disorder

We present a new model of quantum phase transitions in matrix product systems of one-dimensional spin-1 chains and study the phases coexistence phenomenon. We find that in the thermodynamic limit the proposed system has three different quantum phases and by adjusting the control parameters we are able to realize any phase, any two phases equal coexistence and the three phases equal coexistence. At every critical point the physical quantities including the entanglement are not discontinuous and the matrix product system has long-range correlation and N-spin maximal entanglement. We believe that our work is helpful for having a comprehensive understanding of quantum phase transitions in matrix product states of one-dimensional spin chains and of certain directive significance to the preparation and control of one-dimensional spin lattice models with stable coherence and N-spin maximal entanglement. (author)

Various chemically reacting flow problems highlighting chemical and physical fundamentals rather than flow geometry are presently investigated by means of a comprehensive mathematical model that incorporates multicomponent molecular diffusion, complex chemistry, and heterogeneous processes, in the interest of obtaining sensitivity-related information. The sensitivity equations were decoupled from those of the model, and then integrated one time-step behind the integration of the model equations, and analytical Jacobian matrices were applied to improve the accuracy of sensitivity coefficients that are calculated together with model solutions.

This textbook presents finite element methods using exclusively one-dimensional elements. The aim is to present the complex methodology in an easily understandable but mathematically correct fashion. The approach of one-dimensional elements enables the reader to focus on the understanding of the principles of basic and advanced mechanical problems. The reader easily understands the assumptions and limitations of mechanical modeling as well as the underlying physics without struggling with complex mathematics. But although the description is easy it remains scientifically correct. The approach using only one-dimensional elements covers not only standard problems but allows also for advanced topics like plasticity or the mechanics of composite materials. Many examples illustrate the concepts and problems at the end of every chapter help to familiarize with the topics.

We present the results of a model-system study of the competition between superconductivity and disorder in narrow superconducting wires. As one moves from the two-dimensional regime toward the one-dimensional limit, large and systematic reductions in the superconducting transition temperature are obtained. The observed behavior extrapolates to the total destruction of superconductivity in the disordered one-dimensional limit. Our findings are in clear disagreement with a recent theoretical treatment. In addition, the superconducting fluctuations appear to be modified by disorder for the narrowest samples

We investigate many characteristic features of revival and fractional revival phenomena via derived analytic expressions for an autocorrelation function of a one-dimensional Rydberg atom with weighting probabilities modelled by a Gaussian or a Lorentzian distribution. The fractional revival phenomenon in the ionization probabilities of a one-dimensional Rydberg atom irradiated by two short half-cycle pulses is also studied. When many states are involved in the formation of the wave packet, the revival is lower and broader than the initial wave packet and the fractional revivals overlap and disappear with time

Full Text Available Intermittent filamentary dynamics of the vortex matter in superconductors is found in films of YBa_{2}Cu_{3}O_{7-δ} deposited on tilted substrates. Deposition of this material on such substrates creates parallel channels of easy flux penetration when a magnetic field is applied perpendicular to the film. As the applied field is gradually increased, magneto-optical imaging reveals that flux penetrates via numerous quasi-one-dimensional jumps. The distribution of flux avalanche sizes follows a power law, and data collapse is obtained by finite-size scaling, with the depth of the flux front used as crossover length. The intermittent behavior shows no threshold value in the applied field, in contrast to conventional flux jumping. The results strongly suggest that the quasi-one-dimensional flux jumps are of a different nature than the thermomagnetic dendritic (branching avalanches that are commonly found in superconducting films.

In this paper we report on a versatile hydrothermal approach developed to fabricate one-dimensional (1D) composite structures. Sulfur and selenium formed liquid and adsorbed onto microrods as droplets and subsequently reacted with metallic ion in solution to produce nanoparticles-decorated composite microrods. 1D composites including ZnO/CdS, ZnO/MnS, ZnO/CuS, ZnO/CdSe, and FeOOH/CdS were successfully made using this hydrothermal strategy and the growth mechanism was also discussed. This hydrothermal strategy is simple and green, and can be extended to the synthesis of various 1D composite structures. Moreover, the interaction between the shell nanoparticles and the one-dimensional nanomaterials were confirmed by photoluminescence investigation of ZnO/CdS.

Theoretical research on onedimensional charge density wave systems is outlined. A simple coupled electron-photon Hamiltonian is studied including a Green's function approach, molecular dynamics, and Monte Carlo path integral method. As in superconductivity, the nonperturbative nature of the system makes the physical ground states and low energy excitations drastically different from the bare electrons and phonons. Solitons carry quantum numbers which are entirely different from those of the bare electrons and holes. The fractional charge character of the solitons is an example of this fact. Solitons are conveniently generated by doping material with donors or acceptors or by photon absorption. Most predictions of the theory are in qualitative agreement with experiments. The onedimensional charge density wave system has potential technological importance and a possible role in uncovering phenomena which might have implications in relativistic field theory and elementary particle physics

We investigate the thermoelectric properties of one-dimensional (1D) graphene antidot arrays by nonequilibrium Green's function method. We show that by introducing antidots to the pristine graphene nanoribbon the thermal conductance can be reduced greatly while keeping the power factor still high, thus leading to an enhanced thermoelectric figure of merit (ZT). Our numerical results indicate that ZT values of 1D antidot graphene arrays can be up to unity, which means the 1D graphene antidot arrays may be promising for thermoelectric applications. -- Highlights: ► We study thermoelectric properties of one-dimensional (1D) graphene antidot arrays. ► Thermoelectric figure of merit (ZT) of 1D antidot arrays can exceed unity. ► ZT of 1D antidot arrays is larger than that of two-dimensional arrays.

The scattering theory is treated for the one-dimensional Dirac equation with potentials that are bounded, measurable, real-valued functions on the real line, having constant values, not necessarily the same, on the left and on the right side of a compact interval. Such potentials appear in the Klein paradox. It is shown that appropriately modified wave operators exist and that the corresponding S-operator is unitary. The connection between time-dependent scattering theory and time-independent scattering theory in terms of incoming and outgoing plane wave solutions is established and some further properties are proved. All results and their proofs have a straightforward translation to the one-dimensional Schroedinger equation with the same class of step potentials

Full Text Available Brazil has played an important role in the development and use of resonance Raman spectroscopy as a powerful characterization tool for materials science. Here we present a short history of Raman scattering research in Brazil, highlighting the important contributions to the field coming from Brazilian researchers in the past. Next we discuss recent and important contributions where Brazil has become a worldwide leader, that is on the physics of quasi-onedimensional carbon nanotubes. We conclude this article by presenting results from a very recent resonance Raman study of exciting new materials, that are strictly one-dimensional carbon chains formed by the heat treatment of very pure double-wall carbon nanotube samples.

A program for the solution of the time- and space-dependent multigroup diffusion equations and the delayed-neutron precursors concentration equations in onedimensional geometries by the weighted residual method is described. The discretized equations are solved through an iterative procedure with convergence accelerated by the over-relaxation method. The system is perturbed through the variation of the nuclide concentrations in specified regions. Two feedback effects are included, namely, the temperature and the burnup. (Author) [pt

The propagation of finite-amplitude plane sound in one-dimensional layered media is studied by the extended method of transfer matrix formalism. For the periodic layered system consisting of two alternate types of liquid, the energy distribution and the phase vectors of the interface vibration are computed and analyzed. It is found that in the pass-band, the second harmonic of sound wave can propagate with the characteristic modulation

During March 1973 the European American Committee on Reactor Physics proposed a series of simple one-dimensional reactor kinetics problems, with the intention of comparing the relative efficiencies of the numerical methods employed in various codes, which are currently in use in many national laboratories. This report reviews the contributions submitted to this benchmark exercise and attempts to assess the relative merits and drawbacks of the various theoretical and computer methods. (author)

Effects of non-uniform heating in the core and additional forced circulation during decay heat removal operation are studied with a simplified mixed convection loop. The heat transfer coefficient is calculated analytically and measured experimentally. The analytic solution obtained from a one-dimensional heat equation is found to agree well with the experimental results. The effects of the non-uniform heating and the forced circulation are discussed

We consider the energy density and energy transfer in small amplitude, one-dimensional waves on a string and find that the common expressions used in textbooks for the introductory physics with calculus course give wrong results for some cases, including standing waves. We discuss the origin of the problem, and how it can be corrected in a way appropriate for the introductory calculus-based physics course. (letters and comments)

Intermittent filamentary dynamics of the vortex matter in superconductors is found in films of YBa_{2}Cu_{3}O_{7-δ} deposited on tilted substrates. Deposition of this material on such substrates creates parallel channels of easy flux penetration when a magnetic field is applied perpendicular to the film. As the applied field is gradually increased, magneto-optical imaging reveals that flux penetrates via numerous quasi-one-dimensional jumps. The distribution of flux avalanche sizes follows a ...

The motion of an electron in a random one-dimensional spatially correlated potential is investigated. The spatial correlation is generated by a Markov chain. It is shown that the influence of the spatial correlation can be described by means of oscillating vertices usually neglected in the Berezinskii diagram technique. Correlation mainly leads to an increase of the localization length in comparison with an uncorrelated potential. However, there is a region of the parameter, where the localization decreases. (author)

In order to motivate their group members to perform certain tasks, leaders use different leadership styles. These styles are based on leaders' backgrounds, knowledge, values, experiences, and expectations. The one-dimensional styles, used by many world leaders, are autocratic and democratic styles. These styles lie on the two opposite sides of the leadership spectrum. In order to precisely define the leadership styles on the spectrum between the autocratic leadership style and the democratic ...

We consider the low-temperature limit of the long-distance asymptotic behavior of the finite temperature density-density correlation function in the one-dimensional Bose gas derived recently in the algebraic Bethe Ansatz framework. Our results confirm the predictions based on the Luttinger liquid and conformal field theory approaches. We also demonstrate that the amplitudes arising in this asymptotic expansion at low-temperature coincide with the amplitudes associated with the so-called critical form factors. (orig.)

Background Photonic crystals are artificial structures that have periodic dielectric components with different refractive indices. Under certain conditions, they abnormally refract the light, a phenomenon called negative refraction. Here we experimentally characterize negative refraction in a onedimensional photonic crystal structure; near the low frequency edge of the fourth photonic bandgap. We compare the experimental results with current theory and a theory based on the group velocity de...

It is generally thought that adiabatic exchange of two identical particles is impossible in one spatial dimension. Here we describe a simple protocol that permits adiabatic exchange of two Majorana fermions in a one-dimensional topological superconductor wire. The exchange relies on the concept of "Majorana shuttle" whereby a $\\pi$ domain wall in the superconducting order parameter which hosts a pair of ancillary Majoranas delivers one zero mode across the wire while the other one tunnels in ...

Quantum Walks are unitary processes describing the evolution of an initially localized wavefunction on a lattice potential. The complexity of the dynamics increases significantly when several indistinguishable quantum walkers propagate on the same lattice simultaneously, as these develop non-trivial spatial correlations that depend on the particle's quantum statistics, mutual interactions, initial positions, and the lattice potential. We show that even in the simplest case of a quantum walk on a onedimensional graph, these correlations can be shaped to yield a complete set of compact quantum logic operations. We provide detailed recipes for implementing quantum logic on one-dimensional quantum walks in two general cases. For non-interacting bosons—such as photons in waveguide lattices—we find high-fidelity probabilistic quantum gates that could be integrated into linear optics quantum computation schemes. For interacting quantum-walkers on a one-dimensional lattice—a situation that has recently been demonstrated using ultra-cold atoms—we find deterministic logic operations that are universal for quantum information processing. The suggested implementation requires minimal resources and a level of control that is within reach using recently demonstrated techniques. Further work is required to address error-correction.

Several examples are known in which massive arrays of metal atomic chains are formed on semiconductor surfaces that show quasi-one-dimensional metallic electronic structures. In this review, Au chains on Si(557) and Si(553) surfaces, and In chains on Si(111) surfaces, are introduced and discussed with regard to the physical properties determined by experimental data from scanning tunneling microscopy (STM), angle-resolved photoemission spectroscopy (ARPES) and electrical conductivity measurements. They show quasi-one-dimensional Fermi surfaces and parabolic band dispersion along the chains. All of them are known from STM and ARPES to exhibit metal-insulator transitions by cooling and charge-density-wave formation due to Peierls instability of the metallic chains. The electrical conductivity, however, reveals the metal-insulator transition only on the less-defective surfaces (Si(553)-Au and Si(111)-In), but not on a more-defective surface (Si(557)-Au). The latter shows an insulating character over the whole temperature range. Compared with the electronic structure (Fermi surfaces and band dispersions), the transport property is more sensitive to the defects. With an increase in defect density, the conductivity only along the metal atomic chains was significantly reduced, showing that atomic-scale point defects decisively interrupt the electrical transport along the atomic chains and hide the intrinsic property of transport in quasi-one-dimensional systems.

It has been pointed out by Bekenstein and Mayo that the behavior of the black hole's entropy or information flow is similar to information flow through one-dimensional channel. Here I analyze the same issue with the use of gravitational anomalies. The rate of the entropy change (S) and the power (P) of the Hawking emission are calculated from the relevant components of the anomalous stress tensor under the Unruh vacuum condition. I show that the dependence of S on the power is S ∝ P{sup 1/2}, which is identical to that for the information flow in a one-dimensional system. This is established by using the (1+1)-dimensional gravitational anomalies first. Then the fact is further bolstered by considering the (1+3)-dimensional gravitational anomalies. It is found that, in the former case, the proportionality constant is exactly identical to the one-dimensional situation, known as Pendry's formula, while in the latter situation its value decreases. (orig.)

In this paper one deals with the theoretical derivation of energy bands and of related wavefunctions characterizing quasi 1D semiconductor heterostructures, such as InAs quantum wire models. Such models get characterized this time by equal coupling strength superpositions of Rashba and Dresselhaus spin-orbit interactions of dimensionless magnitude a under the influence of in-plane magnetic fields of magnitude B. We found that the orientations of the field can be selected by virtue of symmetry requirements. For this purpose one resorts to spin conservations, but alternative conditions providing sensible simplifications of the energy-band formula can be reasonably accounted for. Besides the wavenumber k relying on the 1D electron, one deals with the spin-like s=±1 factors in the front of the square root term of the energy. Having obtained the spinorial wavefunction, opens the way to the derivation of spin precession effects. For this purpose one resorts to the projections of the wavenumber operator on complementary spin states. Such projections are responsible for related displacements proceeding along the Ox-axis. This results in a 2D rotation matrix providing both the precession angle as well as the precession axis

We present an analytical study of one-dimensional semiconductor superlattices in external electric fields, which may be time dependent. A number of general results for the (quasi)energies and eigenstates are derived. An equation of motion for the density matrix is obtained for a two-band model...

In a recent paper we give the first rigorous derivation of the celebrated Ginzburg-Landau (GL) theory, starting from the microscopic Bardeen-Cooper-Schrieffer (BCS) model. Here we present our results in the simplified case of a one-dimensional system of particles interacting via a delta-potential....

It is demonstrated that in a one-dimensional Ising chain with nearest-neighbor interactions, irradiated by a weak resonant transverse field, a stimulated wave of flipped spins can be triggered by a flip of a single spin. This analytically solvable model illustrates mechanisms of quantum amplification and quantum measurement

A one-dimensional mathematical model is derived for a three-stage pulse-tube refrigerator (PTR) that is based on the conservation laws and the ideal gas law. The three-stage PTR is regarded as three separate single-stage PTRs that are coupled via proper junction conditions. At the junctions there

A one-dimensional mathematical model is derived for a three-stage pulse-tube refrigerator (PTR) that is based on the conservation laws and the ideal gas law. The three-stage PTR is regarded as three separate single-stage PTRs that are coupled via proper junction conditions. At the junctions there

The objective of this paper is to assess the core degradation modeling in COMPASS code by simulating the PHEBUS FPT3 experiment. For the comparison purpose, the numerical simulation by using MELCOR 2.1 have also conducted for the FPT3 experiment. Consequently, COMPASS results of PHEBUS FPT3 have been compared with the experimental data and MELCOR results. For the purpose of COMPASS code validation, the numerical simulation for PHEBUS FPT3 experiment has been conducted. The temperature of the main component has been secured by using COMPASS code for a fuel, cladding, control rod and surrounding structure. And they are compared with that of experimental data as well as MELCOR simulation results. MELCOR are showing that an oxidational reaction starts a little bit earlier time and has the slightly higher value of the accumulated hydrogen mass, while COMPASS code predicts the slightly lower value of the accumulated hydrogen mass.

An analysis of the passive dynamic walking of a compass-gait biped model under the OGY-based control approach using the impulsive hybrid nonlinear dynamics is presented in this paper. We describe our strategy for the development of a simplified analytical expression of a controlled hybrid Poincaré map and then for the design of a state-feedback control. Our control methodology is based mainly on the linearization of the impulsive hybrid nonlinear dynamics around a desired nominal one-periodic hybrid limit cycle. Our analysis of the controlled walking dynamics is achieved by means of bifurcation diagrams. Some interesting nonlinear phenomena are displayed, such as the period-doubling bifurcation, the cyclic-fold bifurcation, the period remerging, the period bubbling and chaos. A comparison between the raised phenomena in the impulsive hybrid nonlinear dynamics and the hybrid Poincaré map under control was also presented.

Full Text Available ) of these atoms. The enthalpy of formation of free atoms, ∆H, can vary between 50 and 300 kJ.mol-1. The ballast is heated by radiation from the walls of the atomizer, and its temperature, Tb(t), is given by the Stefan-Boltzmann equation [4...]: ( )ddt T t E c S V T t T tb b b b b b w b( ) ( ) ( )= −σ ρ 4 4 (6) where σ is the Stefan-Boltzmann constant, Eb is the emissivity of the ballast material, cb is the heat capacity of ballast material, ρb...