Refine
Year of publication
Document Type
- Preprint (1037) (remove)
Language
- English (1037) (remove)
Keywords
- AG-RESY (14)
- Approximation (9)
- Case-Based Reasoning (9)
- RODEO (9)
- Mehrskalenanalyse (8)
- Wavelet (8)
- Boltzmann Equation (7)
- Location Theory (7)
- Numerical Simulation (7)
- Case Based Reasoning (6)
Faculty / Organisational entity
- Kaiserslautern - Fachbereich Mathematik (554)
- Kaiserslautern - Fachbereich Informatik (290)
- Kaiserslautern - Fachbereich Physik (157)
- Fraunhofer (ITWM) (19)
- Kaiserslautern - Fachbereich Maschinenbau und Verfahrenstechnik (9)
- Kaiserslautern - Fachbereich Elektrotechnik und Informationstechnik (6)
- Kaiserslautern - Fachbereich Sozialwissenschaften (1)
- Kaiserslautern - Fachbereich Wirtschaftswissenschaften (1)
Fast reconstruction formulae in x-ray computerized tomography demand the directions, in which the measurements are taken, to be equally distributed over the whole circle. In many applications data can only be provided in a restricted range. Here the intrinsic difficulties are studied by giving a singular value decomposition of the Radon transform in a restricted range. Practical limitations are deduced.
The Trippstadt Problem
(1984)
Close to Kaiserslautern is the town of Trippstadt, which, together with five other small towns forms a local administration unit (Verbandsgemeinde) called Kaiserslautern-Süd. Trippstadt has its own beautiful public swimming pool, which causes problems though; the cost for the upkeep of the pool is higher than the income and thus has to be divided among the towns belonging to the Verbandsgemeinde. Because of this problem the administration wanted to find out which fraction of the total number of pool visitors came from the different towns. They planned to ask each pool guest where he came from. They did this for only three days though because the waiting lines at the cashiers became unbearably long and they could see that because of this the total number of guests would decrease. Then they wondered how to find a better method to get the same data and that was when I was asked to help with the solution of the problem.
In these notes we will discuss some aspects of a problem arising in carindustry. For the sake of clarity we will set the problem into an extremely simplified scheme. Suppose that we have a body which is emitting sound, and that the sound is measured at a finite number of points around the body. We wish to determine the intensity of the sound at an observation point which is moving.
We report on the exchange bias effect as a function of the in-plane direction of the applied field in two-fold symmetric, epitaxial Ni80Fe20/Fe50Mn50 bilayers grown on Cu(110) single crystal substrates. An enhancement of the exchange bias field, Heb, up to a factor of two is observed if the external field is nearly, but not fully aligned perpendicular to the symmetry direction of the exchange bias field. From the measurement of the ex-change bias field as a function of the in-plane angle of the applied field, the unidirectional, uniaxial and four-fold anisotropy contributions are determined with high precision. The symmetry direction of the unidirec-tional anisotropy switches with increasing NiFe thickness from [110] to [001].
We want to study solid objects in real three dimensional space aiming at two issues:; (i1) modelling solids subject to boolean set algebra, including wire models,; (i2) determining the behaviour of moving solids, e.g. when they collide and the resulting points of contact.; ; This research has been initiated by the FORD Motor Company, Cologne. It is motivated by the intention to provide for a model of an automatical car gear, which gives a high precision basis to the optimization of moving tolerances.
Estimation of P(R kl/gleich S) is considered for the simple stress-strength model of failure. Using the Pareto and Power distributions together with their combined form a useful parametric solution is obtained and is illustrated numerically. It is shown that these models are also applicable when only the tails of distributions for R and S are considered. An application to the failure study concerning the fractures is also included.
Stability and Robustness Properties of Universal Adaptive Controllers for First Order Linear Systems
(1987)
The question: What is an adaptive controller? is as old as the word adaptive control itself. In this paper we will adopt a pragmatic viewpoint which identifies adaptive controllers with nonlinear feedback controllers, designed for classes (families) of linear systems. In contrast to classical linear feedback controllers which are designed for individual systems, these non-linear controllers are required to achieve a specific design objective (such as e.g. stability, tracking or decoupling) for a whole prescribed family of linear systems.
Patterns are considered as normalized measures and distances between them are defined as distances of the corresponding measures using metrics in measure spaces. This idea can be applied for pattern recognition if smeared patterns have to be compared with given ideal patterns. Different metrics are sensitive to different characteristics of the patterns - this is demonstrated in discussing examples. Particular attention is paid to a problem of Quality Control for an artificial fabric, where the distance to uniformity is defined and evaluated; the results are now used in industry.
As shown by Krasnosel" skii, the classical Preisach model allows to construct a hysteresis operator Wbetween spaces of real functions of time. This construction, via the definition of a measure mü in the so-called Preisach plane, is recalled. Characterizations in terms of mü are given for several mapping and continuity properties of W in various function spaces, for the invertibility of W and for the corresponding mapping and continuity properties of the inverse.
The performance of a combustion engine is essentially determined by the charge cycle, i.e. by the inflow of fresh air through the inlet pipe into the cylinder after a combustion cycle. The amount of air, exchanged during this process, depends on many factors, e.g. the number of revolutions per minute, the temperature, the engine and valve geometry. In order to have a tool in designing the engine one is interested in calculating this amount. The proper calculation would involve the solution of three-dimensional hydrodynamical equations governing the gas flow including chemical reactions in a complicated geometry, consisting of the cylinder, valves, inlet and outlet pipe. Since this is clearly too ambitious, we consider a simplified model.
As an alternative to the commonly used Monte Carlo Simulation methods for solving the Boltzmann equation we have developed a new code with certain important improvements. We present results of calculations on the reentry phase of a space shuttle. One aim was to test physical models of internal energies and of gas-surface interactions.
Industrial mathematics has many faces; but its essential feature is the cooperation of partners - from industry and from universities - with quite different interest (business versus academic carreer), normally working on different time scales. They measure success in a different way (selling rate against citing index), they have different hierarchies of values and are very often distrusting each other. Industry doubts that mathematicians are willing and/or able to produce something real practical and useful (and the mathematicians should not be too much surprised about this attitude, they very often doubt themselves) - mathematicians are afraid to loose their competence (their ideal of scientific truth, to say it more idealistically), to sell their souls.
Special technological applications like the construction of a dental attachment require structural parts which have very small operall dimensions. Very often these parts are subjected to high loadings. The failure of a small spring was the starting point for an investigation with the aim to design a suitable new spring shape.
We consider universal adaptive stabilization and tracking controllers for classes of linear systems. Under the technical assumption of linear scaling invariance necessary and sufficient conditions for adaptive stabilization are derived. For scalar systems sufficient conditions for adaptive tracking of finite dimensional reference signals are explored.
We present the concept of a universal adaptive tracking controller for classes of linear systems. For the class of scalar minimum phase systems of relative degree one, adaptive tracking is shown for arbitrary finite dimensional reference signals. The controller requires no identificaiton of the system parameters. Robustness properties are explored.
On the Moving Preisach Model
(1990)
Fleeces made from artificial fabric are the basic material for many products, ranging from carpets to napkins. Itturns out that their quality is determined by the distribution of the fibres, which can be measured eithter by the optical transmission properties or by the thickness of the material. In both cases one obtains a 2-dimensional signal and one would like to have an objective quality criterion, based on a suitable analysis of these data, which, moreover, can be automated.; In this paper we propose a solution to this problem, based on multiresolution techniques, which have been developed in image analysis through the last few years. Moreover, we use these techniques to investigate fractal properties of the textures.
We present a deterministic simulation scheme for the Boltzmann Semiconductor Equation. The convergence of the method is shown for a simplified space homogeneous case. Numerical experiments, which are very promising, are also given in this situation. The extension for the application to the space inhomogeneous equation with a self consistent electric field is quoted. Theoretical considerations in that case are in preparation.
This paper contains the basic ideas and practical aspects for numerical methods for solving the Boltzmann Equation. The main field of application considered is the reentry of a Space Shuttle in the transition from free molecular flow to continuum flow. The method used will be called Finite Pointset Method (FPM) approximating the solution by finite sets of particles in a rigorously defined way. Convergence results are cited while practical aspects of the algorithm are emphasized. Ideas for the transition to the Navier Stokes domain are shortly discussed.
Treating polyatomic gases in kinetic gas theory requires an appropriate molecule model taking into account the additional internal structure of the gas particles. In this paper we describe two such models, each arising from quite different approaches to this problem. A simulation scheme for solving the corresponding kinetic equations is presented and some numerical results to 1D shockwaves are compared.
This report contains the following three papers about computations of rarefied gas flows:; ; a) Rarefied gas flow around a disc with different angles of attack, published in the proceedings of the 17th RGD Symposium, Aachen, 1990.; ; b) Hypersonic flow calculations around a 3D-deltawing at low Knudsen numbers, published in the proceedings of the 17th RGD Symposium,; Aachen, 1990.; ; c) Rarefied gas flow around a 3D-deltawing, published in the proceedings of the Workshop on Hypersonic Flows for Reentry Problems,; Part 1, Antibes, France, January 22-25, 1990.; ; All computations are part of the HERMES Research and Development Program.
This article describes the basic concepts of an extensible customizable knowledge-basedgraphical editor and its adoption to the DOCASE methodology and tool environment. Oneaspect in this field is the mapping of conceptual models (expressed in a specific language)to their graphical representations. This also has impacts to the semantic of the user actionsin a graphical editor tool. The ability to extend and customize the editor can be used tobuild specific graphical interfaces to various kinds of tools in the software developmentprocess. Major aspects of ODE are semantics-directed editing besides normal syntax-directed editing, support of abstraction mechanisms, multiple modeless views to attack com-plexity, semantic analization and animation. The result is an highly customizable graphicaleditor construction set that matches requirements of applications in many domains of systemdesign.
Based on the experiences from an autonomous mobile robot project called MOBOT-III, we found hard realtime-constraints for the operating- system-design. ALBATROSS is "A flexible multi-tasking and realtime network-operating-system-kernel". The focusin this article is on a communication-scheme fulfilling the previous demanded assurances. The centralchapters discuss the shared buffer management and the way to design the communication architecture.Some further aspects beside the strict realtime-requirements like the possibilities to control and watch a running system, are mentioned. ALBATROSS is actually implemented on a multi-processor VMEbus-system.
We have presented here a two-dimensional kinetical scheme for equations governing the motion of a compressible flow of an ideal gas (air) based on the Kaniel method. The basic flux functions are computed analytically and have been used in the organization of the flux computation. The algorithm is implemented and tested for the 1D shock and 2D shock-obstacle interaction problems.
The paper presents a parallelization technique for the finite pointset method, a numerical method for rarefied gas flows.; First we give a short introduction to the Boltzmann equation, which describes the behaviour of rarefied gas flows. The basic ideas of the finite pointset method are presented and a strategy to parallelize the algorithm will be explained. It is shown that a static processor partition leads to an insufficient load-balance of the processors. Therefore an optimized parallelization technique based on an adaptive processor partition will be introduced, which improves the efficiency of the simulation code over the whole region of interesting flow situation. Finally we present a comparison of the CPU-times between a parallel computer and a vector computer.
Using particle methods to solve the Boltzmann equation for rarefied gases numerically, in realistic streaming problems, huge differences in the total number of particles per cell arise. In order to overcome the resulting numerical difficulties the application of a weighted particle concept is well-suited. The underlying idea is to use different particle masses in different cells depending on the macroscopic density of the gas. Discrepance estimates and numerical results are given.
The wave equation with a Preisach hysteresis operator can be considered as a one-dimensional projection of Maxwell" s equations in a ferromagnetic medium. An initial-boundary value problem for this equation is solved here with emphasizing the fact that under a bounded forcing term the solutions remain bounded. This is due to the strong dissipation of hysteresis energies. New proofs of hysteresis energy inequalities are given without referring to the structure of hysteresis memory.
The efficient numerical treatment of the Boltzmann equation is a very important task in many fields of application. Most of the practically relevant numerical schemes are based on the simulation of large particle systems that approximate the evolution of the distribution function described by the Boltzmann equation. In particular, stochastic particle systems play an important role in the construction of various numerical algorithms.
Retrieval of cases is one important step within the case-based reasoning paradigm. We propose an improvement of this stage in the process model for finding most similar cases with an average effort of O[log2n], n number of cases. The basic idea of the algorithm is to use the heterogeneity of the search space for a density-based structuring and to employ this precomputed structure, a k-d tree, for efficient case retrieval according to a given similarity measure sim. In addition to illustrating the basic idea, we present the expe- rimental results of a comparison of four different k-d tree generating strategies as well as introduce the notion of virtual bounds as a new one that significantly reduces the retrieval effort from a more pragmatic perspective. The presented approach is fully implemented within the (Patdex) system, a case-based reasoning system for diagnostic applications in engineering domains.
We consider a transmission boundary-value problem for the time-harmonic Maxwell equations neglecting displacement currents. The usual transmission conditions, which require the continuity of the tangential components of the electric and magnetic fields across boundaries are slightly modified. For this new problem we show that the uniqueness of the solution depends on the topological properties of the domains under consideration. Finally we obtain existence results by using a boundary integral equation approach.
We consider a transmission boundary-value problem for the time-harmonic Maxwell equations without displacement currents. As transmission conditions we use the continuity of the tangential parts of the magnetic field H and the continuity of the normal components of the magnetization B=müH. This problem, which is posed over all IR3, is then restricted to a bounded domain by introducing artificial boundary conditions. We present uniqueness and existence proofs for this problem using an integral equation approach and compare the results with those obtained in the unbounded case.
We consider two transmission boundary-value problems for the time-harmonic Maxwell equations without displacement currents. For the first problem we use the continuity of the tangential parts of the electric and magnetic fields across material discontinuities as transmission conditions. In the second case the continuity of the tangential components of the electric field E is replaced by the continuity of the normal component of the magnetization B=müH. For this problem existence of solutions is already shown in [6]. If the domains under consideration are not simply connected the solution is not unique. In this paper, we improve the regularity results obtained in [6] and then prove existence and uniqueness theorems for the first problem by extracting its solution out of the set of all solutions of the second problem. Thus we establish a connection between the solutions corresponding to the different transmission boundary conditions.
In this paper noises and disturbances are treated as distributions of some general class. The problem of sensitivity minimization is considered. A design procedure for the construction of Luenberger observers which estimate the state of a system with a given rate of accuracy has been proposed. The design procedure is applied to identify the first derivatives of an oscillating signal. The constraints on a noise and on a sampling which are necessary to estimate the derivatives to a given accuracy have been obtained.
A multiparameter, polynomial feedback strategy is introduced to solve the universal adapative tracking problem for a class of multivariable minimum phase system and reference signals generated by a known linear time-invariant differential equation. For 2-input, 2-output, minimum phase systems (A,B,C) with det(CB)0, a different polynomial tracking controller is given which does not invoke a spectrum unmixing set.
Several topological necessary conditions of smooth stabilization in the large have been obtained. In particular, if a smooth single-input nonlinear system is smoothly stabilizable in the large at some point of a connected component of equilibria set, then the connected component is to be an unknoted, unbounded curve.
Given a proper antistable rational transfer function g, a balanced realization of g is contructed as a matrix representation of the abstract shift realization introduced in Fuhrmann [1976]. The required basis is constructed as a union of sets of polynomials orthogonal with respect to weights given by the square of the absolute values of minimal degree Schmidt vectors of the corresponding Hankel operators. This extends results of Fuhrmann [1991], obtained in the generic case.
The polynomial approach introduced in Fuhrmann [1991] is extended to cover the crucial area of AAK theory, namely the characterization of zero location of the Schmidt vectors of the Hankel operators. This is done using the duality theory developed in that paper but with a twist. First we get the standard, lower bound, estimates on the number of unstable zeroes of the minimal degree Schmidt vectors of the Hankel operator. In the case of the Schmidt vector corresponding to the smallest singular the lower bound is in fact achieved. This leads to a solution of a Bezout equation. We use this Bezout equation to introduce another Hankel operator which have singular values that are the inverse of the singular values of the original Hankel operator.
Diffeomorphisms are given between different subsets of linear systems of fixed McMillan degree. The sets considered are the set of all systems of fixed McMillan degree, the subset of stable systems, the subset of bounded real systems, the subset of positive real systems, the subset of stable systems with Hankel singular values bounded by one. State space techniques are used in the proofs.
The paper presents theoretical and numerical investigations on simulation methods for the Boltzmann equation with axisymmetric geometry. The main task is to reduce the computational effort by taking advantage of the symmetry in the solution of the Boltzmann equation.; The reduction automatically leads to the concept of weighting functions for the radial space coordinate and therefore to a modified Boltzmann equation. Consequently the classical simulation methods have to be modified according to the new equation.; The numerical results shown in this paper - rarefied gas flows around a body with axisymmetric geometry - were done in the framework of the European space project HERMES.
Elements of the differential topology are used to prove necessary conditions for stabilizability in large by a smooth feedback. Criteria for the smooth feedback stabilizing a smooth nonlinear system locally to have the smooth piecewise smooth extension, which stabilizes the system over a given compact set, have been obtained.
A method of decoupling normalizing transformations has been developed. According to the method only the part of differential equations corresponding to the dynamic on a center manifold has to be modified by means of the normalizing transformations of a Poincare type. The existence of the normalizing transformation completely decoupling the stable dynamic from the center manifold dynamic has been proved. A numerical procedure for the calculation of asymptotic series for the decoupling normalizing transformation has been proposed. The developed method is especially important for the perturbation theory of center manifold and, in particular, for the local stabilization theory. In the paper some sufficient conditions for local stabilization have been given.
On the Mróz Model
(1992)
This report contains a collection of abstracts for talks given at the "Deduktionstreffen" held at Kaiserslautern, October 6 to 8, 1993. The topics of the talks range from theoretical aspects of term rewriting systems and higher order resolution to descriptions of practical proof systems in various applications. They are grouped together according the following classification: Distribution and Combination of Theorem Provers, Termination, Completion, Functional Programs, Inductive Theorem Proving, Automatic Theorem Proving, Proof Presentation. The Deduktionstreffen is the annual meeting of the Fachgruppe Deduktionssysteme in the Gesellschaft für Informatik (GI), the German association for computer science.
Questions arising from Statistical Decision Theory, Bayes Methods and other probability theoretic fields lead to concepts of orthogonality of a family of probability measures. In this paper we therefore give a sketch of a generalized information theory which is very helpful in considering and answering those questions. In this adapted information theory Shannon's classical transition channels modelled by finite stochastic matrices are replaced by compact families of probability measures that are uniformly integrable. These channels are characterized by concepts such as information rate and capacity and by optimal priors and the optimal mixture distribution. For practical studies we introduce an algorithm to calculate the capacity of the whole probability family which is appli cable even for general output space. We then explain how the algorithm works and compare its numerical costs with those of the classical Arimoto-Blahut-algorithm.
The system of shallow water waves is one of the classical examples for nonlinear, twodimensional conservation laws. The paper investigates a simple kinetic equation depending on a parameter e which leads for e to 0 to the system of shallow water waves. The corresponding equilibrium distribution function has a compact support which depends on the eigenvalues of the hyperbolic system. It is shown that this kind of kinetic approach is restricted to a special class of nonlinear conservation laws. The kinetic model is used to develop a simple particle method for the numerical solution of shallow water waves. The particle method can be implemented in a straightforward way and produces in test examples sufficiently accurate results.
SPIN-NFDS Learning and Preset Knowledge for Surface Fusion - A Neural Fuzzy Decision System -
(1993)
The problem to be discussed in this paper may be characterized in short by the question: "Are these two surface fragments belonging together (i.e. belonging to the same surface)?" The presented techniques try to benefit from some predefined knowledge as well as from the possibility to refine and adapt this knowledge according to a (changing) real environment, resulting in a combination of fuzzy-decision systems and neural networks. The results are encouraging (fast convergence speed, high accuracy), and the model might be used for a wide range of applications. The general frame surrounding the work in this paper is the SPIN- project, where emphasis is on sub-symbolic abstractions, based on a 3-d scanned environment.
This paper refers to the problem of adaptability over an infinite period of time, regarding dynamic networks. A never ending flow of examples have to be clustered, based on a distance measure. The developed model is based on the self-organizing feature maps of Kohonen [6], [7] and some adaptations by Fritzke [3]. The problem of dynamic surface classification is embedded in the SPIN project, where sub-symbolic abstractions, based on a 3-d scanned environment is being done.
Case-based problem solving can be significantly improved by applying domain knowledge (in opposition to problem solving knowledge), which can be acquired with reasonable effort, to derive explanations of the correctness of a case. Such explanations, constructed on several levels of abstraction, can be employed as the basis for similarity assessment as well as for adaptation by solution refinement. The general approach for explanation-based similarity can be applied to different real world problem solving tasks such as diagnosis and planning in technical areas. This paper presents the general idea as well as the two specific, completely implemented realizations for a diagnosis and a planning task.
The paper presents a fast implementation of a constructive method to generate a special class of low-discrepancy sequences which are based on Van Neumann-Kakutani tranformations. Such sequences can be used in various simulation codes where it is necessary to generate a certain number of uniformly distributed random numbers on the unit interval.; From a theoretical point of view the uniformity of a sequence is measured in terms of the discrepancy which is a special distance between a finite set of points and the uniform distribution on the unit interval.; Numerical results are given on the cost efficiency of different generators on different hardware architectures as well as on the corresponding uniformity of the sequences. As an example for the efficient use of low-discrepancy sequences in a complex simulation code results are presented for the simulation of a hypersonic rarefied gas flow.
This paper considers the numerical solution of a transmission boundary-value problem for the time-harmonic Maxwell equations with the help of a special finite volume discretization. Applying this technique to several three-dimensional test problems, we obtain large, sparse, complex linear systems, which are solved by using BiCG, CGS, BiCGSTAB resp., GMRES. We combine these methods with suitably chosen preconditioning matrices and compare the speed of convergence.
In these lectures we will mainly treat a billard game. Our particles will be hard spheres. Not always: We will also touch cases, where particles have interior energies due to rotation or vibration, which they exchange in a collision, and we will talk about chemical reactions happening during a collision. But many essential aspects occur already in the billard case which will be therefore paradigmatic. I do not know enough about semiconductors to handle collisions there - the Boltzmann case is certainly different but may give some idea even for the other cases.
Discrete families of functions with the property that every function in a certain space can be represented by its formal Fourier series expansion are developed on the sphere. A Fourier series type expansion is obviously true if the family is an orthonormal basis of a Hilbert space, but it also can hold in situations where the family is not orthogonal and is overcomplete. Furthermore, all functions in our approach are axisymmetric (depending only on the spherical distance) so that they can be used adequately in (rotation) invariant pseudodifferential equations on the frames (ii) Gauss- Weierstrass frames, and (iii) frames consisting of locally supported kernel functions. Abel-Poisson frames form families of harmonic functions and provide us with powerful approximation tools in potential theory. Gauss-Weierstrass frames are intimately related to the diffusion equation on the sphere and play an important role in multiscale descriptions of image processing on the sphere. The third class enables us to discuss spherical Fourier expansions by means of axisymmetric finite elements.
Spline functions that interpolate data given on the sphere are developed in a weighted Sobolev space setting. The flexibility of the weights makes possible the choice of the approximating function in a way which emphasizes attributes desirable for the particular application area. Examples show that certain choices of the weight sequences yield known methods. A pointwise convergence theorem containing explicit constants yields a useable error bound.
Exact Solutions of Discrete Kinetic Models and Stationary Problems for the Plane Broadwell Model
(1993)
The paper presents the shuffle algorithm proposed by Baganoff, which can be implemented in simulation methods for the Boltzmann equation to simplify the binary collision process. It is shown that the shuffle algborithm is a discrete approximation of an isotropic collision law. The transition probability as well as the scattering cross section of the shuffle algorithm are opposed to the corresponding quantities of a hard-sphere model. The discrepancy between measures on a sphere is introduced in order to quantify the approximation error by using the shuffle algorithm.
In this paper, we deal with the problem of spherical interpolation of discretely given data of tensorial type. To this end, spherical tensor fields are investigated and a decomposition formula is described. Tensor spherical harmonics are introduced as eigenfunctions of a tensorial analogon to the Beltrami operator and discussed in detail. Based on these preliminaries, a spline interpolation process is described and error estimates are presented. Furthermore, some relations between the spline basis functions and the theory of radial basis functions are developed.
We discuss how kinetic and aerodynamic descriptions of a gas can be matched at some prescribed boundary. The boundary (matching) conditions arise from requirement that the relevant moments (p,u,...) of the particle density function be continuous at the boundary, and from the requirement that the closure relation, by which the aerodynamic equations (holding on one side of the boundary) arise from the kinetic equation (holding on the other side), be satisfied at the boundary. We do a case study involving the Knudsen gas equation on one side and a system involving the Burgers equation on the other side in section 2, and a discussion for the coupling of the full Boltzmann equation with the compressible Navier-Stokes equations in section 3.
Simulation methods like DSMC are an efficient tool to compute rarefied gas flows. Using supercomputers it is possible to include various real gas effects like vibrational energies or chemical reactions in a gas mixture. Nevertheless it is still necessary to improve the accuracy of the current simulation methods in order to reduce the computational effort. To support this task the paper presents a comparison of the classical DSMC method with the so called finite Pointset Method. This new approach was developed during several years in the framework of the European space project HERMES. The comparison given in the paper is based on two different testcases: a spatially homogeneous relaxation problem and a 2-dimensional axisymmetric flow problem at high Mach numbers.
This paper is concerned with the development of a self-adaptive spatial descretization for PDEs using a wavelet basis. A Petrov-Galerkin method [LPT91] is used to reduce the determination of the unknown at the new time step to the computation of scalar products. These have to be discretized in an appropriate way. We investigate this point in detail and devise an algorithm that has linear operation count with respect to the number of unknowns. It is tested with spline wavelets and Meyer wavelets retaining the latter for their better localisation at finite precision. The algorithm is then applied to the one dimensional thermodiffusive equations. We show that the adaption strategy merits to be modified in order to take into account the particular and very strong nonlinearity of this problem. Finally, a supplementary Fourier discretization permits the computation of two dimensional flame fronts.
Abstract: The classification of quasi - primary fields is outlined. It is proved that the only conserved quasi - primary currents are the energy - momentum tensor and the O(N)-Noether currents. Derivation of all quasi - primary fields and the resolution of degeneracy is sketched. Finally the limits d = 2 and d = 4 of the space dimension are discussed. Whereas the latter is trivial the former is only almost so. (To appear in the Proceedings of the XXII Conference on Differential Geometry Methods in Theoretical Physics, Ixtapa, Mexico, September 20-24, 1993)
Based on experiences from an autonomous mobile robot project called MOBOT -III, we found hard realtime-constraints for the operating-system-design. ALBATROSS is "A flexible multi-tasking and realtime network-operatingsystem-kernel", not limited to mobile- robot-projects only, but which might be useful also wherever you have to guarantee a high reliability of a realtime-system. The focus in this article is on a communication-scheme fulfilling the demanded (hard realtime-) assurances although not implying time-delays or jitters on the critical informationchannels. The central chapters discuss a locking-free shared buffer management, without the need for interrupts and a way to arrange the communication architecture in order to produce minimal protocol-overhead and short cycle-times. Most of the remaining communication-capacity (if there is any) is used for redundant transfers, increasing the reliability of the whole system. ALBATROSS is actually implemented on a multi-processor VMEbus-system.
We study deterministic conditional rewrite systems, i.e. conditional rewrite systemswhere the extra variables are not totally free but 'input bounded'. If such a systemR is quasi-reductive then !R is decidable and terminating. We develop a critical paircriterion to prove confluence if R is quasi-reductive and strongly deterministic. In thiscase we prove that R is logical, i.e./!R==R holds. We apply our results to proveHorn clause programs to be uniquely terminating.This research was supported by the Deutsche Forschungsgemeinschaft, SFB 314, Project D4
We investigate restricted termination and confluence properties of term rewritADing systems, in particular weak termination and innermost termination, and theirinterrelation. New criteria are provided which are sufficient for the equivalenceof innermost / weak termination and uniform termination of term rewriting sysADtems. These criteria provide interesting possibilities to infer completeness, i.e.termination plus confluence, from restricted termination and confluence properADties.Using these basic results we are also able to prove some new results aboutmodular termination of rewriting. In particular, we show that termination ismodular for some classes of innermost terminating and locally confluent termrewriting systems, namely for nonADoverlapping and even for overlay systems. Asan easy consequence this latter result also entails a simplified proof of the factthat completeness is a decomposable property of soADcalled constructor systems.Furthermore we show how to obtain similar results for even more general cases of(nonADdisjoint) combined systems with shared constructors and of certain hierarADchical combinations of systems with constructors. Interestingly, these modularityresults are obtained by means of a proof technique which itself constitutes a modADular approach.
We present a convenient notation for positive/negativeADconditional equations. Theidea is to merge rules specifying the same function by using caseAD, ifAD, matchAD, and letADexpressions.Based on the presented macroADruleADconstruct, positive/negativeADconditional equational specifiADcations can be written on a higher level. A rewrite system translates the macroADruleADconstructsinto positive/negativeADconditional equations.
A Nonlinear Ray Theory
(1994)
A proof of the famous Huygens" method of wavefront construction is reviewed and it is shown that the method is embedded in the geometrical optics theory for the calculation of the intensity of the wave based on high frequency approximation. It is then shown that Huygens" method can be extended in a natural way to the construction of a weakly nonlinear wavefront. This is an elegant nonlinear ray theory based on an approximation published by the author in 1975 which was inspired by the work of Gubkin. In this theory, the wave amplitude correction is incorporated in the eikonal equation itself and this leads to a sytem of ray equations coupled to the transport equation. The theory shows that the nonlinear rays stretch due to the wave amplitude, as in the work of Choquet-Bruhat (1969), followed by Hunter, Majda, Keller and Rosales, but in addition the wavefront rotates due to a non-uniform distribution of the amplitude on the wavefront. Thus the amplitude of the wave modifies the rays and the wavefront geometry, which in turn affects the growth and decay of the amplitude. Our theory also shows that a compression nonlinear wavefront may develop a kink but an expansion one always remains smooth. In the end, an exact solution showing the resolution of a linear caustic due to nonlinearity has been presented. The theory incorporates all features of Whitham" s geometrical shock dynamics.
Based on normalized coprime factorizations with respect to indefinite metrics and the construction of suitable characteristic functions, the Ober balanced canonical forms for the classes of bounded real and positive real are derived. This uses a matrix representation of the shift realization with respect to a basis related to sets of orthogonal polynomials.
The edge enhancement property of a nonlinear diffusion equation with a suitable expression for the diffusivity is an important feature for image processing. We present an algorithm to solve this equation in a wavelet basis and discuss its one dimensional version in some detail. Sample calculations demonstrate principle effects and treat in particular the case of highly noise perturbed signals. The results are discussed with respect to performance, efficiency, choice of parameters and are illustrated by a large number of figures. Finally, a comparison with a Fourier method and a finite volume method is performed.
Whenever new parts of a car have been developed, the manufacturer needs an estimation of the lifetime of this new part. On one hand the construction must not be too weak, so that the part holds long enough to satisfy the customer, but on the other hand, if the construction is too excessive, the part gets too heavy.; One is interested in methods that only need few measured data from the specimen itself, but use data about the material, because constructing and testing of specimen is expensive.
Monte-Carlo methods are widely used numerical tools in various fields of application, like rarefied gas dynamics, vacuum technology, stellar dynamics or nuclear physics. A central part in all applications is the generation of random variates according to a given probability law. Fundamental techniques to generate non-uniform random variates are the inversion principle or the acceptance-rejection method. Both procedures can be quite time-consuming if the given probability law has a complicated structure.; In this paper we consider probability laws depending on a small parameter and investigate the use of asmptotic expansions to generate random variates. The results given in the paper are restrictedto first order expansions. We show error estimates for the discrepancy as well as for the bounded Lipschitz distance of the asymptotic expansion. Furthermore the integration error for some special classes of functions is given. The efficiency of the method is proved by a numerical example from rarefied gas flows.