Random field sampling for a simplified model of melt-blowing considering turbulent velocity fluctuations

In melt-blowing very thin liquid fiber jets are spun due to high-velocity air streams. In literature there is a clear, unsolved discrepancy between the measured and computed jet attenuation. In this paper we will verify numerically that the turbulent velocity fluctuations causing a random aerodynamic drag on the fiber jets -- that has been neglected so far -- are the crucial effect to close this gap. For this purpose, we model the velocity fluctuations as vector Gaussian random fields on top of a k-epsilon turbulence description and develop an efficient sampling procedure. Taking advantage of the special covariance structure the effort of the sampling is linear in the discretization and makes the realization possible

A Dynamic Stochastic Model for DNA Replication Initiation in Early Embryos

Background: Eukaryotic cells seem unable to monitor replication completion during normal S phase, yet must ensure a reliable replication completion time. This is an acute problem in early Xenopus embryos since DNA replication origins are located and activated stochastically, leading to the random completion problem. DNA combing, kinetic modelling and other studies using Xenopus egg extracts have suggested that potential origins are much more abundant than actual initiation events and that the time-dependent rate of initiation, I(t), markedly increases through S phase to ensure the rapid completion of unreplicated gaps and a narrow distribution of completion times. However, the molecular mechanism that underlies this increase has remained obscure.Methodology/Principal Findings: Using both previous and novel DNA combing data we have confirmed that I(t) increases through S phase but have also established that it progressively decreases before the end of S phase. To explore plausible biochemical scenarios that might explain these features, we have performed comparisons between numerical simulations and DNA combing data. Several simple models were tested: i) recycling of a limiting replication fork component from completed replicons; ii) time-dependent increase in origin efficiency; iii) time-dependent increase in availability of an initially limiting factor, e. g. by nuclear import. None of these potential mechanisms could on its own account for the data. We propose a model that combines time-dependent changes in availability of a replication factor and a fork-density dependent affinity of this factor for potential origins. This novel model quantitatively and robustly accounted for the observed changes in initiation rate and fork density.Conclusions/Significance: This work provides a refined temporal profile of replication initiation rates and a robust, dynamic model that quantitatively explains replication origin usage during early embryonic S phase. These results have significant implications for the organisation of replication origins in higher eukaryotes

Modeling Inhomogeneous DNA Replication Kinetics

In eukaryotic organisms, DNA replication is initiated at a series of chromosomal locations called origins, where replication forks are assembled proceeding bidirectionally to replicate the genome. The distribution and firing rate of these origins, in conjunction with the velocity at which forks progress, dictate the program of the replication process. Previous attempts at modeling DNA replication in eukaryotes have focused on cases where the firing rate and the velocity of replication forks are homogeneous, or uniform, across the genome. However, it is now known that there are large variations in origin activity along the genome and variations in fork velocities can also take place. Here, we generalize previous approaches to modeling replication, to allow for arbitrary spatial variation of initiation rates and fork velocities. We derive rate equations for left- and right-moving forks and for replication probability over time that can be solved numerically to obtain the mean-field replication program. This method accurately reproduces the results of DNA replication simulation. We also successfully adapted our approach to the inverse problem of fitting measurements of DNA replication performed on single DNA molecules. Since such measurements are performed on specified portion of the genome, the examined DNA molecules may be replicated by forks that originate either within the studied molecule or outside of it. This problem was solved by using an effective flux of incoming replication forks at the model boundaries to represent the origin activity outside the studied region. Using this approach, we show that reliable inferences can be made about the replication of specific portions of the genome even if the amount of data that can be obtained from single-molecule experiments is generally limited

Computational Methods to Study Kinetics of DNA Replication

New technologies such as DNA combing have led to the availability of large quanti-ties of data that describe the state of DNA while undergoing replication in S phase. In this chapter, we describe methods used to extract various parameters of replica-tion — fork velocity, origin initiation rate, fork density, numbers of potential and utilized origins — from such data. We first present a version of the technique that applies to “ideal ” data. We then show how to deal with a number of real-world complications, such as the asynchrony of starting times of a population of cells, the finite length of fragments used in the analysis, and the finite amount of DNA in a chromosome. Key words: DNA replication, replication fork velocity, origin initiation

Evidence for Sequential and Increasing Activation of Replication Origins along Replication Timing Gradients in the Human Genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics

On adjoint-based optimization of a free surface Stokes flow

This work deals with the optimal control of a free surface Stokes flow which responds to an applied outer pressure. Typical applications are fiber spinning or thin film manufacturing. We present and discuss two adjoint-based optimization approaches that differ in the treatment of the free boundary as either state or control variable. In both cases the free boundary is modeled as the graph of a function. The PDE-constrained optimization problems are numerically solved by the BFGS method, where the gradient of the reduced cost function is expressed in terms of adjoint variables. Numerical results for both strategies are finally compared with respect to accuracy and efficiency

Slender Body Theory for the Dynamics of Curved Viscous Fibers

The paper at hand presents a slender body theory for the dynamics of a curved inertial viscous Newtonian ber. Neglecting surface tension and temperature dependence, the ber ow is modeled as a three-dimensional free boundary value problem via instationary incompressible Navier-Stokes equations. From regular asymptotic expansions in powers of the slenderness parameter leading-order balance laws for mass (cross-section) and momentum are derived that combine the unrestricted motion of the ber center-line with the inner viscous transport. The physically reasonable form of the one-dimensional ber model results thereby from the introduction of the intrinsic velocity that characterizes the convective terms.The paper at hand presents a slender body theory for the dynamics of a curved inertial viscous Newtonian ber. Neglecting surface tension and temperature dependence, the ber ow is modeled as a three-dimensional free boundary value problem via instationary incompressible Navier-Stokes equations. From regular asymptotic expansions in powers of the slenderness parameter leading-order balance laws for mass (cross-section) and momentum are derived that combine the unrestricted motion of the ber center-line with the inner viscous transport. The physically reasonable form of the one-dimensional ber model results thereby from the introduction of the intrinsic velocity that characterizes the convective terms

Asymptotic modeling framework for fiber-flow interactions in a two-way coupling

In this work we describe fiber-flow interactions by help of a two-way coupling approach that is based on slender-body theory and the modeling of exchange functions in terms of drag forces and heat sources. The exchange functions are incorporated in the conservation equations for linear momentum and energy with respect to flow and fibers and satisfy a generalized action-reaction principle