Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

激光聚焦爆炸耦合场的理论建模与跨尺度计算

针对激光聚焦爆炸的电磁-热力耦合效应,在宏观尺度上,把描述激光电磁波散射和传播的Maxwell方程和高温高压气动流场的Euler方程结合起米,利用热力学状态方程(EOS)和电离平衡方程(Saha方程)并通过理论建模和数值仿真,研究和揭示激光聚焦爆炸效应及激光支持吸收波(LSC/LSD)的产生和演化、以及相关的反冲压力和动量耦合等相互作用机制.

Self-organized generation of transverse waves in diverging cylindrical detonations

Self-organized generation of transverse waves associated with the transverse wave instabilities at a diverging cylindrical detonation front was numerically studied by solving two-dimensional Euler equations implemented with an improved two-step chemical kinetic model. After solution validation, four mechanisms of the transverse wave generation were identified from numerical simulations, and referred to as the concave front focusing, the kinked front evolution, the wrinkled front evolution and the transverse wave merging, respectively. The propagation of the cylindrical detonation is maintained by the growth of the transverse waves that match the rate of increase in surface area of the detonation front to asymptotically approach a constant average number of transverse waves per unit length along the circumference of the detonation front. This cell bifurcation phenomenon of cellular detonations is discussed in detail to gain better understanding on detonation physics.

柱面气相散心爆轰波胞格演化规律研究

本文采用二阶精度NND格式，应用改进的二步法爆轰计算模拟，通过求解二维Euler方程对柱面气相散心爆轰波胞格演化过程进行了数值模拟。计算结果表明在传播过程中，空间尺度的扩张导致了散心爆轰波后气流的自然膨胀，使得多波结构的爆轰阵面呈现出显著的胞格自组织特性。根据计算结果与理论分析，本文归纳了五种胞格演化模式，分别命名为内凹波阵面会聚、波阵面扭结、褶皱波面失稳、胞格自合并和三波点滑移，并定义了各种模式的物理特征，分析了其相关的演化机制和规律。

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2°,飞行高度20 km的巡航状态下的静气动弹性特性进行了数值模拟.验证了两种求解器对静气动弹性模拟的准确性.同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟.另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.

NUMERICAL AND EXPERIMENTAL STUDY ON LIQUID-SOLID FLOW IN A HYDROCYCLONE

Hydrocyclones are widely used in industry, of which the geometrical design using CFD techniques is gaining more popularity in recent years. In this study, the Euler-Euler approach and the Reynolds stress model are applied to simulate the liquid-solid flowfield in a hydrocyclone. The methodology is validated by a good agreement between experimental data and numerical results. Within the research range, the simulation indicates that the liquid-solid separation mainly occurs in the conical segment, and increasing conical height or decreasing cylindrical height helps to improve the grade efficiencies of solid particles. Based on these results, two of the same hydrocyclones are designed and installed in series to establish a liquid-solid separation system. Many experiments are then conducted under different conditions, in which the effects of the water cut and the second hydrocyclone on the separation are investigated. The results also confirm that smaller solid particles are more susceptible to the inlet conditions, and the second hydrocyclone plays a more important role as the water cut reduces.

Part I. On the breaking of nonlinear dispersive waves. Part II. Variational principles in continuum mechanics

A model equation for water waves has been suggested by Whitham to study, qualitatively at least, the different kinds of breaking. This is an integro-differential equation which combines a typical nonlinear convection term with an integral for the dispersive effects and is of independent mathematical interest. For an approximate kernel of the form e^(-b|x|) it is shown first that solitary waves have a maximum height with sharp crests and secondly that waves which are sufficiently asymmetric break into "bores." The second part applies to a wide class of bounded kernels, but the kernel giving the correct dispersion effects of water waves has a square root singularity and the present argument does not go through. Nevertheless the possibility of the two kinds of breaking in such integro-differential equations is demonstrated.

Difficulties arise in finding variational principles for continuum mechanics problems in the Eulerian (field) description. The reason is found to be that continuum equations in the original field variables lack a mathematical "self-adjointness" property which is necessary for Euler equations. This is a feature of the Eulerian description and occurs in non-dissipative problems which have variational principles for their Lagrangian description. To overcome this difficulty a "potential representation" approach is used which consists of transforming to new (Eulerian) variables whose equations are self-adjoint. The transformations to the velocity potential or stream function in fluids or the scaler and vector potentials in electromagnetism often lead to variational principles in this way. As yet no general procedure is available for finding suitable transformations. Existing variational principles for the inviscid fluid equations in the Eulerian description are reviewed and some ideas on the form of the appropriate transformations and Lagrangians for fluid problems are obtained. These ideas are developed in a series of examples which include finding variational principles for Rossby waves and for the internal waves of a stratified fluid.

Photon-photon scattering in a plasma channel

The Heisenberg-Euler correction due to photon-photon scattering, a still unverified quantum electrodynamics effect, on electromagnetic wave interaction inside a plasma channel is investigated theoretically. From a signal laser beam in the relativistic overdense plasma channel, photon-photon scattering can produce a detectable output beam of different frequency and polarization. (C) 2003 American Institute of Physics.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Uncovering the Lagrangian from observations of trajectories

We approach the problem of automatically modeling a mechanical system from data about its dynamics, using a method motivated by variational integrators. We write the discrete Lagrangian as a quadratic polynomial with varying coefficients, and then use the discrete Euler-Lagrange equations to numerically solve for the values of these coefficients near the data points. This method correctly modeled the Lagrangian of a simple harmonic oscillator and a simple pendulum, even with significant measurement noise added to the trajectories.

Transmission matrices and lumped parameter models for continuous systems

The use of transmission matrices and lumped parameter models for describing continuous systems is the subject of this study. Non-uniform continuous systems which play important roles in practical vibration problems, e.g., torsional oscillations in bars, transverse bending vibrations of beams, etc., are of primary importance.

A new approach for deriving closed form transmission matrices is applied to several classes of non-uniform continuous segments of one dimensional and beam systems. A power series expansion method is presented for determining approximate transmission matrices of any order for segments of non-uniform systems whose solutions cannot be found in closed form. This direct series method is shown to give results comparable to those of the improved lumped parameter models for one dimensional systems.

Four types of lumped parameter models are evaluated on the basis of the uniform continuous one dimensional system by comparing the behavior of the frequency root errors. The lumped parameter models which are based upon a close fit to the low frequency approximation of the exact transmission matrix, at the segment level, are shown to be superior. On this basis an improved lumped parameter model is recommended for approximating non-uniform segments. This new model is compared to a uniform segment approximation and error curves are presented for systems whose areas very quadratically and linearly. The effect of varying segment lengths is investigated for one dimensional systems and results indicate very little improvement in comparison to the use of equal length segments. For purposes of completeness, a brief summary of various lumped parameter models and other techniques which have previously been used to approximate the uniform Bernoulli-Euler beam is a given.

Advancements in jet turbulence and noise modeling: accurate one-way solutions and empirical evaluation of the nonlinear forcing of wavepackets

Jet noise reduction is an important goal within both commercial and military aviation. Although large-scale numerical simulations are now able to simultaneously compute turbulent jets and their radiated sound, lost-cost, physically-motivated models are needed to guide noise-reduction efforts. A particularly promising modeling approach centers around certain large-scale coherent structures, called wavepackets, that are observed in jets and their radiated sound. The typical approach to modeling wavepackets is to approximate them as linear modal solutions of the Euler or Navier-Stokes equations linearized about the long-time mean of the turbulent flow field. The near-field wavepackets obtained from these models show compelling agreement with those educed from experimental and simulation data for both subsonic and supersonic jets, but the acoustic radiation is severely under-predicted in the subsonic case. This thesis contributes to two aspects of these models. First, two new solution methods are developed that can be used to efficiently compute wavepackets and their acoustic radiation, reducing the computational cost of the model by more than an order of magnitude. The new techniques are spatial integration methods and constitute a well-posed, convergent alternative to the frequently used parabolized stability equations. Using concepts related to well-posed boundary conditions, the methods are formulated for general hyperbolic equations and thus have potential applications in many fields of physics and engineering. Second, the nonlinear and stochastic forcing of wavepackets is investigated with the goal of identifying and characterizing the missing dynamics responsible for the under-prediction of acoustic radiation by linear wavepacket models for subsonic jets. Specifically, we use ensembles of large-eddy-simulation flow and force data along with two data decomposition techniques to educe the actual nonlinear forcing experienced by wavepackets in a Mach 0.9 turbulent jet. Modes with high energy are extracted using proper orthogonal decomposition, while high gain modes are identified using a novel technique called empirical resolvent-mode decomposition. In contrast to the flow and acoustic fields, the forcing field is characterized by a lack of energetic coherent structures. Furthermore, the structures that do exist are largely uncorrelated with the acoustic field. Instead, the forces that most efficiently excite an acoustic response appear to take the form of random turbulent fluctuations, implying that direct feedback from nonlinear interactions amongst wavepackets is not an essential noise source mechanism. This suggests that the essential ingredients of sound generation in high Reynolds number jets are contained within the linearized Navier-Stokes operator rather than in the nonlinear forcing terms, a conclusion that has important implications for jet noise modeling.

基于双面反射镜的N×N光开关的特性分析和实验研究

报道了基于双面反射镜的N×N光开关器件。介绍了使用双面反射镜的2×2, 4×4光开关的集成光路设计和工作原理; 采用Benes网络, 以2×2和4×4光开关为基本单元的N×N光开关器件的整体结构, 并根据“一笔画”原理, 分析了4×4, 8×8和16×16光开关矩阵的可重排无阻塞特性和光开关矩阵的光路选择算法。最后, 基于2×2, 4×4光开关技术制备了16×16光开关矩阵。测试表明, 该器件具有良好的插入损耗、回波损耗、串扰和开关时间等性能, 从而验证了设计思想和工艺的可行性。在基于双面反射镜的光开关矩

Reavaliação deposicional da Bacia de São José de Itaboraí com base em dados geológicos e geofísicos

A Bacia de São José de Itaboraí está localizada no Município de Itaboraí, no Estado do Rio de Janeiro. Ela foi descoberta em 1928, pelo Engenheiro Carlos Euler, que após analisar um suposto caulim encontrado na Fazenda São José pelo seu então proprietário, Sr. Ernesto Coube, verificou que se tratava de calcário. Os Professores Rui Lima e Silva e Othon H. Leonardos, enviados ao local para estudos, encontraram uma grande quantidade de fósseis de gastrópodes continentais, despertando o interesse científico pela região. Os estudos preliminares de campo e análises químicas evidenciaram boas perspectivas de exploração do calcário para a fabricação de cimento do tipo Portland. Por mais de 50 anos, a Companhia Nacional de Cimento Portland Mauá (CNCPM) explorou a pedreira. Desde sua descoberta, a Bacia de São José, paralelamente às atividades de mineração, foi objeto de pesquisas científicas realizadas por geólogos, paleontólogos e arqueólogos. No início da década de 80, a Cia. de Cimento Mauá decidiu abandonar a área em função do esgotamento econômico da reserva de minério. Com a retirada das bombas que impediam a inundação da pedreira, formou-se uma lagoa que passou a impedir o livre acesso aos afloramentos. Desde então as pesquisas sobre a Bacia ficaram concentradas aos materiais coletados no período de exploração de calcário. Material esse distribuído no Museu Nacional (MN), Departamento Nacional da Produção Mineral (DNPM), Instituto de Geociências da UFRJ, entre outros. Em 1990, a área que pertencia a CNCPM foi desapropriada por pressão da comunidade científica. A mesma passou a pertencer ao Município de Itaboraí, que criou o Parque Paleontológico de São José de Itaboraí, por meio da Lei 1.346, de 12 de dezembro de 1995. O objetivo desse trabalho foi gerar novos dados através do método geofísico conhecido como magnetometria. Para isso foram realizados levantamentos de campo utilizando um magnetômetro portátil e GPS, foram analisados e corrigidos dados utilizando softwares específicos, elaborados modelos e criados perfis a partir de descrições de testemunhos de sondagem. Os resultados obtidos visam possibilitar uma nova interpretação da geologia e da estratigrafia da bacia, dando condições para que se possa ter uma atualização dos conhecimentos relacionados à região, após quase meio século de atividade mineradora.