Numerical Method for Cost-Weight Optimisation of Stringer-Skin Panels

The need to integrate cost into the early product definition process as an engineering parameter is addressed. The application studied is a fuselage panel that is typical for commercial transport regional jets. Consequently, a semi-empirical numerical analysis using reference data was coupled to model the structural integrity of thin-walled structures with regard to material failure and buckling: skin, stringer, flexural, and interrivet. The optimization process focuses on direct operating cost (DOC) as a function of acquisition cost and fuel burn. It was found that the ratio of acquisition cost to fuel burn was typically 4:3 and that there was a 10% improvement in the DOC for the minimal DOC condition over the minimal weight condition because of the manufacturing cost saving from having a reduced number of larger-area stringers and a slightly thicker skin than that preferred by the minimal weight condition. Also note that the minimal manufacturing cost condition was slightly better than the minimal weight condition, which highlights the key finding: The traditional minimal weight condition is a dated and suboptimal approach to airframe structural design.

Dielectric constant characterization using a numerical method for the microstrip ring resonator

A new method of dielectric-constant measurement is developed. The dielectric constant epsilon(r) RF/microwave substrate is extracted by combining the microstrip ring resonator measurement with Ansoft HFSS electromagnetic simulation software. The developed method has two advantages: (i) characterization of dielectric constant versus multiple frequency points, and (ii) compatibility with electronics design automation (EDA) software tools. This characterization method can reduce the design cycle of microwave circuits and devices. (C) 2004 Wiley Periodicals, Inc.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

NUMERICAL SIMULATIONS OF BURST PROCESSES IN FIBER-BUNDLES

We discuss a very effective numerical method for simulating fibre-bundle models with equal load-sharing and with local load-sharing. Particular attention is paid to the case of the local load-sharing model, in which the critical load x(c) is defined as the average load per fibre that causes the final complete failure. It is shown that x(c) --> 0 when the size of the system N --> infinity. We also show analytically that the power law of the burst size distribution, D(Delta) alpha Delta(-xi), is approximately correct. The exponent xi in the local load-sharing case is not universal, since it depends on the strength distribution as well on as the size of the system.

A Review of Numerical Modelling of Wave Energy Converter Arrays

Large-scale commercial exploitation of wave energy is certain to require the deployment of wave energy converters (WECs) in arrays, creating ‘WEC farms’. An understanding of the hydrodynamic interactions in such arrays is essential for determining optimum layouts of WECs, as well as calculating the area of ocean that the farms will require. It is equally important to consider the potential impact of wave farms on the local and distal wave climates and coastal processes; a poor understanding of the resulting environmental impact may hamper progress, as it would make planning consents more difficult to obtain. It is therefore clear that an understanding the interactions between WECs within a farm is vital for the continued development of the wave energy industry.To support WEC farm design, a range of different numerical models have been developed, with both wave phase-resolving and wave phase-averaging models now available. Phase-resolving methods are primarily based on potential flow models and include semi-analytical techniques, boundary element methods and methods involving the mild-slope equations. Phase-averaging methods are all based around spectral wave models, with supra-grid and sub-grid wave farm models available as alternative implementations.The aims, underlying principles, strengths, weaknesses and obtained results of the main numerical methods currently used for modelling wave energy converter arrays are described in this paper, using a common framework. This allows a qualitative comparative analysis of the different methods to be performed at the end of the paper. This includes consideration of the conditions under which the models may be applied, the output of the models and the relationship between array size and computational effort. Guidance for developers is also presented on the most suitable numerical method to use for given aspects of WEC farm design. For instance, certain models are more suitable for studying near-field effects, whilst others are preferable for investigating far-field effects of the WEC farms. Furthermore, the analysis presented in this paper identifies areas in which the numerical modelling of WEC arrays is relatively weak and thus highlights those in which future developments are required.

Simulation of Oil Shale Semi-Coke Particle Cold Transportation in a Spouted Bed Using CPFD Method

The spouted bed was widely used due to its good mixing of particles and effective phase transferability between the gas and solid phase. In this paper, the transportation process of particles in a 3D spouted bed was studied using the Computational Particle Fluid Dynamics (CPFD) numerical method. Experiments were conducted to verify the validity of the simulation results. Distributions of the pressure, velocities and particle concentration of transportation devices were investigated. The motion state and characteristics of multiphase flows in the transportation device were demonstrated under various operating conditions. The results showed that a good consistency was obtained between the simulated results and the experimental results. The motion characteristics of the gas-solid two-phase flow in the device was effectively predicted, which could assist the optimal operating condition estimation for the spouted transportation process.

Geometrical dependence in photoionization of H-2(+) in high-intensity, high-frequency, ultrashort laser pulses

The ionization dynamics of H2 + exposed to high-intensity, high-frequency, ultrashort laser pulses is investigated with two theoretical approaches. The time-dependent Schrödinger equation is solved by a direct numerical method, and a simple two-center interference-diffraction model is studied. The energy and angular distributions of the photoelectron for various internuclear distances and relative orientations between the internuclear axis of the molecule and the polarization of the field are calculated. The main features of the photoelectron spectrum pattern are described well by the interference-diffraction model, and excellent quantitative agreement between the two methods is found. The effect of quantal vibration on the photoelectron spectrum is also calculated. We find that vibrational average produces some broadening of the main features, but that the patterns remain clearly distinguishable.

Doppler cooling of gallium atoms: 2. Simulation in complex multilevel systems

This paper derives a general procedure for the numerical solution of the Lindblad equations that govern the coherences arising from multicoloured light interacting with a multilevel system. A systematic approach to finding the conservative and dissipative terms is derived and applied to the laser cooling of p-block elements. An improved numerical method is developed to solve the time-dependent master equation and results are presented for transient cooling processes. The method is significantly more robust, efficient and accurate than the standard method and can be applied to a broad range of atomic and molecular systems. Radiation pressure forces and the formation of dynamic dark states are studied in the gallium isotope 66Ga.

Two-dimensional hygrothermal transfer in porous building materials

A two-dimensional mathematical model for evaluating the simultaneous heat and moisture migration in porous building materials was proposed. Vapor content and temperature were chosen as the principal driving potentials. The numerical solution was based on the control volume finite difference technique with fully implicit scheme in time. Two validation experiments were developed in this study. The evolution of transient moisture distributions in both one-dimensional and two-dimensional cases was measured. A comparison between experimental results and those obtained by the numerical model proves that they are fully consistent with each other. The model can be easily integrated into a whole building heat, air and moisture transfer model. Another main advantage of the present numerical method lies in the fact that the required moisture transport properties are comparatively simple and easy to determine.

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

Design and testing of transonic linear cascade tunnel with optimized slotted walls

In linear cascade wind tunnel tests, a high level of pitchwise periodicity is desirable to reproduce the azimuthal periodicity in the stage of an axial compressor or turbine. Transonic tests in a cascade wind tunnel with open jet boundaries have been shown to suffer from spurious waves, reflected at the jet boundary, that compromise the flow periodicity in pitch. This problem can be tackled by placing at this boundary a slotted tailboard with a specific wall void ratio s and pitch angle a. The optimal value of the s-a pair depends on the test section geometry and on the tunnel running conditions. An inviscid two-dimensional numerical method has been developed to predict transonic linear cascade flows, with and without a tailboard, and quantify the nonperiodicity in the discharge. This method includes a new computational boundary condition to model the effects of the tailboard slots on the cascade interior flow. This method has been applied to a six-blade turbine nozzle cascade, transonically tested at the University of Leicester. The numerical results identified a specific slotted tailboard geometry, able to minimize the spurious reflected waves and regain some pitchwise flow periodicity. The wind tunnel open jet test section was redesigned accordingly. Pressure measurements at the cascade outlet and synchronous spark schlieren visualization of the test section, with and without the optimized slotted tailboard, have confirmed the gain in pitchwise periodicity predicted by the numerical model. Copyright © 2006 by ASME.

Generalized Langevin equation: An efficient approach to nonequilibrium molecular dynamics of open systems

The generalized Langevin equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general nonequilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e. g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly nontrivial, especially if the memory kernels decay nonexponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this paper, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.

Hydrodynamic Loading on a Bottom Hinged Oscillating Wave Surge Converter

This paper presents experimental and numerical studies into the hydrodynamic loading of a bottom-hinged large buoyant flap held rigidly upright in waves. Possible applications and limitations of physical experiments, a linear potential analytical method, a linear potential numerical method, a weakly non-linear tool and RANS CFD simulations are discussed. Different domains of applicability of these research techniques are highlighted considering the validity of underlying assumptions, complexity of application and feasibility in terms of resources like time and computing power needed to obtain results. Conclusions are drawn regarding the future extension of the numerical methods to the case of a moving flap.

Relativistic breather-type solitary waves with linear polarization in cold plasmas

Linearly polarized solitary waves, arising from the interaction of an intense laser pulse with a plasma, are investigated. Localized structures, in the form of exact numerical nonlinear solutions of the one-dimensional Maxwell-fluid model for a cold plasma with fixed ions, are presented. Unlike stationary circularly polarized solitary waves, the linear polarization gives rise to a breather-type behavior and a periodic exchange of electromagnetic energy and electron kinetic energy at twice the frequency of the wave. A numerical method based on a finite-differences scheme allows us to compute a branch of solutions within the frequency range Ωmin<Ω<ωpe, where ωpe and Ωmin are the electron plasma frequency and the frequency value for which the plasma density vanishes locally, respectively. A detailed description of the spatiotemporal structure of the waves and their main properties as a function of Ω is presented. Small-amplitude oscillations appearing in the tail of the solitary waves, a consequence of the linear polarization and harmonic excitation, are explained with the aid of the Akhiezer-Polovin system. Direct numerical simulations of the Maxwell-fluid model show that these solitary waves propagate without change for a long time.